NASA  CONTRACTOR 
REPORT 


SPECTROSCOPIC  MEASUREMENTS 

OF  WATER  VAPOR  PLASMAS 

AT  HIGH  RESOLUTION  —  THE  OPTICAL 

TRANSITION  PROBABILITIES 

FOR  OH(A  -  X  11)  , 

i  .  .  . . . 

1  T.-: 

by  Lothar  Klein 


Prepared  by 

GREYRAD  CORPORATION 

Princeton,  N.J. 

for  Lewis  Research  Center 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION  •  WASHINGTON,  D.  C.  •  AUGUST  1972 


THIS  DOCUMENT  IS  BEST 
QUALITY  AVAILABLE.  THE 
COPY  FURNISHED  TO  DTIC 
CONTAINED  A  SIGNIFICANT 
NUMBER  OF  PAGES  WHICH  DO 
NOT  REPRODUCE  LEGIBLY. 


1,  Report  No.  2.  Government  Accession  No, 

CR-2069 

3.  Recipient's  Catalog  No. 

4.  Title  and  Subtitle  SPECTROSCOPIC  MEASUREMENTS  OF  WATER 
VAPOR  PLASM;\S  AT  HIGH  RESOLUTION  -  THE  OPTICAL 
TRANSITION  PROBABILITIES  FOR  OH(A^S  -  X^II) 

5.  Report  Date 

August  1972 

6.  Performing  Organization  Code 

7.  Author(s) 

Lothar  Klein 

8.  Performing  Organization  Report  No. 

None 

10.  Work  Unit  No. 

9.  Performing  Organization  Name  and  Address 

Greyrad  Corporation 

Princeton,  New  Jersey 

11.  Contract  or  Grant  No. 

NAS  3-13424 

13.  Type  of  Report  and  Period  Covered 

Contractor  Report 

12.  Sponsoring  Agency  Name  and  Address 

National  Aeronautics  and  Space  Administration 

Washington,  D.C.  20546 

14.  Sponsoring  Agency  Code 

15.  Supplementary  Notes 

Project  Manager,  Erwin  A.  Lezberg,  Physics  and  Chemistry  Division,  NASA  Lewis  Research 
Center,  Cleveland,  Ohio 


16.  Abstract 

Emission  and  absorption  spectra  of  water-vapor  plasmas  generated  in  a  wall -stabilized  arc  at 
atmospheric  pressure  and  4  A  current,  and  at  0.  03  atm  and  15  to  50  A,  were  measured  at  high 
spatial  and  spectral  resolution.  The  gas  temperature  was  determined  from  the  shape  of  Doppler- 
broadened  rotational  lines  of  OH  since  the  measured  resolving  power  of  the  spectrometer  was 
400,  000.  The  observed  non- equilibrium  population  distributions  over  the  energy  levels  of  atoms 
are  interpreted  in  terms  of  a  theoretical  state  model  for  diffusion- controlled  arc  plasmas.  Ex¬ 
cellent  correlation  is  achieved  between  the  measured  and  the  predicted  occupation  of  the  hydrogen 
energy  levels.  It  is  shown  that  the  population  distribution  over  the  non- predissociating 
rotational- vibrational  levels  of  the  A  S  state  of  OH  is  close  to  an  equilibrium  distribution  at  the 
gas  temperature,  although  the  total  density  of  this  state  is  much  higher  than  its  equilibrium 
density.  The  reduced  intensities  of  the  rotational  lines  originating  in  these  levels  yielded 
Boltzmann  plots  that  were  strictly  linear,  but  the  slope  did  not  correspond  to  the  measured  gas 
temperature.  This  is  attributed  to  an  exponential  dependence  of  the  transition  moment  on  the 
rotational  energy  term.  Appropriate  correction  factors  were  calculated  and  compared  with 
theoretical  predictions. 


17.  Key  Words  (Suggested  by  Author(s} )  18.  Distribution  Statement 

Hydroxyl  ultraviolet  bands;  Transition  prob-  Unclassified  -  unlimited 

abilities;  Line  strengths;  Spectroscopic  tem¬ 
perature  measurements;  Non- equilibrium  arc 

plasmas;  Plasma  spectroscopy  _ 


19.  Security  Classif.  (of  this  report} 

20.  Security  Classif.  (of  this  page) 

21.  No.  of  Pages 

22.  Price* 

Unclassified 

Unclassified 

150 

$3.00 

For  sale  by  the  National  Technical  information  Service,  Springfield,  Virginia  22151 


DSC  QUAIJTY  X] 


FOREWORD 


The  material  presented  in  this  report  was  included  in  the  doctoral 
dissertation  submitted  by  the  author  to  the  National  University  of  Buenos 
Aires,  His  present  address  is  the  Hasma  Spectrum  Corporation,  28?  Schenck 
Avenue,  Great  Neck,  New  York  11021. 

The  high-resolution  plasma  spectrometer  was  built  in  collaboration 
with  H,  J.  Babrov,  now  at  Indiana  University  at  South  Bend,  while  he  and 
the  avithor  were  with  the  Control  Instrument  Division  of  the  Warner  &  Swasey 
Company.  The  grating  drive  was  designed  by  W.  L,  Ronai  amd  S.  A,  Dolin, 
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SUMMARY 


Hiech-order  spectroscopy  has  been  applied  to  quantitative  measurements  of 
the  spectra  of  waier-vapor  plasmas  generated  in  a  watll -stabilized  arc.  Because 
of  the  high  spectreil  and  spatial  resolution  achieved  in  end-on  measurements,  the 
emission  and  absorption  coefficients  pertaining  to  homogeneous  arc  regions  were 
obtained  directly  from  measured  line  spectra,  leading  to  the  occupation  of  the 
upper  aind  the  lower  state  for  the  transition.  The  electron  density  was  derived 
from  the  profile  of  the  Stark -broadened  Balmer  lines.  The  gas  temperature  was 
determined  from  the  half-width  of  a  Doppler-broadened  rotational  line  of  OH. 

The  measured  resolving  power  of  the  spectrometer  was  of  the  order  of  ^0,000  in 
these  measurements. 

Significant  departures  from  local  thermodynamic  equilibrium  (LTE)  were  ob¬ 
served  in  a  water-vapor  arc  at  reduced  pressure,  and  ailso  in  the  peripheral  zone 
of  a  low -current  arc  at  atmospheric  pressure  with  water  vapor  in  an  argon  matrix. 

The  "distribution  temperature”  characterizing  the  population  distribution  over  highly 
excited  energy  levels  of  atoms  is  not  only  much  lower  than  the  electron  temperature^ 
but  also  substantially  lower  than  the  gas  temperature.  The  noted  departures  from 
LTE  occur  because  of  the  lack  of  equipartition  of  kinetic  energy  between  the 
electrons  and  the  heavy  particles  in  the  plasma,  and  because  of  diffusion  of  ground 
state  neutrals  from  the  arc  peripheiy  into  the  arc  center  and  ambipolar  diffusion 
of  electrons  and  ions  in  the  opposite  direction.  Based  on  these  mechanisms,  the 
various  interaction  processes  in  a  hydrogen  plasma  are  analyzed  quantitatively 
on  the  microscopic  level,  and  it  is  shown  that  the  occupation  of  the  excited  ener¬ 
gy  states  is  predominantly  determined  by  step-wise  excitation  induced  by  electron 
collision.  Excellent  agreement  is  obtained  between  the  predicted  and  the  measured 
number  densities  of  the  excited  states  of  hydrogen  in  the  water-vapor  plasma  at 
reduced  pressure. 

Ionization  of  atoms  and  dissociation  of  molecules  such  as  OH  occur  by  similar 
processes  in  the  plasma,  although  the  latter  is  dominated  by  heavy  particle  colli¬ 
sion,  and  the  relevant  temperature  is  then  the  gas  temperature.  A  close  analogy 
exists  between  the  population  distribution  over  the  auto -ionizing  levels  of  oxygen 

and  the  predissociating  levels  of  OH,  The  lower  rotational -vibrational  levels 
2 

of  the  A  E  excited  electronic  state  of  OH  are  effectively  uncoupled  from  the  disso¬ 
ciated  state  because  of  the  step-wise  excitation  mechanism  in  conjunction  with 
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the  near^-eaiiality  of  the  transition  rates.  The  relative  population  distribution 

over  these  levels  is  therefore  close  to  an  equilibrium  distribution  at  the  ^as 

2 

temperature^  even  i^hen  the  total  density  of  the  A  E  state  is  orders  of  magnitude 
higher  than  its  equilibrium  density. 

The  reduced  intensities  of  the  non-predissociating^ rotational  lines  belonging 

2  2 

to  the  various  bands  of  the  OH  (A  E  X  H  )  system  yielded  Boltzmann  plots  that 
were  strictly  linear,  although  the  slope  did  not  correspond  to  the  measiired  gas 
temperature.  This  is  attributed  to  an  exponential  dependence  of  the  transition 
moment  on  the  rotational  energy  term«  Appropiate  correction  factors  were  calcu- 
lated  and  compared  with  theoretical  predictions. 
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1.  INTRODUCTION 


The  use  of  quantitative  spectroscoi^  as  a  diagnostic  tool  in  nany  areas 
of  current  interest  has  been  obstructed  by  the  paucity  of  accurate  data  on 
transition  probabilities  especially  for  molecular  lines  and  by  the  lack  of 
appropiate  state  models  to  characterise  plasmas  in  absence  of  local  thermodynamic 
equilibrium  (LTE).  The  plasmas  obtained  in  stabilised  arcs  at  atmospheric 
pressure  have  been  shown  to  be  in  LTE  (l),  and  have  been  used  in  the  past  to 
determine  transition  probabilities  for  atomic  lines «  The  initial  concern  of 
this  investigation  was  to  apply  the  same  approach  to  molecular  species  that 
axe  still  stable  at  arc  temperatures,  and  a  water  vapor  plasma  generated  in  a 
wall-stabilised  arc  was  used  to  determine  the  transition  probabilities  of  ro¬ 
tational  lines  belonging  to  the  ultraviolet  bands  of  OH^ 

In  LTE  detailed  baletncing  exists  between  the  collisionail  processes  in  the 
plasma  and  the  pertinent  rates  are  much  larger  than  for  the  equivalent  radiative 
processes.  Although  the  latter  are  not  in  detailed  balance  (i,e.  photons  can 
escape  the  plasma),  the  species  concentration  and  the  population  distributions 
over  energy  levels  in  the  plasma  are  then  practically  the  saune  as  in  a  system 
in  complete  thermodynamic  equilibrium  at  the  same  temperature,  pressure,  and 
chemical  composition  at  standard  conditions*  The  occupation  of  energy  levels  can 
thus  be  calculated  which  makes  LTE  plasmas  a  convenient  source  for  the  experimental 
determination  of  transition  probabilities,  since  the  number  density  of  emitters 
must  be  known  in  order  to  derive  the  opticad  transition  probability  from  a 
spectroscopic  measurement  of  the  integrated  line  radiance* 
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In  Section  2G  Me  have  calculated  the  equilibrium  composition  of  a  Kater 
vapor  plasma  at  various  presstires  as  a  function  of  temperature „  Provided  the 
plasma  is  in  LTE5,  it  is  then  sufficient  to  determine  experimentally  the  temper¬ 
ature  of  an  isothermal  plasma  region  in  order  to  infer  the  transition  probabilities 
from  the  istensities  of  the  emitted  optic ally- thin  lines «  In  principle ^  a  suitable 
"thermometer  line"  of  knoMn  transition  probability  can  b©  used  to  determine  the 
plasma  temperatursp  HoMever^  as  argued  in  Section  2D  (see  also  Bef,  l),  the 
generalis^ed  Fowler-Milne  method  is  preferable  because  over  an  extended  temper¬ 
ature  range  a  one-to-one  correspondence  is  established  betHeea  the  measured 
line  intensity  and  the  calculated  occupation  of  the  upper  energy  level  for  the 
transition  0 

Because  of  the  quantization  of  molecular  energy  states  corresponding  to 
nuclear  motion p  transitions  between  two  electronic  states  of  a  molecule  fora  an 
entire  band  system «  Rigorously^  a  transition  probability  or  a  related  quantity 

-  oscillator  strength  or  transition  moment  ^  -  can  only  he  assigned  to  an  indi¬ 
vidual  rotational  line.  It  is  still  useful,  however,  to  consider  the  relation¬ 
ships  between  the  transition  moments  for  the  rotational  lines  in  a  band,  and 
the  various  bands  in  a  systsn, which  can  be  derived  from  theoretical  assumptions 
and  lead  essentially  to  a  factoring  of  the  transition  moment.  The  approach 
followed  here  was  to  determine  first  the  transition  probabilities  for  the  indi¬ 
vidual  rotational  lines  and  then  to  use  the  concept  of  Franck-Condois  factors 
("band  stengths")  and  Hoenl -London  factors  ("line  strengths")  as  a  first-order 
approximation  for  the  representation  of  the  measured  data  in  terms  of  the  transition 
moments  by  introducing  appr opiate  correction  factors. 

*)  square  of  the  matrix  element  for  electric  dipole  transition 
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A  repiresentation  of  rotational  line  intensities  in  terns  of  "band  transition 
probabilities"  or  "band  oscillator  stirengths"  does  not  appear  to  be  useful  since 
a  single  frequency  cannot  be  assigned  to  a  beuid.  The  introduction  of  an  average 
band  frequency  laplies  a  weighting  of  the  vairious  rotational  lines  which  leads 
to  different  values  depending  on  the  particulau:  population  distribution  over  the 
rotational  energy  levels. 

Almost  all  previous  measuarements  of  the  spectra  of  stabilised  ares  were 
made  by  viewing  the  arc  side-on.  In  order  to  obtain  the  spectral  data  pertaining 
to  Isothexnal  axe  regions  (radial  zones  for  ares  of  eylindzleal  symnetzy), 
spatial  scanning  of  the  arc  and  an  Abel  inversion  of  the  experinental  lateral 
intensity  profile  are  required.  This  is  a  viable  approach  when  the  integrated 
intensities  of  optlcally-thin  lines  are  measured  (see  Section  4).  Severe 
problems  arise,  however,  in  the  interpretation  of  line  absorptanees  and  line 
shapes  measured  across  spatial  t^perature  gradients. 

These  difficulties  are  avoided  by  viewing  the  are  in  the  axi€d.  direction 
(end-on  mode),  since  then  an  essentially  homogeneous  plasma  sample  is  seen  by 
the  spectrometer,  provided  the  arc  axis  is  accurately  aligned  parallel  to  the 
optical  axis  of  the  system,  and  approplate  stops  are  used  to  limit  the  field  ^ 

of  view  and  the  aperture  of  the  spectrometer.  Even  when  extreme  care  has  been 
exercised  in  setting-up  the  arc,  a  direct  proof  of  correct  arc  alignment  is  best 
obtained  by  an  actual  measurement.  In  Section  JD  we  describe  such  an  experiment 
irtilch  edso  demonstrates  the  relative  esse  with  which  information  on  line  shapes 
can  be  obtained  by  high-resolution  measurements  in  the  end-on  node. 

A  scanning  spectrometer  was  built  whose  optic£tl  characteristics  match  the 
constraints  of  the  arc,  and  an  echelle  grating  was  used  to  achieve  high  spectral 
resolution  at  small  apertures.  The  perfoxnance  of  the  instrument  is  described 
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in  some  detail  since  high-order  spectroscopy  has  apparently  not  been  applied 
befoire  to  quantitative  measiarements  of  arc  plaismas. 

Provided  the  emission  and  absorption  coefficients  for  an  optically-thick 
line  are  invariant  over  the  spectral  bandpass  (because  of  high  spectral  reso¬ 
lution)  and  a^.ong  the  line  of  sight  of  the  spectrometer  (because  of  high  spatial 
resolution),  the  equation  of  radiative  tramsfer  yields  a  simple  analytical  solu¬ 
tion,  and  the  measured  emission  and  absorption  can  be  easily  interpreted  in 
terms  of  the  emission  and  absorption  coefficients  within  the  plasma.  Ihe  line 
profile  for  the  emission  coefficient  (i.e.  the  optically  thin  line)  can  then 
be  obtained  from  the  measured  line  profile  independent  of  any  a-prioid  assump¬ 
tion  regarding  the  broadening  mechnism,  and  one  can  also  determine  the  occupa¬ 
tion  of  the  upper  and  the  lower  energy  level  for  the  transition. 

As  shown  theoreticadly  in  Section  2A  and  demonstrated  experimentally  in 
Section  5^  this  information  can  often  be  obtained  by  high-resolution  emission 
measurements  alone  in  lieu  of  an  actual  absorption  measurement.  This  method 
is  rot  to  be  confused  with  the  check  for  self -absorption  customarily  used  in 
plasma  spectroscopy  and  based  on  calculating  the  Planck  function  at  the  electron 
temperature  in  order  to  estimate  the  absorptsince  from  the  peak  intensity  of  the 
line.  Our  results  for  the  Balmer  lines  demonstrate  how  misleading  such  a  pro¬ 
cedure  can  be  if  applied  to  non-equilibritim  plasmas. 

By  measuring  the  line  profile  one  can  determine  the  broadening  mechanisms, 
i.e.  obtain  information  on  the  environment  of  the  emitter.  Thus  we  derived  the 
kinetic  temperature  of  the  heayy^  particles  in  the  plasma  from  the  shape  of  the 
Doppler-broadened  rotational  lines.  The  shape  of  the  hydrogen  lines  is  deter¬ 
mined  predominantly  by  Stark  broadening  and  can  then  be  used  to  obtain  the  elec¬ 
tron  density  in  the  plasma.  In  this  work  the  shape  of  the  Stark-broadened  H 

Y 

and  lines  as  predicted  by  Griem’s  theory  (2-4),  has  been  directly  verified 
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Moreover,  the  calculated 


for  electron  densities  as  low  as  6  x  10  cm  ^ . 
shape  of  the  line  wings  is  probably  also  correct  at  electron  densities  which 
are  still  lower  by  an  order  of  magnitude* 

As  the  experimental  program  progressed,  it  became  obvious  that  our  initial 
optimism  regjaTding  the  attainment  of  equilibrium  in  the  arc,  was  ill-founded. 

Since  the  spectroscopic  observables  could  not  be  correlated  with  the  LTS  model, 
the  arc  plaisma  can  no  longer  be  characterized  by  a  few  parameters  -  temper¬ 
ature,  pressure,  and  chemical  composition  at  standard  conditions.  We  were 
faced  essentially  by  a  twofold  problem*  First  we  had  to  determine  experimen¬ 
tally  the  plasma  state,  i,e.  the  concentration  of  the  various  species  (molecules, 
atoms,  ions)  euid  their  distribution  over  energy  levels.  Secondly,  the  different 
interaction  processes  in  the  plasma  had  to  be  analyzed  in  order  to  approximate 
the  complex  plasma  state  by  a  simplified  state  model  and  to  elucidate  the 
mechainisms  brinering  about  the  observed  departuros  from  equilibrium. 

The  refined  diagnostic  information  made  available  by  high-resolution  measure¬ 
ments  in  emission  and  absorption  -  which  might  be  largely  redundaiit  for  LTB- 
plasmas  -  becomes  important  In  the  absence  of  equilibrium.  If  the  optical  transition 
probabilities  for  lines  belonging  to  a  certain  species  are  known,  the  spectros¬ 
copic  observables  cam  be  converted  into  the  population  distribution  over  energy 
levels  for  that  species,  without  involving  a-priori  assumptions  concerning  the 
plasma  state.  The  combined  emission-absorption  method  permits  one  to  extend  the 
range  of  these  measurements  to  the  lower  energy  levels.  In  emission  the  same 
data  can  either  not  be  obtained  at  all  (metastable  levels),  or  one  has  to  resort 
to  quantitative  measurements  in  the  vacuum  ultraviolet  which  are  notoriously 
difficult.  But  it  is  precisely  the  occui>ation  of  the  lower  energy  levels  which 
is  often  of  crucial  iwrortance  for  the  interpretation  of  the  plasma  state. 
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The  population  distribution  over  energy  levels  is  conveniently  expressed 
in  terms  of  the  occupation  of  energy  states ,  for  which  we  use  the  notation  n*. 
Thus,  if  n(i)  is  the  occupation  of  energy  level  i  of  statistical  weight  g(i)j 

n*(i)  =  n(i)/g(i)  (cm"^)  (l-l) 

Numerically  the  occupation  (also  called  population  density  or  number 
density)  is  expressed  in  cm"^.  The  representation  of  the  population  distri¬ 
bution  in  the  form  of  a  plot  of  log  n-^(i)  vs*  is  here  also  called  a 

"Boltzmaum  plot*' ,  but  this  does  not  imply  that  its  slope  is  necessarily  re¬ 
lated  to  a  temperature  in  the  true^  i.e.  thermodynamic  sense.  In  non- 
equilihrium  situations^  the  ''distribution  temperature"  introduced  by 
formally  equating  the  source  function  to  the  Planch  function  (Eq.  2-7);, 
does  not  refer  to  any  meaningful  temperature^  but  is  only  a  convenient 
parameter  to  characterize  the  population  distribution  between  two  energy  levels. 

As  shown  in  this  work,  the  velocity  distribution  of  electrons  and  of 
heavy  particles  in  arc  plasmas  is  often  still  Maxwellian  even  when  considerable 
departures  from  equilibrium  occur,  including  a  lack  of  equipartition  of  kinetic 
energy  between  free  electrons  and  heavy  particles.  One  can  then  properly  speak 

of  an  electron  temperature  and  a  gas  temperature  T  ,  However,  the  author 

e  g 

finds  the  widespread  use  of  such  expressions  as  excitation,  rotational,  or 

vibrational  "temperatures”  objectionable,  because  they  invoke  the  wrong  physical 

picture.  If  the  departures  from  LTE  are  relatively  mild  (as,  for  instance, 

in  the  low-current  arc  at  atmosi^eric  pressure),  the  population  distribution 

over  some  electronic  levels  will  be  approximately  determined  by  T  and  equi- 

e 

valently  over  some  rotational  levels  by  T^,  but  then  the  concept  of  "excitation 
temperature"  or  "rotational  temperature"  is  redundant,  at  best, 

*)  Our  results  are  mostly  expressed  in  terms  of  logarithms  to  base  10  ("log") 
natural  logarithms  are  designated  by  "In".  ’ 


If  the  departures  from  equilibrium  are  more  severe  (e,g,  in  the  arc  at 
reduced  pressure )>  the  "distribution  temperatures"  for  energy  levels  corres¬ 
ponding  to  a  particular  degree  of  freedom  are,  in  general,  different  for  the 
various  energy  levels  and  also  for  the  various  species*  In  some  cases  the 
"distribution  temperature"  is  practically  constant  over  a  range  of  energy 

levels,  leading  to  a  Boltzmann  plot  that  is  linear  by  sections*  However, 

a 

this  appears  to  be  only  a  spurious  effect  and  is  attributed  to/specific  re¬ 
lationship  between  the  pertinent  population  and  depopulation  rates  as  discussed 
in  Section  lOB  in  connection  with  the  observed  linearity  of  the  Boltzmann  plot 
for  predissociating  rotational  levels. 

In  order  to  illustrate  better  the  practical  application  of  our  method 
we  have  described  in  some  detail  how  the  Boltzmann  plots  for  hydrogen  and 
oxygen  were  obtained*  The  experimental  data  for  the  various  bands  of  OH  were 
represented  similarly  by  what  we  call  here  an  "apparent  Boltzmann  plot",  because 
the  reduced  intensity  shown  on  the  ordinate  is  not  simply  proportional  to  the 
occupation  of  the  rotational  energy  level  but  contains  also  the  -  unknown  - 
dependence  of  the  transition  moment  on  the  rotational  energy  term.  In 
Section  6D  we  discuss  the  functional  form  of  this  dependence  and  establish 
essentially  an  upper  bound  for  its  magnitude. 

In  the  second  part  of  this  report  the  observed  population  distributions 
are  interpreted  in  terms  of  a  theoretical  state  model  for  the  plasma.  The 
basic  concepts  used  in  the  theory  of  rate  processes  in  plasmas  are  described 
in  the  literature  (see  e.g.  Refs.  6  and  8).  Griem  has  used  theoretical 
estimates  to  derive  general  criteria  for  the  validity  of  LTE  in  plasmas 

Z) •  hydrogen  plasmas^  Bates  et  al.  have  obtained  the  general  steady- 

state  solution  of  the  rate  equations  corresponding  to  the  collisional- 
radiative  model  (^). 
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In  the  approach  followed  here  we  proceed  by  steps  using  the  detailed 
spectroscopic  information  avadlable  to  establish  a  correlation  between 
theoretical  deductions  and  experimentad  evidence.  It  is  readily  seen  that 
water  vapor  plasmas  are  eminently  suited  to  study  non-equilibrium  effects. 
Conclusions  of  general  validity  are  derived  using  as  representative  examples 
the  two  plasma  regimes  investigated:  (a),  a  water  vapor  plasma  in  an  argon 
matrix  generated  in  a  low-current  arc  at  atmospheric  pressure,  and  (b),  a 
plasma  of  pure  water  vapor  generated  in  an  high-current  arc  at  reduced  pressure. 

First,  we  attempt  to  simplify  the  complex  plasma  state  by  estimating  the 
relative  importance  of  the  various  elementary  processes  by  which  transitions 
between  energy  states  occur.  Hydrogen  is  chosen  for  these  model  calculations 
because  of  the  simplicity  of  its  energy  level  scheme  and  because  the  relevant 
rate  coefficients  are  at  least  approximately  known.  If  the  ground  state  is 
strongly  overpopulated  with  respect  to  its  occupation  in  equilibrium  with  free 
electrons  and  ions,  it  is  shown  that  the  population  and  depopulation  of  energy 
levels  is  determined  predominantly  by  step-wise  excitation  induced  by  electron 
collision.  The  occupation  of  energy  levels  falls  off  rapidly  with  increasing 
principal  quantum  number  which  leads  to  progressively  lower  "distribution  temp¬ 
eratures".  The  electron  temperature  becomes  then  a  rather  elusive  quantity. 

(in  a  sense,  this  situation  is  just  the  opposite  of  what  is  predicted  for  the 
collisional-radiative  model).  V/e  then  apply  the  step-wise  excitation  model 
to  the  calculation  of  the  occupation  of  the  hydrogen  energy  levels  in  the 
water  vapor  arc  at  reduced  pressure. 

The  observed  non-equilibrium  distributions  over  the  energy  levels  of  oxygen 
and  of  OH  can  also  be  explained,  at  least  qualitatively,  by  essentially  the 
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same  model,  although  encounters  with  heavy  particles  are  more  effective  than 
electron  collisions  in  inducing  transitions  involving  nuclear  motion.  The 
respective  Boltzmann  plots  are  similar  hecause  ionization  of  oxygen  and  disso¬ 
ciation  of  OH  can  both  occur  by  an  Atiger  process. 

As  shown  by  this  study  of  non-equilibrium  plasmas,  a  close  anailogy  exists 
between  ionization  by  electron  collision  and  dissociation  by  heavy  particle 
collision,  and  both  processes  are  doDiinated  by  step-wise  excitation.  The 
absence  of  equilibrium  between  bound  and  ionized  -  or  dissociated  —  states 
is  then  appropriately  discussed  in  terms  of  the  occupation  of  the  last  bound 
state.  This  concept  is  developed  in  Section  7  with  reference  to  the  ionization 
equilibrium.  Since  the  Saha  equation  is  only  a  specialized  form  of  the  general 
mass  action  eqiiation,  identical  arguments  apply  also  to  the  dissociation  equi¬ 
librium  . 

The  principle  of  detailed  balancing  or  microscopic  reversibility  (2)  applies 
rigorously  to  the  rate  coefficients  for  a  particular  process  and  its  inverse, 
but  not  to  the  process  itself.  Thus,  if  the  excitation  rate  for  atoms  or 
molecules  in  the  electronic  ground  state,  the  recombination  rate  for  an  electron- 
ion  pair,  or  the  dissociation  rate  for  a  molecule  such  as  OH,  are  smaller  than 
the  probability  for  a  competing  irreversible  process  -  diffusion  -  the  species 
concentration  in  an  isothermal  volume  element  of  the  plasma  can  deviate  con¬ 
siderably  from  the  equilibrium  concentration.  Based  on  the  particulate  natiire 
of  the  plasma,  Ke  use  gas-kinetic  considerations  to  estimate  the  effect  of 
diffusion  processes  and  the  lack  of  equipairtition  of  kinetic  energy  between 
electrons  and  heavy  particles* 
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2,  THEORY 


A,  Emission-Absorption  Measurements  of  Hcmiogeneous  Plasmas  at  High 
Spectrzd  Resolution 

If  scattering  can  be  disregarded,  the  change  of  intensity  of  a  besun  of 
monochromatic  light  of  wavelength  X.  incident  on  a  thin  layer  of  pleisma  of 
thickness  dx  is  given  by  the  difference  between  the  spontaneous  emission 
(emission  coefficient  €^)  amd  the  net  absorption  '  (net  absorption  coefficient  kj^) 
by  the  plasma » 

dlj^(x)  =  ^  ^ 

Introducing  the  source  function i 

5  ej^(x)/k'(x) 

Eq.  2-1,  the  equation  of  radiative  transfer,  can  be  written* 

dlx(x)  =  k^(x)  [S^(x)  -  l^(x)Jdx  (2-2) 

In  an  homogeneous  plasma  (e.g.  the  wall-stabilized  arc  in  end-on  measure¬ 
ments),  the  source  function  and  the  net  absorption  coefficient  are  Independent 

of  position:  Sj^(x)-S^{  k’(x)  =  k‘.  With  the  boundary  conditions  for  a  plasma 
of  length  Li 

lj^(x<0)  =  1^(0)  ,  Ij^(x>L)  .  I^(L) 

Bq,  2-2  can  easily  be  integrated,  yielding: 


I^{L)  .  S,[l  -  exp  (-kj^  l)J  +  I  (O)  exp  (-k*  L) 


(2-3) 


=  1®“  + 

“  X  "^X 


■*)  net  absorption  equals  absorption  minus  induced  emission,  which  should  be 
properly  considered  negative  absorption. 
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In  an  emission  measurement!  we  measure  •  In 

technique  used  for  absorption  measurements,  the  radiation  from  the  background 
source  is  selectively  modulated.  By  synchronous  rectification  of  the  ac  sig¬ 
nal  from  the  detector,  can  then  be  measured  without  interference  by  the 

dc  signal  due  to  plasma  emission.  the  intensity  of  the  light  from 

the  bctckground  source,  is  measured  independently,  and  we  obtain  the  absorp- 


tance 


[l  -  exp  (-kj^  L)] 


(2-4) 


Finally,  a  conciirrent  measurement  of  emission  and  absorption  yields  the 
source  functloni 

If  local  thermodynamic  equilibrium  obtains,  the  source  function  equals 
the  Planck  function  by  Kirchhoff's  law.  Ignoring  the  small  variation  of  X 
over  the  linewidth,  is  thus  constant  within  a  line.  This  is  still  true  in 
collision-dominated  non-equilibrium  plasmas  (i.e,  the  line  shapes  for  emission 
and  absorption  are  also  the  same),  but  S  is  then  given  by  the  principle  of 
detailed  balancing  as: 

S(line)  =  -  1  )  "^  (2-6) 

where  j  is  used  to  denote  the  upper,  and  i  the  lower  energy  level  for  the 
transition.  In  non-equilibrium,  a  distribution  temperature  can  be  defined 
by  formally  equating  the  source  function  to  the  Planck  function: 

exp  [  (E^  -  E^)/kTjj  j  3  n*(i)  /  n*(j)  (2-?) 
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Because  of  the  near-constancy  of  S  over  the  width  of  the  line,  low-resolution 
instruments  can  be  used  to  determine  the  distribution  temperature,  provided 
the  spectral  bandpass  remains  the  same  in  the  emission  and  absorption  measure¬ 
ments  (lO).  When  a  spectrometer  of  sufficiently  high  resolving  power  is  available 
such  that  a  record  of  the  true  line  shape  can  be  obtained,  the  spectral  radiance 
ovBT  the  line  (i,e,  the  profile  of  the  unabsorbed  line  corresponding  to  spontaneous 
emission)  can  be  constructed  point -by- point  from  the  spectral  scan  in  emission 

of  an  optically  thick  line,  since,  as  follows  from  Eqs.  2-3,  2-4,  and  2-5: 

jEm 

\  =  In  [(  1  -  la  [(i  -  (2-8) 

The  integral  of  N  over  the  line  is  related  to  the  occupation  of  the  upper 
energy  level  via  the  transition  probabilities  (Eq,  2-11).  Under  favorable 
circumstances,  an  emission-absorption  measurement  at  high  spectral  resolution  yields 
the  occupation  of  both,  the  upper  and  the  lower  energy  level  for  the  transition. 

The  precision  of  this  me-thod  depends  on  the  magnitude  of  the  absorptance.  At 
small  absorption  the  optically  thin  approximation  is  recovered,  since  N  wl®* 
for  0.  The  other  limiting  case  is  given  by  a  line  which  is  black  in  the 
center,  and  then  only  the  distribution  temperature  can  be  determined,  since 
w  for  a.^-^  1. 

A  spectral  scan  of  the  line  in  both  modes  is  not  required,  and  in  plasma 
measurements  the  source  function  can  be  obtained  by  determining  only  the 
absorptance  at  the  line  center.  (The  situation  is  reversed  for  most  flames 
where  considerations  of  noise  favor  absorption  measurements).  Often  the  source 
function  can  be  e^rtracted  from  high-resolution  measurements  in  emission  alone. 
Consider  the  case  of  two  adjacent  lines,  one  of  which  is  strong  and  the 
other  weak,  and  whose  relative  line  strengths  are  accurately  known.  Since 
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Russel**Saunder8  coupling  is  well  obeyed  by  the  li^t  atons  of  interest  in 

plasna  spectroscopy,  atcmic  aultiplets  can  easily  be  found  that  meet  these 

requirements.  The  case  of  the  ultraviolet  bands  of  OH  is  also  favorable 

because  of  the  presence  of  main  and  satellite  branches  in  the  rotational 

structure  which  form  several  convenient  line  pairs.  Their  relative  strength  is 

given  by  the  Hoenl-London  factor.  If  R  is  their  line  strength  ratio,  then 

N*  “  R  N",  The  strong  line  will  be  self-absorbed  in  the  center  if  we  find 
o  o 

I'  <  R  I”,  (single  primes  refer  to  the  strong,  double  primes  to  the  weak  line)*\ 
o  o 

Assuming  the  weak  line  to  be  optically  thin  (l^  w  H^),  we  can  write,  using 
Eq,  2-8i 

H  ~  K/K  ~  A  -  a;)  (2-9) 

It  should  be  noted  that  can  be  directly  derived  fro*  the  measured  ratio 
I^/l^  only  when  the  true  emission  profiles  have  been  resolved.  Once  has 
been  found  (e.g,  from  a  plot  of  the  function  on  the  right  side  of  Eq,  2-10), 
one  can  improve  the  accuracy  by  calculating  now  a”  and  iterating.  However, 
this  is  rarely  required  if  lines  are  chosen  for  which  R  >  10,  Since 
the  distribution  temperature  between  the  occupation  of  the  two  energy  levels 
follows  immediately. 

Once  the  source  function  has  been  determined,  the  experimental  profile  of 
can  be  converted  into  the  profile  of  by  using  Eq,  2-8;  this  is  not  neces¬ 
sary,  however,  if  only  the  half  width  -  or  any  other  fractional  width  -  is  re¬ 
quired,  since,  as  follows  from  the  same  equation i 

h*  /  “  [1  -  (1  -  h*  *  "o  /  ” 


*)  Here  euad  in  what  follows  the  subscript  for  the  line  center  has  been 
simplified,  and  the  emission  superscript  has  been  dropped. 
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The  half  width  hw  corresponding  to  N  ,  the  profile  of  the  optically-thin 

\ 

line,  can  thus  be  obtained  by  measuring  directly  the  width  of  the  recorded 
emission  profile  at  the  fractional  height  given  by  Eq.  2-10  with  n  =»  2, 


Optical  Transition  Probabilities  for  the  Rotational  Lines  of  OH  (A^S-X^H ) 


The  integrated  spectred  rauiiance  of  am  optically-thin  line,  corresponding 
to  a  transition  from  energy  level  j  to  energy  level  i,  emitted  by  an  homogeneous 
plasma  sample  of  length  L,  is  given  by: 


f  N,  dX 
line  ^ 


L  /  e.  dX 
line  ^ 


L  n(j)  A(j,i)  (2-11) 


in  energy  per  second,  per  square  centimeter,  and  per  steradian.  n(j)  is  the 
population  density  of  the  upper  energy  level  J,  and  A( j,i)  the  Einstein  spontaneous 
transition  probability. 

Consider  the  energy  levels  i  and  j  to  be  formed  by  a  group  of  closely 
spaced  energy  sublevels  i  and  j  ,  respectively.  These  sublevels  do  not 
necessarily  correspond  to  degenerate  states  of  equal  statistical  weight:  however, 
it  is  assumed  that  the  sublevels  j  are  so  close  together  that  their  relative 
occupations  are  given  simply  by  their  statistical  weights: 


n(.is) 


n(j)  g(.lj  /  E  g(.1  ) 


(2-12) 


Nj.(.1,i) 


2  Nji  ,1  ) 

L  s  r 

s,r 


--  L  n(j)  E  g(j^)  A(jg,i^)  /  2  g(j^)  (2-13) 

SfT  s 


By  equating-  Eqs.  2-11  and  2-13  we  obtain 


A(,i,l)  =  E  g(j  )  A(.i  i  )  /  2  g(.i  ) 

s,r  ®  ®  ""  s  ^ 


(2-14) 
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If  the  sublevels  j  correspond  to  truly  degenerate  states  of  the  energy 
s 

level  j,  then  gij)  =  1  for  all  s,  andi 
s 

A(j,i)  *  E  /  g(j)  (2-15) 

s,r 

For  electric  dipole  transitions  the  transition  probability  is  related  to  the 
square  of  the  matrix  element  of  the  dipole  moment  (tramsition  moment)  byt 

A(j,i)  “  2  |R(4,i-)l^  /  g(^)  (2-16) 

3  h  s.r  ®  ^ 

Because  of  the  more  complicated  structure  of  molecules,  errors  l>ave  often 

been  made  when  these  concepts  have  been  applied  to  molecular  spectroscopy*  The 

conceptual  framework  for  a  correct  quantitative  interpretation  of  intensities 

in  diatomic  molecular  spectra  has  recently  been  reviewed  in  detail  by  Tatum  (ll<)  * 

2  +  2 

The  ultraviolet  bands  of  OH  belong  to  a  E  -  transition*  The  selec¬ 

tion  rules  for  J  (A  J  =*  0,  +  1,  but  not  J»0<->J«0)  and  p  (only  states  of  different 
parity  combine)  are  strictly  obeyed.  For  Hund*s  case  (b)  we  also  have  the  selec¬ 
tion  rule  for  K  equivalent  to  the  selection  rule  for  J,  The  transitions  that 
satisfy  all  selection  rules  form  the  six  main  branches  of  each  vibrational  band* 
Since  the  X  II  state  is  intermediate  between  case  (a)  and  case  (b),  the  selection 
rules  for  K  are  not  strictly  valid  and  are  violated  by  the  six  wesdc  satellite 
branches*  The  energy  level  diagram  and  the  transitions  corresponding  to  the  12 
branches  are  shown  schematically  in  Fig*  1*  The  energy  levels  of  the  electronic 
ground  state  have  been  labeled  according  to  Hund*s  case  (b)*  Although  this  is 
a  H  state,  A-type  doubling  does  not  lead  to  doublet  lines  because  of  the  +  - 

selection  rule  (see  Fig.  l). 
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The  integrated  spectral  radiance  of  a  line  of  a  branch  is  then  given  according 
to  Eq,  2-11  byj 


NT(n*v'K'J',n"v''K"J") 


h  c  L 
4  IT  X 


n(n'v'K*J')  A(n*v'K’J’,n"v"K"J“)  (2-1?) 


where  n'  amd  n"  in  the  index  are  used  to  denote  the  upper  {hZ),  and  the  lower 
electronic  state  (X  H),  respectively.  Since  the  spin  splitting  is  small, 

Eq.  2-12  is  used  to  relate  the  occupation  of  each  sublevel  J’  to  the  occupation 


of  the  rotational  level  K' i 


n(n’v*K’J')  =  n(n'v*K') 


2J*  +  1 

(2S  +  1)(2K’  +  1) 


(2-18) 


The  spin  quantum  number  S  is  not  indexed  because  of  the  selection  rule* 

S  *  S  **  S.  The  statistical  weight  of  the  tJ*  sublevel  is  2J®  1  (No  yV""type 

doubling),  and  the  transition  probabilities  in  Eq.  2-1?  are  then  given  byj 


A-,  5  A(n'v'K'J",n"v"K"J")  =  — ^  E  |  R(n‘v*K' J‘m' ,n’VK"J'’m")p/(2J'+l)  (2-19) 

3nX-^  n*  ,ni'«  ' 

as  follows  from  Bq,  2-16,  m  is  the  magnetic  quEuitum  number. 

Based  on  the  Born-Oppenheimer  approximation,  the  matrix  element  can  be 
separated  into  a  radial  and®  angular  part.  The  integration  over  the  rotational 
coordinates  cam  be  carried  out,  amd  we  obtain i 


64  TT  S(K'J’,K"J’') 

- r  |R(n'v'J' ,n"v"J")p  - 

3  h  X-*  (2J'  +  l) 


(2-20) 


The  Hoenl -London  factors  or  rotational  line  strengths  S(K' J* ,K"J'' )  obey  the 


sua  rule I 

E  S(K'J',K”J") 
K".J" 


2J'  +  1  }  E  S(K'’J',K"J") 

K',J' 


2J"  +  1  (2-21) 
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The  Hoenl -London  factors  given  in  the  literature  have  to  be  used  uith 
caution  since  they  are  often  not  properly  nomaliaed,  thus  Dieke  and  Cross- 
Mhite’s  (12)  HL-factors  for  the  OH  bands  have  to  be  divided  by  8  in  order  to 
satisfy  Eqs.  2-21.  The  HL-factors  in  Table  4  of  Ref.  12  also  have  been  calcu¬ 
lated  with  the  sane  average  value  for  the  eouiding  constant  (a  ■  -7. 55)* 
dependence  of  a  on  rotational  quantum  number  for  the  v“0  level  of  the  electro- 
ni-j  ground  state  is  shoun  in  Fig,  2  for  K  <  20,  For  higher  K  values •  the 

interaction  between  K  and  S  has  to  be  taken  Into  accotxnt,  and  the  eqna- 
tioi.  used  for  this  computation  (Eq,  14  of  Ref,  12)  is  then  no  longer  valid.  The 
eff.Mct  on  the  HL-factors  of  the  increase  of  a  with  K  is  practically  negligible 
foi  the  strong  branches,  but  significant  for  the  satellite  branches . (see  Sec¬ 
tion  6g), 

The  validity  of  factoring  the  matrix  elements  in  Eq.  2-20* 

iR(n*v*J',n*'v"J'')|  =  J*  (v*  J* )  R^Cr)  ’F"(v"J’' )  dr  (2-22) 

in  order  to  assign  a  transition  probability  to  a  band  and  derive  Franck-Condon 

factors  for  the  different  bands  of  a  system,  has  been  critically  examined  by 

James  (l^),  and  by  Anketell  and  Learner  (14),  with  special  consideration  of  the 

OH  molecule.  Because  of  the  vau?iation  of  R  with  internucleao:  distance  (i.e,  its 

e 

dependence  also  on  v  and  J),  it  can,  in  general,  not  validly  be  replaced  by  an 
average  value  that  can  be  taken  out  of  the  integral.  Secondly,  as  denoted  in 
Eq,  2-22,  the  vibrational  wavef unctions  are  also  J-dependent  if  the  effect  of 
vibration-rotation  Interaction  is  taken  into  account. 

Although  rigorously  the  matrix  element  will  thus  be  different  for  each 
rotational  line,  its  dependence  on  J  is  a  second  order  effect,  and  it  appears 
plausible  that  the  matrix  element  can  still  be  represented  by  a  factor  that  is 
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constant  for  a  vibrational  band,  amd  a  correction  factor  that  depends  on  the 
rotational  tramsition  involved j 


iR(n'v'J',n"v"J'')l 


■(n'v'.n*’v")|^  T(J',J") 


(2-23) 


This  representation  of  the  matrix  element  is  useful  only  if  T(J*J")  cam 
be  approximated  by  a  simple  mathematical  expression  containing  the  rotational 
quantum  numbers.  This  is  indeed  the  case  for  the  rotational  lines  of  OH,  as 
shown  in  Section  6D  and  further  discussed  in  Section  12, 

Formally,  one  can  also  define  an  average  transition  probability  for  a 


vibrational  band* 


A(n*v' ,n”v") 


64  TT^  2 

- R(n'v',n'V)r 

3  h 


(2-24) 


This  concept  Is  not  useful,  however,  because  it  is  not  possible  to  assign 
unambiguously  an  averasre  wavelength  \  to  a  band.  The  same  reservation  applies 
to  the  equivalent  quantity:  "band  oscillator  strength". 

The  measured  integrated  radiance  of  a  rotational  line  is  expressed  here  in 
terms  of  the  state  occupation  of  the  upper  rotational  level  K',  the  "P-famtor" , 
and  the  (rotational)  "line  strength  factor".  As  follows  from  Eqs.  2-17,  2-18, 


and  2-20: 


hcL  , 

- -  n*(n'v'KV'  P^'r.'v'J*  ,n"v"J")  S(K' J*  ,K"J")A^ 

4Tr 


(2-25) 


where  the  state  occupation  'defined  by  Ec ,  l-l)  is  given  by: 


n*(n'v*K‘ ) 


and:  P(n'v' J' ,n"v"J" ) 


n(n' v’ K' ) 

(2S  +  l)(2K'  +  i) 


=  n(n'v'K'J'm* ) 


64  TT^  I  /  ,2 

-  |R(n'v‘J’ ,n"v"J")| 

3  h 


(2-26) 


(2-27) 
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The  "P-f actor"  defined  by  Eq.  2-2?  is  formally  equivalent  to  a  transition 
probability  times  the  cube  of  a  wavelength  (compare  Eqs.  2-24  and  2-2?),  and  is 
introduced  here  because  it  is  dimensionally  more  convenient  than  the  square  of 
the  matrix  element. 

Customarily,  the  classical  expression  involving  the  oscillator  strength  is 
used  to  relate  the  occupation  of  the  lower  level  for  a  transition  to  the  inte¬ 
grated  line  absorption  coefficient.  Here  we  use  the  source  function  to  establish 
the  relation  between  the  equivalent  quantities  obtained  by  an  absorption  and  by 
an  emission  measurement,  respectively! 

/  k*  dX  =  S'^  fe.dk  -  S'^  N^/L  (2-28) 

line  ^  line  ^ 

as  follows  from  the  definition  of  the  source  function  S  and  from  the  invariance 
of  S  over  the  line  (see  preceding  Section  2A).  By  substituting  the  expressions 
for  S  (Eq.  2-6)  auid  for  (Eq,  2-25),  we  obtain! 

/  K  dX  «  — -  n»(n"v"K")  P(n'v' J*  ,n"v"J")  S(K' J' ,K"J")  X  -  /  dX  (2-29) 

line  ^  3'^'=  line  ^ 


In  Eq.  2-29  we  have  written  n*(n"v"K")  insteaui  of  the  correct  factor 
[^n*(n"v"K")  -  n*(n'v*K')J  because  the  electronic  energy  term  for  the  transition 
OH  (A  S  -  X  H  )  Is  so  large  that  in  all  practical  situations  the  state  occupation 
of  the  lower  level  is  orders  of  magnitude  higher  than  the  state  occupation  of 
the  unper  level,  i.e.  stimulated  emission  is  insignificant  compared  to  absorption. 

If  the  integrated  absorption  coefficient  is  related  to  the  level  occupation 
(and  the  oscillator  strength),  the  ^-type  doubling  of  the  spin  sub-levels  J" 
has  to  be  taken  into  account! 


n*(n"v"K") 


n(n"v"K") 

2(2S  +  1)(2K"  +  1) 


n(n"v"rj"p"m")  (2-30) 
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Many  years  ago,  Condon  and  Shortley  (^)  have  already  stressed  the  symmetry 
achieved  by  representing  the  intensities  and  integrated  absorption  coefficients 
of  atomic  lines  in  terms  of  state  occupations  and  transition  moments  (called  "line 
strengths”  by  these  authors),  instead  of  in  terms  of  level  occupations  and  transi-’ 
tion  probabilities  (or  oscillator  strengths). 

It  should  be  noted  here  parenthetically  that  in  absence  of  LTE  the  level 
occupation  is  mostly  not  a  useful  quantity.  The  population  distribution  over 
excited  energy  levels  is  then  conveniently  discussed  in  terms  of  the  Boltzmann 
plot,  i*e,  the  representation  of  the  log  of  the  state  occupation  vs,  the  energy 
term.  The  slope  of  this  plot  is  determined  by  the  "distribution  temperature” 

(Eq.  2-6),  an  experimentally  measurable  pearamater.  The  pertinent  state  occu¬ 
pations  and  not  the  species  concentration  and  the  -undefined  -  partition  functions 
enter  the  mass  action  equations  for  non-equilibrium  plasmas  (see  Section  ?). 

Even  the  rate  coefficients  are  more  conveniently  expressed  in  terms  of  the 
state  occupations  (Section  8a). 

in  the  case  of  atomic  lines  the  two  representations  involve  essen¬ 
tially  only  a  matter  of  form,  Eqs.  2-25  ®nd  2-29  -  based  on  Condon  4  Shortley 's 
representation  -  appear  definitely  more  appropiate  in  the  case  of  molecular 
lines.  This  is  due  to  the  fact  that  -  in  addition  to  the  practical  advantage 
of  avoiding  the  complications  due  to  ,4 -type  doubling  -  the  radial  part  of  the 
transition  moment  can  often  be  conveniently  expressed  in  the  form  given  by 
2-23.  Attempts  to  introduce  a  wavelength-dependent  electronic  transition 
probability  or  oscillator  strength  (S^)  appear  rather  artificial. 

We  shall  often  use  the  quantity  formally  equivalent  to  the  radial  part  of 
the  transition  moment  of  the  vibration-  and  rotation-less  molecule i 

P°  -  P(A^L,0.0;x2n.0,l)  (2-31) 
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P®  is  essentially  the  "P-f actor"  corresponding  to  the  P-1  lines  of  the 
o 

(0,0)  band.  Bennet  and  Dalby  have  inf  erred  an  absolute  "band  oscillator  strength" 
for  the  (0,0)  band  -  f(0,0)«8x  10*^  -  from  lifetime  measurements  with  a 
wide  spectral  bandpatss -  (15)  •  Based  on  this  value ,  one  can  estimate  to  be 
of  the  order  of  3*3  x  10  s 


C.  Equilibrium  Composition  of  Water  Vapor  Plasmas 

In  the  temperature  range  of  Interest  (2000  -  6000K),  ionisation  is  negli¬ 
gible.  and  six  components  of  a  water  vapor  plasma  have  to  be  considered.  With 
their  partial  pressures  expressed  as« 

p(H-)  =  u  j  p(0,)  »  V  I  p(H-O)  -  w  t  p(0H)  -  X  I  p(0)  -  y  ;  p(H)  -  z  | 


we  can  write* 

„  1/2  „-l  „  „ 

^  1/2  „-l  _  „ 

X  U  W  IV^ 

z 

y  v*^^^  * 

u+v+w+x+y+z  •  P 

2u  +  z  -  4v  +  x  +  2y 

where  the  last  equation  follows  from  the  condition  that  the  total  concentra¬ 
tion  of  hydrogen  atoms  must  be  twice  the  total  concentration  of  oxygen  atoms. 

This  system  of  equations  was  solved  for  the  six  unknowns  on  an  electronic 
computer  for  various  pressures  and  temperatures,  using  two  sets  of  input  data, 
with  the  equilibrium  constants  taken  from*  (a)  the  composition  data 

for  a  water  vapor  plasma  at  1  atm  given  by  Burhorn  and  Wienecke  (l6)j  and 
(b).  the  data  compiled  in  the  JANAF  tables  (l7) «  Sizable  discrepancies  appeared 
in  the  output  data  that  were  largest  for  the  datum  of  major  Interest,  the  con¬ 
centration  of  OH.  amounting  to  as  much  as  60%  in  some  cases. 
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This  disagreement  was  traced  to  a  discrepancy  in  the  equilibrium  constant  for 
the  dissociation  of  OH,  Since  the  basic  input  data  are  not  indicated  in  Ref,  17, 
we  resorted  to  a  direct  calculation  of  this  equilibrium  constant. 

We  used  the  value  for  the  dissociation  energy  of  OH  accepted  at  present 
=»  35,501  cm  ^),  and  considered  the  effect  of  vibration-rotation  interaction 
in  the  calculation  of  Q(oh),  the  internal  partition  function  of  OH,  However, 
over  the  temperature  range  of  Interest  the  small  increase  with  temperature  of 
the  correction  factor  for  Q(0H)  due  to  vibration-rotation  interaction  (see  Eq, 
7-39  in  Ref,  18),  is  approximately  cancelled  by  the  increase  of  Q(o),  The 
contribution  of  the  excited  electronic  states  to  the  partition  function  is 
also  negligible,  and  K^COH)  can  then  be  represented  by  an  expression  in  closed 
form,  accurate  to  within  5%  for  T  <  6000  K  .  With  partial  pressures  in  atm 

io"22180/T  ^2-32) 

The  agreement  between  the  calculated  K^COH)  and  the  value  listed  in  Ref,  17 

was  within  1%  at  3OOO  K,  and  still  within  m  at  6000  K,  Therefore  we  conclude 
that  the  JANAP  tables  are  essentially  correct  whereas  the  calculation  of  Ref,  I6 
appear  to  be  in  error.  With  the  equilibrium  constants  taken  from  the  JANAP 
tables,  the  equilibrium  composition  of  water  vapor  plasmas  was  then  computed 
(in  steps  of  100  K)  at  the  following  pressures!  0,02>  0,03i  0,04|  0,05  1  atm. 

The  results  for  1  atm  and  0.02  atm  are  shown  grai*iically  in  Figs.  3  and  4. 
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D,  The  Fowler-Mllne  Method  In  Plasma  Spectroscopy 

In  equlllbriuB  plasmas  the  population  density  of  excited  states  relatlre 
to  the  species  concentration  Is  a  monotone-increasing  function  of  temperature. 
The  concentration  of  a  free  radical  such  as  OH  Increases  initially  with  temp¬ 
erature,  but  then  decreases  as  the  dissociation  process  into  atoms  begins  to 
predominate.  At  still  higher  temperatures  the  dissociation  loss  also  overrides 
the  increase  in  the  Boltsmsnn  factor  for  excited  states.  The  number  density  of 
emitters  reaches,  therefore,  a  maximum  value  n(nax)  at  a  specific  temperature  - 
the  "normal"  temperature  -  beyond  which  it  decreases  again.  For  LTE-plasmas 
of  given  chemical  composition  at  standard  conditions,  n(max)  and  the  normal 
temperature  can  both  be  calculated.  If  the  plasma  is  optically  thin,  n(max) 
for  a  certain  energy  level  corresponds  to  a  maximiua  in  line  Intensity  for  tran¬ 
sitions  originating  in  this  level.  This  is  the  basis  of  the  Fowler-Milne  method 
introduced  by  Larenz  in  the  dieignostics  of  arc  plasmas  (l2).  The  method  is 
appealing  because  no  prior  knowledge  of  the  transition  probability  is  required* 
in  fact,  it  should  be  possible  to  derive  both,  the  temperature  and  the  transi¬ 
tion  probability  from  such  a  measurement. 

The  Fowler-Milne  method  has  been  used  in  the  past  to  determine  the  temp¬ 
erature  profile  of  optically  thin  arc  plasmas,  if  the  Abel  inversion  of  the 
experimental  data  obtained  by  side-on  measurements  showed  that  the  radial  inten¬ 
sity  distribution  was  peaked  off  the  arc  center.  The  usefulness  of  this  pro¬ 
cedure  is  impaired  in  practice  because  small  experimental  scatter  in  the  meas¬ 
ured  intensities  or  deviations  from  radial  arc  syametxy  may  induce  a  spurious 
off-axis  peak  in  the  inverted  data  and  lead  to  temperatuxres  which  are  utterly 
wrong. 

The  noted  experimental  difficulties  can  be  avoided  in  end-on  measurements. 
Variation  of  the  arc  current  provides  a  simple  means  of  changing  the  core 
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teaperature  while  changes  in  the  intensity  of  a  spectral  line  are  directly 
nonitored  by  the  recorder  display.  Provided  the  arc  can  be  run  at  a  current 
large  enough  to  exceed  the  nomal  temperature,  it  is  then  possible  to  adjust 
the  arc  current  accurately  and  unambiguously  for  peak  intensity  of  the  spectral 
line.  As  opposed  to  the  customary  side^on  measurements,  the  temperature  regime 
of  interest  (normal  temperature  for  the  studied  spectral  feature)  is  shifted 
by  this  procedure  from  the  aurc  perijiiery,  where  measurements  are  less  accurate, 
to  the  arc  core. 

In  Fig,  5  w®  show  the  calculated  occupation  of  two  rotational  levels  of 
OH  in  a  water  vapor  plaisma.  It  cem  be  seen  that  over  a  wide  temperature  range 
the  number  density  of  emitters  (i.e.  the  line  intensity)  is  almost  constant. 
This  im^ies,  of  course,  that  the  Fowler-Milne  method  is  not  a  particularly 
good  (i,e,  sensitive)  method  to  determine  the  plasma  temperature. 

In  the  ease  of  Interest  this  argiment  has  to  be  inverted,  howeveri  if  the 
arc  current  is  adjusted  to  such  a  value  that  the  peak  temperature  in  the  arc 
core  is  close  to  the  normal  temperature  for  the  spectral  line  whose  transits.  & 
probability  we  want  to  measure,  the  line  intensity  will  be  constant  over  a  re¬ 
latively  wide  cross  section  of  the  arc  and,  moreover,  will  be  insensitive  to 
small  fluctuations  of  arc  temperature. 
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3.  APPARATUS  AND  EXPERIRENTAL  TECHNIQUES 
A,  High-Resolution  Scanning  Spectrometer 

The  requirements  for  a  spectrometer  matching  the  constraints  of  arc  sources 
used  in  plasma  spectroscopy  are  quite  different  from  the  traditional  design  cri¬ 
teria  for  commercial  spectrtmeters .  The  focal  length  of -these  instruments  is 
relatively  sraadl  as  compaured  to  the  size  of  the  grating?  but  in  end-on  measure¬ 
ments  of  arc  plasmas,  for  example,  f-numbers  of  the  order  of  f/lOO  axe  dictated 
by  the  requirement  to  isolate  a  narrom  beam  of  li^t  emerging  from  a  homogeneous 
plasma  region.  Thus  in  a  spectrometer  with  a  large  aperture  only  a  small  frac¬ 
tion  of  the  grating  surface  will  be  illuminated,  which  is  not  only  wasteful  but 
also  severely  limits  the  attainable  resolution*  Concurrently,  very  narrow  field 
stops  are  required,  and  the  entrance  slit  has  to  be  masked  down  to  a  fraction  of 
a  millimeter.  This  obviates  the  need  to  correct  for  image  curvature  in  a  Czemy- 
Tumer  instrument,  which  is  unavoidable  in  high  resolution  work  with  long  slits, 
Finally,  the  convenience  of  a  direct  wavelength  read-out  over  a  wide  spectral 
range,  although  important  for  a  general  purpose  spectrometer,  has  to  be  sacri¬ 
ficed  if  a  grating  is  to  be  used  in  several  higher  orders,  so  as  to  achieve  maxi¬ 
mum  resolution. 

He  have  built  a  Czemy-Turner  spectrometer,  specifically  designed  for  high- 

resolution  plausma  spectroscopy  and  whose  optical  parameters  -  f/40  al>erture  and 

5  meter  focal  length  -  match  the  characteristics  of  our  wall -stabilized  arc. 

For  work  in  the  near  UV  and  visible,  the  600  lines/mm  grating  blazed  at  54®  is 

used  in  higher  orders.  With  the  grating  fully  illuminated,  the  theoretical  re- 

o 

solving  power  e.g,  at  3G00  A  in  the  9th  grating  order  is  830,000.  The  actual 
performance  of  a  scanning  spectrometer  of  such  high  resolving  power  is  usually 
not  determined  by  the  optical  components,  however,  as  would  be  the  case  for  a 
spectrograph  using  photographic  detection,  but  rather  by  the  mechanical  quality 
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of  the  drive  Bechaolan,  'Rius  any  tmerenness  in  the  rotational  motion  of  the 

grating  ■*  doe  to  slip-stick  effects  of  the  bearing,  translational  errors  of  run 

of  the  pushrod  and  other  effects  -  has  to  be  small  on  an  angular  scale  equivalent 

0 

to  the  maximum  theoretical  resolution.  At  JOOO  A  (9th  order)  this  turns  out  to 
be  an  angle  as  small  as  0.3  seconds  of  arc. 

We  use  a  high-precision  grating  drive,  folloMing  a  design  by  MeCubbin  et 
al.  (go).  A  spectrum  is  scanned  by  rotating  the  grating  via  a  lever  am  moved 
tangentially  by  an  advancing  50  mm  micrometer  screw ,  The  lever  am  is  held  by  a 
rigidly  mounted  rolled  spring,  whose  pulling  force  is  independent  of  displace¬ 
ment,  against  a  pushrod  which  has  quartz  flats  attached  to  both  ends,  ‘Hie  quartz 
plates  press  agadnst  hardened  steel  balls  on  the  tip  of  the  micrometer  screw  and 
on  the  lever  am. 

The  angle  of  rotation  of  9°  provided  by  the  fine  grating  drive  covers  the 
blaze  of  the  grating  at  higher  orders  (i.e,  the  angular  range  of  maximum  grating 
efficiency).  Synchronous  electric  motors  and  a  gear  train  are  used  to  provide 
three  scanning  speeds  (ratio  IjIOiIOO),  the  slowest  speed  of  rotation  being 
13  seconds  of  arc  per  minute.  A  motor-driven  "coarse  drive",  and  an  awgMi ay 
scede  are  available  for  positioning  the  grating  at  €uiy  angle.  For  accurate  wave¬ 
length  settings  we  resort  to  the  mlcrcmieter  scale  of  the  fine  drive  and  to  refe¬ 
rence  lines  of  known  wavelength.  Using  the  technique  of  multl^e  overlaiiping 
orders,  the  few  strong  lines  in  the  Hg  I  speetrim  provide  enough  wavelength 
markers.  For  rigidity,  the  front  plate  containing  entrance  and  exit  slit,  the 
grating  drive  assembly,  and  the  two  Gzerny-Tumer  mirrors  are  mounted  on  a  heavy 
"ship-builder* s"  45-cn-wide  steel  channel.  The  simple  expedient  of  embedding  the 
steel  chemnel  in  sand  provides  isolation  against  mechanical  vibrations  transmitted 
throu^  the  floor.  The  instrument  cover  and  the  exposed  part  of  the  steel  channel 
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were  lined  first  with  urethane  foaa  for  themal  insulation,  and  then  with  black 
cotton  velvet  to  nininize  light  scattering, 

A  solenoid-actuated  fold-mirror  is  placed  behind  the  exit  slit  and  the  dis¬ 
persed  light  can  be  directed  to  either  of  two  photomultipliers  so  that  simulta¬ 
neous  measurements  can  be  made  over  a  wide  spectral  range  without  detector  inter 

chsmge,  an  important  feature  for  jdasna  spectroscopy.  A  dry-ice-cooled  RCA  7102 

o 

(S-1  response)  is  used  for  measurements  above  8000  A,  corresponding  to  low  gra¬ 
ting  orders  (2nd  and  3rd),  and  glass  filters  are  suitable  for  discriminating 
against  overlapping  orders.  For  the  near  UV  and  visible,  an  EMI  9558  QA  photo¬ 
multiplier  with  an  S-20  response,  especially  selected  for  low  dark  current  and 
high  cathode  sensitivity,  is  used  in  conjunction  with  a  simple  quartz  prism- 
monochromator  of  the  Bunsen-Kirchhoff  type,  as  order  sorter,  Ihe  collimator 
and  telescope  lens  are  also  mmde  of  quartz,  and  two  indexed  telescoping  mounts 
are  provided,  so  that  the  change  of  focal  length  with  wavelength  can  be  com¬ 
pensated.  Die  bandpass  defined  by  a  slit  of  width  s  in  the  focal  plane  must 
be  narrower  thw  the  free  spectral  range  in  the  grating  order  used.  For  a 
grating  used  in  higher  orders  near  its  blaze  *  blaze  wavelength  in  the 
first  order),  s  is  then  determined  byt 


s 


(3-1) 


where  D  is  the  linear  dispersion  of  the  priaa  monochrcwator.  In  order  to  mini¬ 
mize  stray  light  usually  a  much  narrower  slitwidth  than  given  by  Eq.  3-1  was 
used. 

Since  in  our  measurements  only  a  small  area  of  the  photomultiplier  cathode 
was  illuminated,  a  permanent  focusing  magnet  was  used  to  improve  the  performance 
of  the  detector,  eis  recommended  by  Knight  et  al.  (2i,).  Ibis  simple  expedient 


proved  to  t>e  effective:  the  gain  was  enhanced  hy  an  order  of  magnitude  and 
the  signal-to-noise  ratio  hy  a  factor  better  than  three.  The  signal  from  the 
photomultiplier,  modulated  by  a  1  KHz  light -chopper,  is  synchronously  recti¬ 
fied  by  a  lock-in  amplifier  and  displayed  by  a  strip  chart  recorder. 

Once  the  spectrometer  has  been  carefolly  aligned,  focusing  of  the  telescope 

Gzerny-Tumer  mirror  is  the  critical  instrument  adjustment  for  maximizing  the 

resolution.  This  was  accomplished  by  using  a  low-pressure  mercury  discharge 

as  the  source  and,  with  the  aid  of  a  microscope,  yisually  adjusting  for  optimum 

0 

sharpness  of  the  hyperfine  components  of  the  Hg  I  5wl  A  line  in  the  plane  of 
the  exit  slit. 

0 

The  hyperfine  structure  of  Hg  I  4358  A  is  more  convenient  for  resolution 
checks.  In  Fig,  6  a  typical  strip  chartreeord  of  this  line  is  shown,  obtained 
at  the  slowest  scanning  speed  with  a  slit  height  of  4  mm,  slit  width  of  0,030  mm, 
and  an  aperture  of  f/80.  The  theoretical  resolving  power  in  the  6th  order  for 

the  grating  set  at  an  angle  of  51®40’  and  illuminated  by  an  f/80  light  bean,  is 

.  0 
360,000  which  corresponds  to  a  resolution  of  12  mA,  Since  the  doublet,  whose 

o 

peak  separation  is  13  mA,  is  better  resolved  than  given  by  the  Rayleigh  criterion, 
the  theoretical  resolution  has  been  achieved.  In  Fig,  7  we  have  plotted  the 
measured  chart  positions  of  the  resolved  hyperfine  lines  vs,  the  wavelength 
increment,  and  the  straight  line  corresponding  to  the  calculated  chart  disper¬ 
sion  has  also  been  drawn.  The  insignificant  data  scatter  demonstrates 

the  smoothness  of  the  grating  drive.  Rather  unexpectedly,  almost  the  same 
precision  could  be  achieved  when  the  line  was  scanned  in  the  reverse  mode  (i.e, 
when  the  lever  arm  is  pulled  by  the  spring  instead  of  being  pushed). 
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After  several  noaths  of  intensive  nse  of  the  instroment,  a  sli^t  degra¬ 
dation  of  resolution  was  noticed,  and  could  he  traced  to  pitting  of  the  quarts 
flats  attached  to  the  pushrod.  These  can  be  quickly  replaced.  Because  of  the 
proximity  of  peissing  trains  (the  spectrometer  w«is  located  only  5P  meters  away 
from  an  active  railroad  track),  we  checked  for  a  possible  transmission  of  floor 
vibrations.  However,  while  sitting  on  the  peak  of  a  hyperfine  line  of  Hg  4358 
no  signal  variation  was  detected  when  a  high-speed  train  rushed  by.  The  slmide 
precautions  taken  in  order  to  isolate  the  spectrometer  from  environmental  effects 
appear,  therefore,  to  be  entirely  adequate. 

The  experimental  set-up  for  concurrent  emission  and  absorption  measurements 
is  shown  schematically  in  Fig.  8,  So  as  to  ensure  that  only  light  emitted  by  a 
homogeneous  plasma  sample  falls  on  the  detector,  we  had  to  limit  the  entrance 
slit  of  the  spectrometer  to  a  "pinhole"  of  0,1  x  0,25  ■*,  and  the  aperture  to 
less  than  f/lOO.  The  cathode  spot  of  the  xenon  aircs  used  as  background  sources 
for  absorption  measurements,  does  not  maintain  the  required  spatial  stability 
over  such  small  dimensions,  but  constantly  "wanders".  The  behavior  of  the 
available  popular  types  of  xenon  lamps  (Hanovia  491C39  and  959C98)  was  essen¬ 
tially  identical,  and  we  assume  that  this  is  an  intrinsic  chaxaeterlstic  of 
free-burning  arcs.  Moreover,  the  spectral  radiance  of  the  xenon  arcs  in  the 
UV  was  too  low  for  absorption  measurements  of  the  strong  OH  lines  at  high  spec¬ 
tral  resolution.  Because  of  its  stability  and  its  higher  intensity  in  the  OV, 
if  available,  another  wall -stabilized  arc  in  the  end-on  node,  run  at  high 
ciurrent  and  filled  with  argon  at,  say,  2  atm  pressure,  would  have  provided  am 
ideaO.  background  source  for  these  measurements. 

Different  masks  on  mirror  which  focuses  the  source  on  the  entrance 
slit  at  unit  magnification,  were  calibrated  amd  used  as  aperture  stops.  For  a 
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unlfora  source,  the  signal  mss  foujtd  to  he  exactly  proportional  to  the  lllonl- 
nated  area  of  Wider  apertures  were  required  for  intensity  calUnratiens  in 
the  UV. 

B,  Intensity  Galihration 

A  high-current  strip  lampwith  quarts  windou  (GB  30A/T24/13)  was  our  stand¬ 
ard  source  for  the  intensity  calihration.  The  teaperature  of  the  tungsten  strip 

was  determined  with  an  optical  pyrometer  recently  calibrated  by 

the  National  Btireau  of  Standards,  If  is  the  neastzred  tiNiperaiiire  in  K,  wo 

determine  the  true  temperature  from  the  following  relation  between  the  radiance 

o 

of  a  blackhody  and  of  the  striplanp  at  ■  653O  kt 

‘  ef  (3^2) 

The  effective  enissivity  is  given  by  the  enlssivity  of  tiuogsten  times  the 
transnittapce  of  the  quarts  window t 

-  0.92  0-3) 

Using  Wien's  law,  which  in  this  range  is  accurate  to  OJLSt,  we  have  with 
T»  s  10*5  j, 

/  9.569  (3.^) 

Eq.  3-2  and  were  solved  by  iteration  based  on  DeVos'  (^)  graphs  of 
S^(j(T,X),  which  had  to  be  extrapolated  to  3OOO  K,  These  grai^is  were  also  used 
to  calculate  the  radiance  of  the  strip  lamp  at  various  wavelengths. 

A  step-down  transformer  with  the  line  voltage  stabilised  by  a  Sola  trans¬ 
former  supplied  ac  power  to  the  striplanp.  This  simple  power  supply  was  ade¬ 
quate}  after  a  short  warm-up  period  the  excursions  in  temperature  fr«a  a  median 
"^true  *  ^  were  less  than  5  K,  This  is  only  slightly  higher  than  the  precision 

of  the  pjrrometer  reading.  The  striplamp  was  run  at  the  highest  permissible 
temperature  in  order  to  increase  the  radiance  in  the  DV,  An  available  carbon  arc 
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was  also  tested  Imt  finally  discaxded  becanse  its  neasured  radiance  In  the  UV  - 
the  only  spectral  region  where  its  nse  conld  be  adv2uitageous  -  was  significantly 
higher  than  calculated,  possibly  due  to  interference  by  CN  bands. 

Since  the  intensl'ty  of  the  strip  lamp  is  so  such  weaker  in  the  UV  than  at 
longer  wavelengths,  stray  light  constitutes  a  serious  problem  for  accurate 
calibrations  in  the  spectral  range  of  the  OH  bands.  In  our  spectrometer  spurious 
light  of  significant  intensity  can  cone  only  fxoa  overlapping  orders.  The  wave¬ 
length  of  the  stray  li^t  is  thus  known,  and  its  magnitude  can  be  determined 
with  suitable  filters.  TTie  stray  is  caused  by  random  scattering  hy  the  optical 
components  of  the  order  sorter,  and  our  tests  verified  this,  since  the  magnitude 
of  the  stray  was  proportional  to  the  width  of  the  exit  slit  of  the  order  sorter. 
When  working  in  the  UV,  spurious  light  was  minimised,  therefore,  by  using  slit- 

widths  much  narrower  than  determined  by  Eq.  3-I,  and  by  using  also  a  glass  filter 

0 

with  a  wide  bcuidpass  in  the  UV  and  a  cut-off  around  4100  A  (Coming  7-54) .  With 

0 

these  precautions  the  residual  stray  light  at  2810  A  (the  shortest  wavelength  at 
which  measurements  were  made  in  this  investigation)  was  determined  to  be  below  5^, 
an  acceptable  level . 

A  stray  light  problem  appeared  again  in  the  infrared  towards  the  limit  of 
the  useful  range  of  a  photomultiplier  with  an  S-20  response,  caused  in  this  ease 
by  the  decrease  of  quantum  efficiency  of  the  detector.  This  stray  could  be  easily 
eliminated  by  resorting  to  available  glass  filters  with  cut-on  in  the  visible. 
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C ,  Flasna  Source 


In  prior  related  work  we  developed  a  wall -stabilised  arc  which  includes 
severs^,  new  design  features  in  order  to  increase  the  versatilitj'  of  this  source 
zind  extend  its  useful  range  to  lower  pressures  and  also  to  the  study  of  gases 
which  before  have  been  difficult  to  maintain  in  a  (Quiescent  state,  e.g.  water 
vapor  plasmas.  Ihis  device  has  been  described  elsewhere  (23).  The  copper 
segments  forming  the  arc  channel  were  redesigned  in  order  to  reduce  the  channel 
diameter  from  6,3  to  4  mm.  This  was  necessary  in  order  to  stabilize  the  arc 
column  which  contracts  at  the  lower  arc  currents,  an  effect  especially  pco^ 
nounced  for  molecular  gases.  Leads  were  attached  to  the  end  plates  of  the 
arc-sections  facing  the  electrodes,  and  the  field  strength  in  the  sample  sec¬ 
tion  could  then  he  determined  by  measuring  the  voltage  drop  (see  Section  9C). 

When  using  a  power  supply  with  a  "drooping"  V-I  characteristic,  a  stable 
arc  can  be  maintained  only  if  the  slope  of  the  V-I  characteristic  of  the  load 
(arc  +  series  resistor  of  resistance  fi)  is  positivei  the  constraint  being,  that 
the  total  voltage  drop  cannot  exceed  the  open  circuit  voltage,  i.e. 
dVf  dV  ^ 

-  -  -  +  B  +  I_>0  I  O.C.V.  >  V  +  I  R  (3-5) 

dl  dl  dl  ^ 

At  high  currents,  the  V-I  characteristic  of  the  arc  is  positive,  and  a 
very  stable  arc  can  be  mintained  with  a  battery  of  Miller  rectifiers  (23). 

At  low  currents  the  slope  becomes  negative,  however,  and  a  series  resistor  has 
to  be  included  in  the  circuit;  but  stabilization  at  4  Amps,  became  critical 
when  a  water-cooled  resistor  (dR/dl  **  O)  was  used*  A  substantial  improvement 
of  arc  stability  was  achieved  by  using  instead  a  battery  of  incandescent  light 
bulbs  as  a” dynamical”  series  resistor  (dR/dl  >  O),  Thus,  for  example,  the  in- 
tensity  of  Ar  I  4300  A  remained  constant  to  within  over  a  run,  and  the 
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run-to-run  repeatability  at  constaat  corrent  waa  of  the  order  of  By 
taking  recoorse  to  a  very  simple  strategem,  the  stability  of  our  are  matches 
the  best  performance  of  arcs  reported  in  the  literature  and  achieved  by 
using  elaborate  feed-back  controls. 

D.  Accurate  Arc  Alignment  for  End-On  Measurements 

In  end-on  measurements  the  arc  axis  must  be  aligned  parallel  to  the 

optical  axis  of  the  system,  so  that  the  spectrometer  views  a  homogenous  plasma 

sample  when  appropiate  field  and  aperture  stops  are  used.  We  achieve  this  with 

the  aid  of  a  small  laser  positioned  in  front  of  the  exit  slit.  Proof  of  correct 

eurc  alignment  can  be  obtained  hy  comparing  the  results  of  end-on  measurements 

with  the  equivalent  inverted  data  derived  frffla  a  spatial  scan  in  the  side-on 

mode  (23).  this  test  is  quite  laborious,  however,  and  a  different  method  was 

used  that  requiires  only  measurements  in  the  end-on  mode. 

The  arc  is  filled  with  argon  at  atraos|iierie  pressure,  the  spectrometer  is 

o 

focused  on  the  arc  center,  and  an  argon  ion  line  -  e,g.  Ar  II  4348  A  -  and  the 

adjoining  continuum  are  measured  at,  say,  three  different  settings  of  the  arc 

current.  Both  spectral  features  are  optically  thin.  For  the  data  reduction 

we  use  three  graphs  derived  from  the  calculated  composition  of  an  argon  plasma 

at  1  atm,  and  representing  as  a  function  of  temperature  (a)  the  integrated 

o 

spectral  radiance  of  Ar  II  4348  A,  (b)  the  spectral  radiance  of  the  continuum 
0 

at  4350  A  (both  per  unit  optical  path  length),  and  (c)  the  ratio  of  these  two 
queuitities, 

Graih  (c)  is  used  to  derive  the  "core  temperature",  i.e,  the  temperature 
of  the  plasma  sample  corresponding  to  the  arc  axis,  from  the  measured  signal 
ratio,  (ihe  instrument  response  remains  invariant  over  the  small  wavelength 
range  covered,  and  the  intensity  ratio  is  independent  of  optical  pathlength). 
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Once  the  core  temperature  has  been  established,  the  effective  optical  path 
length  L*  is  obtained  from  the  other  t¥o  plots,  Ijrpical  results  of  these 
measurements  are  shomn  in  Table  I. 


TABLE  I 

CORE  TEMPERATURE  AMD  EFFECTIVE  OPTICAL  PATH  LENGTH  FOR  AN  ARGON  ARC  AT  1  ATM 


I 

(Amp) 

T 

(K) 

L»(Ar  11) 
(cm) 

L*(cont. ) 
(cm) 

T(24) 

(i^ 

35 

11,950 

10.9 

10.7 

11,640 

50 

12,750 

11.2 

10,7 

12,320 

75 

13,500 

13.5 

13.8 

13,230 

At  given  pressure  and  chemical  composition  of  the  plasma,  the  core  temp¬ 
erature  of  a  wall-stabilized  arc  is  uniquely  determined  by  its  diameter.  Our 
measurements  can  thus  be  compaired  with  the  equivalent  data  reported  by  Bues 
et  al.  (^)  for  a  mall -stabilized  atrc  of  almost  the  same  relevant  dimension 
(4.10  mm  arc  diameter  in  Ref,  24  and  3.97  nn  in  our  arc).  The  good  agreement 
between  both  sets  of  temperatures  is  yet  improved  when  corrections  are  made 
for  the  slight  difference  in  arc  diameter. 

More  significant  is  the  agreement  between  the  derived  effective  optical 
path  lengths  and  the  actual  arc  dimensions^  which  were  10. 7  cm  for  the  length 
Of  the  are  channel  (test  section),  and  14.8  cm  for  the  total  arc  length  (dis¬ 
tance  between  electrodes),  Otir  measurements  are  precise  enough  to  bring  out 
the  fact  that  at  the  higher  currents  the  hot  are  core  extends  beyond  the  con¬ 
fining  channel  into  the  electrode  chambers. 
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Based  on  these  results,  it  can  be  conelnded  thatt 

(a)  The  axe  plasma  is  in  local  thermodynamic  equilibrium. 

(b)  The  correct  quantities  have  been  used  for  the  transition  probabili'^ 

8  —1 

of  Ar  II  4348  (a^  m  1.37  x  10  s  ),  and  for  the  representation  of 
the  continuum  intensity  (Eq.  5'"5). 

(c)  Good  precision  has  been  achieved  in  absolute  intensity  calibration. 

(d)  The  arc  axis  is  coincident  with  the  optical  axis, 

(e)  The  field  and  aperture  stops  used  are  appropiate  for  limiting  the 
li^t  flux  falling  on  the  detector  to  the  lif^t  emitted  by  the 
homogeneous  axe  core. 

Portuitious  agreement  appears  unlikely  in  view  of  the  relatively  wide 
temperature  range  covered  by  the  three  concurrent  measurements  and  the 
different  temperature  dependence  of  the  intensities!  The  intensity  of  Ar  II  4348 
increases  by  a  factor  of  about  20  when  the  current  is  raised  froa  35  amperes  to 
75  amperes,  the  intensity  of  the  continuum  only  by  a  factor  of  4, 

As  a  final  check,  we  made  use  of  the  high'-resolution  capabilities  of  our 

o 

instrument  to  derive  the  Stark  broadening  paurameter  for  the  Ar  II  4348  A  line 
from  the  measured  line  shapes,  as  shown  in  Table  II 

TABIiE  II 

MSASUHSD  STARK  BROADENING  PARAMETER  FOR  AR  II  4348  A 


I 

T 

n 

e 

P 

0 

1 — 1 

XI 

(U 

CO 

<< 

<1 

(Amp) 

(K) 

1 

g 

(-A) 

(.A) 

35 

11,950 

6.33 

163 

0.87 

142 

22.4 

50 

12,750 

9.54 

231 

0.93 

215 

22.5 

75 

13.500 

12.8 

302 

0.96 

290 

22.6 
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The  neasored  line  shape  is  a  eoDTOlation  of  the  Stark  profile,  which  is 
approxiaately  Lorentsian,  and  of  two  Caossians  (Doppler  and  instruaent  broad¬ 
ening).  Ihe  Doppler  half  width  followed  from  the  known  teaperatnre,  amd  the 
instmaent  broadening  was  deterained  experiaentally  hy  a  concurrent  scan  at 

the  saae  instmaent  setting  (i.e,  saae  grating  order,  slit  widths,  and  aperture) 

e 

of  the  hyperfine  stmctore  of  Hg  I  4358  A,  generated  by  a  low-pressure  aercury 
arc.  Based  on  the  Voigt  profile,  a  correction  factor  P  was  calculated  in  order 
to  derive  AXg  froa  the  measured  (full)  halfwidth.  Since  P  is  close  to 

unity,  the  precision  of  aeasuring  Alg  is  not  sensibly  affected  by  the  correction 
for  other  broadening  aechtmisas. 

The  Stark  broadening  parameter  is  obtained  by  dividing  "the  equili- 

briua  electron  density  calculated  for  an  argon  plasma  at  1  ata  pressure  at  the 

measured  teaperature.  The  results  are  shown  in  the  last  column  of  Table  II,  eon- 

venlently  nonudized  for  n  -  10  cm  It  follows  from,  theory  that  the  Stark 

broadening  parameter  shoidd  be  essentially  constant  (i.e,  a  function  of  n  alone), 

o 

and  increase  only  sli^tly  with  temperature,  Ihis  is  borne  out  by  the  results 
of  our  aeasureaents. 

Bo  precise  experiaental  values  for  the  Stark  broadening  parameter  of  Ar  II 
o 

4348  A  are  presently  available  with  which  to  eoapare  our  results.  It  is  only 
hnown  from  aeasureaents  by  other  investigators  that  Griea's  calculated  values 
for  the  argon  ion  lines  are  probably  too  low  by  a  factor  of  the  order  of  three. 

Our  Stark  broadening  parameter  la  3.8  tines  larger  than  calculated  by  Grien.(4). 

It  appears  that  these  results,  beyond  validating  our  approach,  clearly  de¬ 
monstrate  that  by  high-resolution  measurements  of  a  hcanogeneous  juad  tlne- 
Indejiendent  (steady-state)  plasma  sample,  impoartant  information  can  be  obtained 
with  higher  precision  and  considerably  less  labor  than  has  heretofore  been 


possible. 


k,  SIDB-OH  HBASURE»ENTS  AND  ABEL  INVERSION 


Die  are  «as  set  up  with  its  axis  parallel  to  the  entraace  slit,  and  spatially 
scanned  at  the  location  of  the  alddle  slndoN.  The  lateral  intensity  profiles  of 
the  Ej^-2  (0,0)  and  Rj^-20  (O,©)  lines  were  recorded  while  viewing  a  water  vapor 
are  at  redneed  pressure.  The  use  of  these  adjacent  lines  Is  advantageous  for 
detetnlnlng  the  tenperature  fron  their  intensity  ratio,  because  of  the  relatively 
wide  span  in  upper  energy  levels,  and  because  the  two  lines  are  of  coiparable  in¬ 
tensity  in  the  tenperature  range  of  our  neasurenents.  Possible  errors  due  to 
self -absorption  are  thus  nlnialsed.  Their  relative  transition  probabilities  were 
set  equal  to  their  Hoenl-London  factors,  which  is  net  quite  correct  (see  Section 
6d).  The  snail  absolute  error  in  the  derived  tenperatures  is  inconsequential  in 
this  context. 

In  order  to  convert  the  lateral  intensity  profile  into  the  radial  intensity 
distribution,  we  used  two  nunerical  nethods  for  the  Abel  inversion,  one  due  to 
Bockasten  (25)^  and  the  other  to  Barr  (26).  In  all  cases,  each  half  of 
the  Intensity  profile  waua  divided  into  20  sones,  corresponding  to  Increnents  of 
0,2  nn.  The  spatial  resolution  of  each  neasurenent  was  0.1  an  as  detemined  by 
the  width  of  the  entrance  slit.  The  natrix  nultiplicstion  was  performed  on  a 
high-speed  digital  computer. 

The  center  of  the  distribution  was  detexwined  as  follows t  Before  the  aetual 

o 

run,  a  lateral  intensity  profile  for  the  continuun  at  3^83  A  enitted  by  sn  argon 
plasna  at  1  atm  pressure,  was  obtzdned.  By  "folding"  this  profile  -  which  is 
narrow  and  quits  syanetrie  -  the  center  can  be  picked  out.  The  intensity  profiles 
for  the  rotational  lines,  obtained  next  in  a  run  with  water  vapor  at  redneed 
pressure,  is  wuch  broader  because  of  are  bulging,  and  shows  slight  assywetrles. 

It  is  thus  considerably  more  difficult  to  detemine  the  center.  Agreement  with 
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the  center  found  in  the  argon  run  is  obtained  if  the  profiles  are  folded  to 
natch  the  wings.  Since  this  is  also  the  procedure  used  by  other  investigators, 
this  was  the  criterion  adopted  to  fix  the  distribution  center. 

The  results  of  the  aeasurenents  of  a  water  vapor  arc  at  50  Aaps.  and  22  Torr 
are  shown  in  Figs.  9  and  10,  In  Pig,  9  the  radial  distribution  of  the  enission 
coefficient  for  the  Rj^«20  line  is  shown  as  derived  from  the  sane  input  data  by 
using  Barr's,  and  by  using  Bockasten's  nethod.  These  data  have  been  coabined 
with  the  equivalent  distribution  of  the  P^-2  line  and  converted  into  the  radial 
temperature  profile  (Fig,  10),  A  smooth  distribution  is  obtained  by  using  Barr's 
method,  but  the  other  method  leads  to  considerable  data  scatter.  This  is  attrib¬ 
uted  not  only  to  the  absence  of  a  built-in  data  smoothing  scheme  (Barr  uses  a 
least-squares  fit),  but  also  to  the  fact  that  Bockasten  fits  the  experimental 
data  to  an  analytical  expression  which  does  not  reproduce  well  the  points  of 
zero  slope.  In  effect,  the  scatter  is  largest  at  the  center  of  the  distribution 
and  at  the  off-axis  peak,  which  are  precisely  the  regions  of  greatest  interest. 

Because  of  the  presence  of  a  pronounced  off-axis  peak  -  a  situation  rarely 
encountered  for  the  atomic  lines  used  in  plasma  spectroscopy  -  our 
measurements  provide  an  excellent  case  for  testing  different  methods 
for  the  Ahel  inversion.  Direct  comparison  with  end-on  measurements  is 
also  feasible. 

The  arc  was  set  up  again  in  the  end-on  mode  and  properly  aligned  with  the 
optical  axis  of  the  system.  The  arc  was  then  spatially  scanned  along  its  hori¬ 
zontal  diameter  and  the  line  intensities  were  recorded  in  intervals  of  0,25  ®m, 

Ihe  dimensions  of  the  entrance  slit  were  0.1  mm  wide  by  0.25  ran  high,  and  the 
beam  aperture  was  stopped  down  to  f/300.  Thus  in  each  measurement  an  area  ele- 
merit  of  the  order  of  only  0,1  aua  was  "seen"  by  the  spectrometer.  Since  this 
amounts  to  less  than  of  the  total  cross  section  of  the  arc^  a  sufficiently 
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high  spatial  resolution  has  heen  achieved. 

Two  ty^eal  Intensity  distributions  are  shown  in  Fig,  11,  These  neasnre- 
aents  not  only  substantiate  the  assumption  of  cylindrical  are  syametry,  but 
also  basically  agree  with  the  radial  Intensity  distribution  derived  by  side-><» 
neasurements.  Side-on  aeasureaents  are  affected,  however,  by  "arc  bulging" 
at  the  location  of  the  small  gap  between,  arc  sections  (23),  This  effect  bo- 
ecaies  very  pronounced  at  :redueed  pressure,  and  the  arc  expands  to  about  twice 
the  channel  diameter.  Thus  in  side-on  measurements  the  intensity  peaks  appear 
slightly  shifted  towards  the  arc  perijhery. 

The  agreement  between  the  radial  temperature  profiles  obtained  by  end-on 
measurasents  (Fig.  12),  and  derived  from  side-on  measurements  (Fig,  10 )  can 
only  be  called  fair,  however.  More  credence  is  given  to  the  results  of  the 
end  -on  mode  because  the  two  lines  were  measured  concurrently  for  each  are  po¬ 
sition,  whereas  in  the  side-on  mode  the  lateral  intensity  distribution  for 
each  line  was  measured  separately,  for  reasons  of  expediency,  A  slight  increase 
of  are  bulging  while  the  11^-20  line  was  being  measured,  could  have  caused  a 
spurious  inereause  of  the  measured  intensities,  leading  to  higher  values  for  the 
derived  temperature  profile.  The  amplification  of  errors  produced  hy  the  Abel 
inversion  may  also  contribute  to  the  distortion  of  the  temperature  profile  in 
the  side -on  mode. 
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5.  DIAGNOSTIC  £ND-ON  MEASUREMENTS 


A,  Choice  of  Arc  Regiae 

In  the  measurements  at  atmospheiric  pressure,  the  sample  section  of  the  euc 
WE«  filled  with  argon  saturated  with  water  vapor  at  ro«i  temperature  by  bubbling 
argon  through  an  150  cm  high  column  filled  with  water  and  containing  a  fine  wire 
maze.  The  intensity  of  a  strong  rotational  line  waa  monitored,  while  the  arc 
current  was  decreased  by  steps.  The  Q^-6  (0,0)  line  became  measurable  only  at 
10  Amps.,  and  its  intensity  increased  rapidly  with  decreasing  arc  current,  whereas 
the  underlying  continuum  -  essentially  proportional  to  the  square  of  the  electron 
density  -  showed  opposite  behavior.  Thus  at  4,1  Amps,  the  line  was  8  tines 
stronger  than  at  8  Amps.,  but  the  continuum  was  only  I/3  as  strong.  These  data 
correspond  to  the  arc  core.  When  the  arc  at  4.1  Amps,  was  laterally  displaced, 
so  that  the  spectrometer  viewed  a  region  close  to  the  arc  perii*ery,  the  line 
intensity  increased  further  by  a  factor  of  the  order  of  3,  whereas  the  intensity 
of  the  continuum  decreased  to  I/6,  When  in  what  follows  Ke  refer  for  short  to 
the  are  at  atmosirfieric  pressure",  we  mean  the  4.1  Amp.  arc  burning  in  argon 
saturated  with  water  vapor. 

Our  attempts  to  run  a  low-current  are  with  pure  water  vapor  were  unsuccess¬ 
ful.  Although  the  condensation  problem  was  solved  by  running  the  are  "hot" 

(i.e.  the  inlet  temperature  of  the  pressurized  cooling  water  was  close  to  100°C), 
the  voltage  required  to  maintain  a  quiescent  arc  exceeded  oar  capability  of 
900  V.O.C. 

By  reducing  the  pressure  to  about  20  Torr,very  stable  are  operation  was 
achieved  with  pure  water  vapor  as  the  sample  gas.  The  plasma  expands,  filling 
the  arc  channel,  and  the  OH  bands  are  strongly  excited  at  relatively  high  currents 
corresponding  to  a  positive  V-I  arc  characteristic. 
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In  order  to  optlnize  the  arc  reglae  for  the  planned  aensojrenents,  the  intensity 
of  the  0^-6  (0,0)  line  was  a^n  nonitored,  while  the  are  current  was  varied  at 
constant  pressure.  The  plot  of  line  intensity  vs,  enzTent  corresponding  to  a 
pressure  of  18  Torr,  shows  a  broad  naxlnum  centered  at  about  37  Anps,  dropping 
to  $0^  of  peak  value  at  10  Aaps.  and  60  Anps.,  respectively.  It  was  also  noted 
that  the  line  intensity  depended  but  weakly  on  the  are  pressure,  whereas  in  LTE 
the  number  density  of  eaitters  should  increase  as  the  square  of  the  pressure. 

After  this  preliminary  work  had  circumscribed  the  arc  regimes  to  be  used, 
two  types  of  spectroscopic  measurements  were  made.  In  what  we  call  "Diagnostic 
Measurements"  our  aim  was  to  characterise  the  state  of  the  plasma  as  complete¬ 
ly  as  possible  by  resorting  to  emission  and  absorption  measurements  and  to  line 
shape  studies  of  suitable  atomic  lines  and  of  sene  rotational  lines.  In  another 
series  of  measurements  we  recorded  the  intensities  of  all  useful  -  i.e,  not 
overlapped  -  rotational  lines  of  the  OH  bands.  Since  it  took  about  six  hours 
to  complete  a  scan  of  the  entire  OH  spectrum  in  the  UV,  these  measurements  nay 
appropiately  be  labelled  "Long  Runs".  Adequate  arc  stability  could  be  maintained 
for  such  extended  time  periods,  bat  the  endurance  of  the  experimenter  was  severe¬ 
ly  tested.  After  two  complete  scans  had  been  recorded,  the  data  acquisition  tine 
was  shortened  in  subsequent  "Long  Runs"  by  starting  the  scan  with  the  (0,0)  band. 
The  short  wavelength  range  -  i.e,  the  region  of  the  (l,0)  and  (2,1)  bands  -  was 
skipped  because  it  required  repeated  adjustments  of  the  order  sorter  emd  also  a 
painstakingly  careful  intensity  calibration. 

In  a  third  series  of  measurements,  12  selected  line  groups  were  measured 
which  comprise  54  rotational  lines  belonging  to  the  five  OH  bands  studied.  Lines 
of  different  branches  from  low  to  high  K-numbers  were  included.  It  took  a  rela¬ 
tively  short  time  to  scan  these  line  groups,  and  measurements  were  made  at  various 
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arc  currents  in  an  attempt  to  apply  the  Fowler-Milne  method  not  only  to  one  line, 
but  to  several  lines  belonging  to  different  bands.  The  later  data  analysis  re¬ 
vealed,  however,  that  this  method  is  not  applicable  because  of  severe  departures 
from  LTE.  Thus  no  new  information  was  obtained  by  these  measurements,  except 
that  of  corroborating  the  results  of  the  "Long  Runs”, 

B.  Measurement  of  Electron  Densities 

The  theory  of  Stark  broadening  of  the  hydrogen  Balmer  lines  provides  the 

theoretical  basis  for  a  convenient  method  to  determine  the  electron  densities 

in  ilasmas  (2-4),  Excellent  agreement  between  measured  and  predicted  line 

shapes  has  been  reported  for  experiments  in  which  n  could  be  independently  de- 

e 

termined  (22,28).  Ihe  theoretical  predictions  have  not  yet  been  reliably  con¬ 
firmed,  hOKever,  for  <  5  x  10^^  cm"^.  It  is  precisely  in  the  lower  range, 
corresponding  to  departures  from  LTE  in  plasmas,  where  this  method  would  be  most 
useful  because  the  electron  density  can  then  no  longer  be  inferred  unambiguously 
from  other  measurements. 

Since  Hp  -  which  is  generally  preferred  -  showed  substantial  self -absorption 
in  most  of  our  measurements,  we  used  the  and  lines  of  the  Balmer  series  as 
diagnostic  lines.  Because  of  the  peculiar  shape  of  whose  center  has  then 
the  appearance  of  a  narrow  cusp,  its  peak  intensity  cannot  be  determined  with 
precision  at  low  electron  densities,  Moderat®  self -absorption  which  has  an  al¬ 
most  insignificant  effect  on  the  integrated  line  intensity,  will  "clip"  the  line 
peak.  It  will  be  further  distorted  by  an  experimental  slit  function^  small  com¬ 
pared  to  the  half  width  of  the  line  but  wide  compared  to  the  narrow  peak.  The 
theoretical  uncertainties  are  also  largest  close  to  the  center,  where  the  line 
is  broadened  by  electron  impact.  The  theory  of  quasi-statie  broadening  by  ions. 


which  determines  the  shape  of  the  line  wings,  is  more  firmly  established.  It 
is  therefore  not  surprsing  that  techniques  essentially  based  on  a  normalization 
at  the  line  peak,  such  as  inferring  the  electron  density  from  the  half  width  or 
from  other  calculated  fractional  intensities  (22)f  yielded  values  that  obviously 
were  too  high* 

In  order  to  compare  the  recorded  line  shape  with  the  theoretical  Stark  pro¬ 
file,  we  normalize  the  integrated  line  intensity* 

1.  S  v^Ak)  d^Ax.)  -  1  (5-1) 

^  line 

o  o 

where  AX  is  the  distance  from  the  line  center  in  A,  and  A  «  u  x  A,  is  the  meas¬ 
ured  area  under  the  line.  The  signal  (relative  line  intensity)  V  is  measured 
in  arbitrary  chart  units  u.  If  other  broadening  mechanisms  can  be  disregarded, 
a  one-to-one  correspondence  can  be  established  between  the  experimental  line  pro¬ 
file  -  the  integrand  in  Eq,  5-I  -  and  the  normalized  Stark  profile  for  this  line* 
The  latter  is  conveniently  expressed  in  terms  of  the  a-paraneters* 

00 

r  3(a)  do  =  1  {  a  =AX/F^  (5-2,3) 

With  wavelength  in  A  units,  the  normalized  field-strength i  F  **  1.253  lO"*^  n 

o  e 

In  Griem's  book  (4)  ealculated  values  for  the  reduced  Stark  broadening 
parameters  S(a)  for  the  Badaer  lines  are  tabulated  for  various  electron  densities 
and  temperatures.  The  S(a)  corresponding  to  interaediate  electron  densities 
were  found  by  log-log  Interpolation.  A  correction  for  temperature  was  not  re- 
ouired  because  the  temperature  dependence  of  these  paraaeters  is  weak,  and  the 
S(a)  for  10,000  K  and  20,000  K  are  given  by  Gldea. 
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starting  with  an  asstmed  we  calculate  the  line  width  w  *  2  AX  corres¬ 
ponding  to  the  value  for  a  given  in  Griem*s  tables,  and  measure  the  relative 
intensity  V(Ax)  on  the  experimental  line  nrofile.  For  each  assumed  value  n 

e 

for  the  electron  density  we  thus  obtain  a  set  of  parameters  , 

and  by  iterating  over  n  we  find  the  best  fit  between  the  measured  line  shape 

e 


and  the  Stark  profile,  according  to  the  condition* 

(s/v)  =  const  =“  s/v  «  F 

'  n  o' 


(5-4) 


The  index  n  refers  to  the  position  on  the  line  profile,  aind  is  conveniently 
expressed  as  the  (apparent)  fractional  line  width*  n  =  V(A\)  / 


TABLE  III 


COMPARISON  OF  MEASURED  LIbffl  SHAPES  WITH  THEORETICAL  STARK  PROFILES 

a)  Balmer  H, 


a)  n  (assumed) 

0 

=•  7  X  10 

14  -3 

1  cm 

j  F^A  » 

0.74 

V/o 

0.19  0,39 

0.53 

0.71 

0.78 

0.88 

1 

s/v 

0.63  0.64 

0.66 

0.67 

0.70 

0.73 

0.51 

b)  (assumed) 

=*  8  X  10 

14  -3 

cm  ^ 

V 

U 

0.80 

v/v 

M  O 

0.17  0.37 

0,51 

0.68 

0.78 

0.89 

1 

s/v 

0.71  0.70 

0.71 

0.71 

0.75 

0.75 

0.51 

c)  n^  (assumed) 

=  1  X  10 

cm 

5  F^A  = 

0.94 

0.28  0.43 

0,61 

0.71 

0.77 

1 

S/v 

0.91  0.8i 

0.78 

Balmer  ' 

0.79 

H 

Y 

0.58 

0.51 

n^  (assumed)  = 

7.5  X  10-"^  cm  ' 

’  >  V*  * 

0.48 

v/v 

0.07  0.14 

0.30 

0.38 

0.43 

0.52 

0.60 

s/v 

0.44  O.JWj. 

0.44 

0.44 

0.46 

0.44 

0.42 
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The  results  of  the  line-shape  measurements  corresponding  to  the  arc  at 

30  Amps,  at  reduced  pressure,  are  presented  in  Table  III,  In  the  case  of  H  , 

5 

the  agreement  of  the  S/V  values  for  small  n,  corresponding  to  the  line  wings, 

14  -3 

clearly  favors  an  electron  density  around  8  x  10  cm  ^  over  values  which  are 

about  20%  higher  or  lower.  The  shape  of  the  wings  of  the  concurrently  measured 

14  -3 

H  line  is  consistent  with  a  theoretical  Stark  profile  for  n  «  7.5  x  10  cm 
Y  e 

The  narrow  peak  of  this  line  is  depressed  due  to  self -absorption  and  distorted 
by  instrumental  broadening,  leading  to  a  value  for  S^/V^  that  is  too  low. 

The  opposite  effect  occurs  for  the  line  whose  narrow  central  dip  predicted 
by  theory  is  smeared  out  due  to  instrumental  broadening. 

These  results  could  not  be  verified  by  an  independent  determination  of  the 
electron  density.  In  our  measurements  of  the  center  region  of  the  low -current 
arc  at  atmospheric  pressure,  however,  the  intensity  of  the  continuum  could  be 
used  for  such  a  check.  In  this  plasma,  the  concentration  of  hydrogen  and  oxygen 
atoms  is  too  low  to  be  of  significance,  and  the  continuum  radiation  is  then  prac¬ 
tically  identical  to  that  of  a  pure  argon  plasma  at  the  same  temperature.  The 
spectral  radiance  of  the  continuum  at  wavelength  is  thus  represented  by  a  simple 
formula 

(cent.)  «  G 

G  is  a  known  constant,  the  ^-factor  -  which  is  only  weakly  temperature  dependent  - 

accounts  for  depart,ures  from  hydrogenic  behavior  and  quantum-mechanical  (Gaunt) 

corrections.  Based  on  the  ^-factors  deterrdned  by  Wende  (30),  and  with  T^..  9300K 

1 S  —3 

we  obtain  n^  *  6.3  x  10  ^  cm  from  the  measured  continuum  intensity  at  two 
o  o 

wavelengths!  3090  A  and  4370  A. 


”e  ^a.T) 
^2 


(5-5) 
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In  this  plasaa  the  higher  excited  states  of  argon  are  in  equilibrim  with 
ions  and  free  electrons.  The  Saha  equation  can  then  be  used  to  deteraine  the 
electron  density  from  measured  line  intensities,  since  the  concentration  of 
argon  ions  and  of  electrons  are  practically  equal.  Thus  n*(j),  the  measured 

O 

number  density  of  the  upper  state  for  Ar  I  43OO  A  yielded 1  n 
The  electron  density  determined  from  the  shape  of  and  in  this  experiment 
was:  n  ="  (6,0  +0,5)  10  ^  cm  Within  the  probable  experimental  error,  all 
three  methods  yield  essentially  the  same  electron  density, 

C,  Measured  Population  Distributions  over  Atomic  Energy  Levels 

In  addition  to  the  Balmer  series  of  hydrogen  from  to  (the  higher 
members  were  too  weak  for  precise  measurements),  several  0  I  lines  were  meas¬ 
ured,  The  transition  probabilities  were  taken  from  the  critical  compilation 
by  Glennon  and  Wiese  (^),  Concurrent  emission  and  absorption  measurements 
of  a  line  yielded  the  population  distribution  between  the  two  energy  levels 
involved^  characterized  hy  the  distribution  temperature  (see  Section  2k), 
However^  most  of  these  measurements  were  energy  limited^  and  in  the  case  of 
the  oxygen  lines  wider  slits  had  to  he  used  than  permissible  for  a  measure¬ 
ment  of  the  true  line  shape.  Only  when  the  measured  absorptance  was  below^ 
say^  0.30  did  we  derive  the  integrated  line  radiance  by  correcting  the 
emission  line  for  self-absorption  in  order  to  obtain  n(j)^  the  occupation 
of  the  upper  energy  level.  It  is  then  preferable  to  use  multiplets  for  a 
concurrent  measurement  of  n(j)  and  n(i).  The  line  strength  ratio  of  the 
multiplet  components  is  accurately  known  from  theory^  and  a  weak  line  in 
the  multiplet  can  then  he  used  to  obtain  n(j)_,  and  a  strong  line  -  for 
which  a  is  significantly  higher  -  to  derive  T  leading  to  n(i). 
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a)  The  Low-Current  Are  at  Atmos|^eric  Pressure!  Ihe  ahsorptances  of  the 


lines  belonging  to  the  oxygen  triplet  at  7773  A  (multiplet  #  1  in  Moore’s  no¬ 
tation,  Ref.  32)  were  of  a  convenient  magnitude  (a  “  0.57  for  the  strongest, 
and  a  *  0,28  for  the  weakest  line)  for  a  concurrent  meastirement  of  n(j)  and  n(i). 
In  the  peripheral  zone  of  the  arc,  the  absorptance  of  was  high  enough  for  a 
determination  of  the  occupation  of  the  first  excited  state  of  hydrogen.  The 
other  measured  B2LLmer  lines  were  practically  optically  thin. 

The  results  for  hydrogen  and  oxygen  are  presented  in  Pig,  13  in  the  form 
of  a  Boltzmann  plot.  The  downward  curvature  of  the  plots  corresponding  to  the 
higher  energy  levels  in  the  off-axis  meaisurements ,  is  a  clear  indication  of 
departures  from  equilibrium. 

In  order  to  determine  the  occupation  of  excited  Ar  I  levels,  we  measured 

o 

the  intensity  of  Ar  I  4-300  A,  and  of  three  other  Ar  I  lines  in  the  same  wave¬ 
length  region.  The  differences  in  energy  terms  are  so  small  that  these  addi¬ 
tional  lines  are  only  useful  for  a  check  of  internal  consistency.  Based  on 
Wende's  value  of  3.1  *  10^  s”^  for  the  transition  probability  of  Ar  I  4300  (33). 
and  using  T  =■  9300  K  for  the  small  temperature  correction,  we  can  express  the 
results  of  these  intensity  measurements  at  the  arc  center  in  the  form  of  a  set 
of  calculated  transition  probabilities  for  these  linesi  Ar  I  4333*  ^.9  (4,8)  } 

Ar  I  4335j  2.9  (2.9)  *  and  Ar  I  4345*  2.3  (2,2)  x  10^  s“^.  Ihe  agreement  with 

Wende’s  values,  indicated  in  parentheses,  is  excellent. 

In  the  temperature  range  of  interest,  the  fractional  ionization  of  argon, 
oxygen,  and  hydrogen,  is  almost  identical i 

[Ar'^J  /  [Arj  =■  F(Tg)  /  n^  «  [o'^]  /  [o]  «  [h'*’]  /  [hJ  (5-6) 
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The  equilibrium  exposition  of  an  argon  plasma  is  thus  not  disturbed  by  the 
additions  of  small  amounts  of  water  vapor.  As  seen  previously,  in  our  measure¬ 
ments  of  the  arc  center  three  independent  methods  had  yielded  essentially  the 
same  electron  density,  which  is  practically  equal  to  the  concentration  of  argon 

ions.  The  occupation  of  a  highly  excited  state  of  the  ion  was  determined  by  a 

o 

measurement  of  the  absolute  intensity  of  the  strong  Ar  II  4348  A  line.  The 

recorded  line  was  barely  above  the  noise  level,  but  in  this  temperature  range 

its  intensity  is  a  very  steep  function  of  temperature  and  relatively  large  errors 

in  the  ^measured  intensity  do  not  affect  significemtly  the  derived  temperature. 

We  thus  obtained  -  9300  K  for  the  population  distribution  between  the  ion 

ground  state  and  the  upper  state  of  this  line.  This  is  close  to  the  temperature 

of  9000  K  obtained  from  the  Boltzmann  plot  for  the  population  distribution  over 

the  hydrogen  and  oxygen  energy  levels  (Fig,  13), 

Ihe  population  density  of  the  argon  ground  state  is  practically  equal  to 

the  total  number  of  particles  at  P  «  1  atm  and  T  «  5OGO  K  (see  Section  5E), 

as  given  by  the  ideal  gas  law.  Thus  one  derives  a  distribution  temperature  of 

8800  K  with  respect  to  the  population  distribution  between  the  ground  state 

o 

and  the  upper  energy  level  for  the  Ar  I  4300  A  line. 

Based  on  these  measurements  it  can  be  concluded  that  in  the  center  of  the 
arc  the  excited  states  of  the  atoms  are  in  Saha  equilibrium  with  their  ions  and 
the  free  electrons  at  an  electron  temperature  of  about  93OO  K,  If  small  depar¬ 
tures  occur,  they  are  not  significantly  larger  than  the  probable  measurement  error. 
The  ground  state  of  argon  appears  to  be  slightly  overpopulated,  however. 

The  electron  density  in  the  plasma  region  corresponding  to  the  arc  perij^ery 

was  determined  by  the  method  of  Section  5B,  using  H  as  the  diagnostic  line.  A 

Y 

good  fit  was  obtained  for  the  wings  of  the  measured  line  shape  and  the  theoretical 
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In  order 


14  -3 

Stark  profile  corresponding  to  an  electron  density  of  2.5  10  cm  . 

to  obtain  the  distribution  temperature  over  highly  excited  states  of  the  argon 
atom,  we  calculate  first  via  the  Saha  eq^uation  the  state  population  n*(oo)  corres¬ 
ponding  to  the  energy  level  of  the  atom  which  coincides  with  the  reduced  ioniza¬ 
tion  energy  Section  7)  write  the  Saha  equation,  with  [Ar+]  =  rig: 

log  “  2  log  n  -  log  Q.(Ar+)  -  1.5  log  Tg  -  15.684  (5-7) 


Eq,  5”?  was  evalaated  with  T  «  78OO  K,  The  temperature  dependence  of  the 
partition  function  for  argon  ions  is  so  wesdc,  that  Q(Ar‘^)  can  he  taken  as  constant 
over  the  relevant  temperature  range.  We  thus  obtain  a  distribution  temperature 

=  3200  K  for  the  population  distribution  between  the  energy  levels  joining 
onto  the  continum  and  the  upper  state  of  the  Ar  I  4300  A  line^  whereas  the 
population  distribution  between  the  latter  and  the  argon  ground  state  is 
characterized  by  a  distribution  temperature  of  7600  K.  The  distribution 
over  the  argon  energy  levels  follows  therefore  the  pattern  found  for  hydrogen 
and  oxygen  (see  Fig.  13). 

The  magnitude  of  the  distribution  temperature  characterizing  the  population 
distribution  over  the  higher  excited  atomic  states  is  of  crucial  importance 
for  the  interpretation  of  the  mechanism  bringing  about  the  noted  departures 
from  LTE.  Even  assuming  an  error  of  30^  in  the  determination  of  the 
electron  density  -  which  appears  unduly  pessimistic  -  the  distribution  tem¬ 
perature  found  in  the  case  of  argon  would  still  be  below  3600  K.  The  measured 
population  distribution  between  the  4th  and  5th  energy  level  of  hydrogen 
corresponds  to  a  distribution  temperature  of  2900  K.  These  distribution  tem¬ 
peratures  are  definitely  lower  than  the  kinetic  temperature  of  heavy  particles 
(see  Section  9C) . 
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1d)  The  High-Ciirrent  Arc  at  Reduced  Press-ure  (Arc  core  at  30  Amps,  and  22  Torr) : 
The  ahsorptance  of  the  first  two  Balmer  lines  was  too  high  for  a  direct  determination 
of  n(3)  and  n(4),  hut  the  distribution  temperatures  T£,(2-3)  and  Tp(2--4)  were 
obtained  with  good  precision  by  the  emission-absorption  method.  The  meas-ured 
ahsorptance  of  Hy-  was  of  a  convenient  magnitude  to  determine  the  occupation  of 
both,  the  upper  and  the  lower  energy  level  for  the  transition.  With  n(2)  thus 
measured,  the  previously  determined  distribution  temperatures  yielded  the  values 
for  log  n*(3)  and  log  n*(4)  shown  in  parentheses  in  Table  IV.  The  ahsorptance  of 
Hg  could  not  be  measured  accurately,  but  was  certainly  less  than  0.20.  The  cor¬ 
rection  for  self-absorption  of  the  integrated  line  intensity  amoimts  then  only  to 
a  few  percent,  and  this  line  can  be  considered  optically  thin. 


TABLE  IV 

THE  BALMER  HYDROGEN  LINES 


line  \  Principal  Q,N.  log  n*(j)  a  log  n*(i) 

O 

(a)  upper(j)  l0¥er(i)  (k) 


6563 

3 

2 

(10.73) 

1  6400 

4861 

4 

2 

(  9.80) 

0.81  5300 

4340 

5 

2 

9.23 

0.32  4800 

4102 

6 

2 

8.83 

The  gas  temperature  T  5700  K  was  derived  from  a  measurement  of  rotational 
line  intensities,  and  we  assume  that  the  chemical  composition  of  the  plasma  is 
given,  effectively,  by  2  H  +  0,  The  occupation  of  the  hydrogen  ground  state 
is  then  obtained  from  the  ideal  gas  law,  since  the  population  of  excited  states 
is  insignificant.  The  contribution  of  other  species  (e.g.  OH  and  0^)  to  the 
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plasma  composition  can  certainly  be  disregarded  in  this  context,  but  it  is  quite 
possible  that  the  ratio  of  hydrogen  to  oxygen  atoms  is  less  than  2:1  in  the  arc 
center. 

The  measured  number  densities  sore  shown  in  Fig.  14  in  the  form  of  a  Boltzmann 
plot*  It  is  readily  seen  that  the  population  distribution  over  the  hydrogen 
energy  levels  is  far  from  equilibrium. 

The  absorption  measurements  of  the  oxygen  atoms  (Table  V)  shoved  the  three 
lines  belonging  to  multiplet  #  1  of  0  I  to  be  black  (a  =  l) .  This  followed 
already  from  the  emission  measurements:  The  three  lines  were  of  the  same  in¬ 
tensity  although  their  line  strength  ratio  is  1:1.67:2.33.  Since  E.(^l)  = 
Ej_(#10),  and  the  occupation  of  the  latter  was  obtained  by  an  emission- 
absorption  measurement  of  multiplet  #  10^  we  could  derive  n(i)  of^l. 


TABLE  V 

THE  OXYGEN  LINES 
Multiplet  upper  term  lower  term 

#  (22)  X(a)  Ej(cm  Pg.  E^(cin“^)  Pg,  log  n*(j)  a  Tp(K)  log  n*(i) 


a)  Quintets  (Spin  =•  2) 


10 

6158 

102865 

86631 

1 1  1 

8.62 

0.198 

5370 

10.51 

1 

7772 

86631 

S  73768  ^s 

«  1 

4800 

12.18 

3 

3947 

99095 

73768 

9.11 

b)  Triplets  (Spin  =■  l) 

5 

4368 

99680 

76795 

9.59 

35 

7947 

113714 

101135 

8.24 

0.30 

5700 

9.64 

55 

7476 

127282 

^P  113910  ^P 

7.46 

<0.10 

high 

c)  Singlets  (Spin  =«  0) 

38 

7157 

116631 

102662 

8.65 

0.53 

5650 

10.19 

65 

6653 

130943 

^P  115918  ^P 

7.94 

<0.10 

high 

(Pg.  “  parentage  of  term) 
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The  Boltzmann  plot  for  oxygen  is  shown  in  Fig,  15 .  It  should  be  noted 
that  the  high  distribution  temperature  for  tho® energy  levels  of  0  I  that  lie 
above  the  lowest  ionization  potential  [E(j)  >  is  also  borne  out  by  the 

absence  of  measurable  absorption  for  multiplets  #  55  and  ^  65* 

B,  High-Resolution  Measurements  of  the  ft^-6  (0»0)  Line  and  its  Satellite 
As  noted  before,  the  available  background  sources  were  unsuitable  for 
absorption  measurements  of  the  rotational  lines  of  OH,  and  we  resorted  to  the 
”two  line  method**,  whose  theoretical  basis  has  been  presented  in  Section  2A, 

The  0^-6  (0,0)  line  was  chosen  because  it  is  one  of  the  strongest  lines  in  the 
ultraviolet  spectrum  of  GH  and  because  the  wavelength  sei>aratlon  between  main 
and  satellite  line  and  their  line  strength  ratio  are  both  of  a  convenient  magni¬ 
tude.  The  calculated  strength  ratio  of  12.6  was  experimentally  verified  in 
run  E  (see  Table  VIII),  where  self -absorption  was  negligible. 

In  order  to  achieve  the  required  high  spectral  resolution,  some  compromise 
had  to  be  made  concerning  the  spatial  resolution.  An  aperture  of  f/l20  was 
chosen  as  compared  to  f/300  used  in  the  ** pinpoint"  measurements  described  in 
Section  4,  In  these  measurements  of  the  arc  at  reduced  pressure,  the  regions 
of  observation  corresponded  to  arc  zones  of  relative  uniformity  -  at  least  with 
respect  to  the  intensity  of  the  rota&onal  lines  -  by  focusing  on  the  arc 
center  (r  »  O),  and  on  the  location  of  the  off-axis  peak  (r  =■  1.1  mm)  shown  in 
Fig,  11,  We  also  made  measurements  close  to  the  arc  periphery  (r  =  1.7  mm). 
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although  in  this  case  only  average  values  over  an  inhomogeneous  region  were  ob¬ 
tained. 

At  the  narrow  slit-widths  required  for  these  high-resolution  measurements, 
the  calibration  signal  from  the  strip  lamp  is  so  weak;,  that  it  is  completely 
buried  in  noise.  Hence  advantage  is  taken  of  the  fact  that  the  meaisured  inte¬ 
grated  line  intensity  of  an  isolated  line  is  independent  of  the  slit  function. 

Thus  we  used  of  the  Qj_-6  line  as  determined  with  wide  slits  (lOOp  for 

the  entrance^  and  400p  for  the  exit  slit)^  in  turn^  for  the  radiance  calibra¬ 
tion  of  the  narrow  slit  (25/25p)  measurements.  By  focusing  on  the  arc  center 
we  ensured  that  the  slightly  increased  field  of  view  in  the  calibration 
measurements  -  determined  by  the  wider  entrance  slit  -  still  corresponded  to 
a  relatively  uniform  arc  zone. 

The  spectral  resolution  was  determined  experimentsdly  by  a  concurrent  scan 

o 

of  the  hyperfine  structure  of  Hg  I  31^5  A  under  identiceil  conditions  as  to  aperture 

and  slit  widths.  In  order  to  improve  the  signal- to -noise  ratio_,  the  slit  height 

was  increased  to  5  this  can  affect  the  resolution  only  adversely.  The 

attained  resolution,  as  determined  by  the  half  widths  of  the  hyperfine  components, 

o 

closely  checked  the  calculated  resolution  of  8  mA  (i,e,  a  resolving  power  of  al^ 
most  400,000), 

In  Fig.  16  we  show  the  spectrum  recorded  in  scan  #  2  of  Table  VI,  The  cal¬ 
culated  Gaussian  of  the  same  half  width hw*  and  normalized  at  the  line  peak,  fits 

the  recorded  line  shape  of  the  Q^-6  line  very  well,  even  out  into  the  wings.  Hence 
the  contribution  of  a  Lorentzian  to  the  line  shape  has  to  be  small,  and  it  can  be 
concluded  that  pressiire  broadening  (including  Stark  broadening  by  chairged  par¬ 
ticles)  is  legligible  as  compared  to  Doppler  broadening.  In  all  experiments, 


including  the  measurements  at  atmospheric  pressure,  the  recorded  profiles  of 

this  line  were  adways  true  Gaussians,  This  result  cannot  he  simply  generalized 

for  all  rotational  lines  of  this  band  system,  however*  Pressure  broadening  may 

be  more  important  for  the  higher  rotational  lines  (34) 

Thus  the  only  broadening  mechanisms  that  have  to  be  considered  in  the 

case  of  the  Q^-6  line  are  Doppler  and  instrument  broadening.  However,  since 

hw*  is  about  5  times  larger  than  the  spectral  slit  widths  the  correction  for 

.  *) 

instrument  broadening  is  only  of  the  order  of  5^  \  i.e.  of  the  same  order  as 

the  estimated  measurement  error*  The  correction  for  self*-absorption  is  con¬ 
siderably  more  important.  First  the  peak  absorptance  was  obtained  from  the 
measured  peak  intensity  ratio  /  I”  ,  and  then  was  used  to  determine  the 
fractional  line  height  corresponding  to  the  true  half  width  hw  (see  Section  2A), 
The  results  of  our  measurements  are  summarized  in  Table  VI. 

It  can  be  seen  that  hw ,  the  true  half  width  of  the  unabsorbed  line , 

as  obtained  from  the  recorded  line  shape,  practically  equals  hw*  of  the  satellite 
line  (which  is  not  self-absorbed,  and  therefore  hw*  =»  hw).  Since  the  weaker 
line  is  noisier,  hw  used  to  calculate  the  Doppler  temperature,  i.e, 

the  kinetic  temperature  of  heavy  particles  or  gas  temperature  T 

The  source  function  -  derived  from  the  measured  peak  intensity  and  the 

calculated  of  the  0^-6  line  -  yielded  Tp(exc),  the  distribution  temperature 

2  2 
for  the  A  S  excited  electronic  state  and  the  X  II  ground  state  of  OH, 


*)  The  square  of  the  half  width  of  a  convolution  of  two  Gaussians  equads 
the  sum  of  the  squares  of  the  half  width  of  the  components* 
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TABLE  VI 


RESULTS  OP  HIGH-RESOLUTION  MEASUREMENTS  OP  THE  (^-6(0,0)  LINE  AND  ITS  SATELLITE 


Sean  # 

1 

2 

3 

4 

5 

Arc  current  (Abp.) 

5© 

50 

50 

15 

15 

Flasna  Pressure  (Terr) 

20 

20 

20 

18 

18 

Radial  are  zone  (an  off  center) 

0 

1.1 

1.7 

0 

1.7 

/"o 

11.55 

9.7 

7.65 

8.1 

6.1 

Peak  absorptance 

0.18 

0.43 

0.68 

0.63 

0.82 

Distribution  Teaperature  Tjj(exc)  (k) 

10,500 

9,700 

8,200 

9,250 

6,900 

0^-6 t  Apparent  halfuidth  ha* 

(■A) 

k? 

43.5 

44.5 

47 

40.5 

True  halfwidth  ha 

(«A) 

^3.5 

40 

36.5 

39.5 

31 

Q2j^-6t  Apparent  halfaidth  ha* 

(hA) 

45 

40 

37 

41 

31.5 

Dopider  teaperature  T 

s 

(K) 

6,500 

5,500 

4,600 

5,400 

3,300 

(K) 

5,600 

4,900 

4,200 

4,700 

3.750 

1,16 

1.12 

1.09 

1.15 

0.9 

p(0,0)  X  10^ 

("ea") 

0.36 

0.32 

0.30 

0.40 

log  P*  +  log  n*(A,0,0) 

14.6? 

14.89 

14.61 

14.95 

14.41 

log  P®  +  log  n*(X,0,l) 

16.59 

16.97 

17.07 

17.13 

17.34 

log  (l(0H) 

3.27 

3.14 

3.03 

3.14 

2.95 

log  P°  +  log  [OHJ 

19.86 

20.11 

20.10 

20,27 

20.29 
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In  order  to  determine  the  distribution  over  rotational  states,  the  in¬ 
tegrated  intensities  of  the  adjacent  P^-2  and  R^-20  lines  of  the  0,0-band 
were  also  measured  (**wide  slit”  measurements,  see  Section  6a),  Since  the 
(undistorted)  line  shapes  and  the  source  functions  for  all  three  lines  are 
practically  indentical,  the  peak  absorptMce  of  the  line  was  used  to 

correct  these  for  self-absorption,  yielding  (i.e.  the  integrated  radiance 
of  the  optically  thin  line).  We  thus  obtained  an  “apparent  rotational  Boltz¬ 
mann  plot  "  for  each  of  the  five  scans,  Ibese  plots  are  strictly  linear, 
which  is  not  at  all  the  case  if  the  uncorrected  are  used.  The  consistency 
between  the  derived  Nj^(P^-2)  and  indicative  of  the  precision  of 

our  method  of  correcting  for  self  absorption.  The  correction  factors  covered 
a  relatively  wide  ranges  thus  for  scan  #  2t  Q^-6:  =»  1,19  =»  1,05 

and  for  scan  #  5*  =»  1.51  j  P^-2*  =•  1.12 

In  all  scans  except  #  5t  the  gas  temperature  is  higher  than  T*  .  derived 

rot 

from  the  “apparent  Boltzmann  plot”.  This  result  can  clearly  not  be  attributed 
to  an  underestimate  of  instrumental  broadening,  because  then  one  would  expect 
the  ratio  T^/t^^^  to  increase  with  decreasing  hw,  whereas  the  opposite  is 
found. 

As  discussed  in  Section  6D,  the  relative  strengths  of  the  rotational  lines 
of  a  band  are  not  simply  given  by  the  Hoenl-London  factors  -  which  were  used 
here  for  the  apparent  Boltzmann  plot  -  but  probably  show  also  a  slight  exx>onential 
dependence  on  the  rotational  energy  term.  It  follows  that  the  difference  between 
the  reciprocals  of  the  "apparent  rotational  temperature”  and  the  gas  temperature 
should  be  a  constant  (Eq,  6-6),  In  effect,  when  we  evaluate  p(0,0)  «  l/kT^^^  -  l/kT 
we  obtain  for  the  first  four  scans  an  approximately  constant  value  of 
p(0,0)  «  0.36  X  10  ^  ”cin”  (This  unit  appears  because  the  energy  term  is  expressed 
in  cm  ^) 
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The  diverging  result  for  scan  #  5  is  attributed  to  the  spurious  effect 

of  plasma  inhomogeneity.  This  measurement  corresponds  to  a  region  of  relatively 

steep  radial  temperature  gradients,  amd  in  the  averaging  involved  in  the  wide- 

slit  measurements  used  to  derive  T*  . ,  the  regions  of  higher  temperature  will 

rox 

be  more  strongly  weighted*  The  probable  error  in  the  Doppler  temperature  is 
also  largest  for  scan  #  5  because  of  the  large  correction  for  self -absorption* 
the  higher  level  of  experimental  noise,  and  the  greater  sensitivity  of  T 
(which  is  proportional  to  the  square  of  the  half  width)  to  small  measurement 
errors* 

By  extrapolating  the  Boltzmann  plot  to  K  «  0,  we  obtain  n*(A,0,G),  i*e* 

2 

the  occupation  of  the  vibrational -rotational  ground  state  of  the  A  S  electronic 

state  of  OH,  times  the  molecular  constant  called  here  and  defined  by  Eq*  2-31. 

The  measured  distribution  temperature  Tj^(exc) yields  then  P^  n*(X,0,l).  The 

occupation  of  the  ground  state  of  OH  (lowest  vibrational -rotational  level  of 
2 

the  XU  electronic  ground  state)  is  related  to  [[OHJ,  the  total  concentration 
of  OH  in  the  plasma,  via  the  partition  function* 

Because  of  the  large  gap  between  the  energy  terms  for  the  X  and  A  electronic 
states,  the  contribution  of  the  A  state  to  the  partition  function  is  negligible 
even  at  the  highest  Tj^(exc)  measured  here*  The  population  distribution  over  the 
vibrational -rotational  levels  of  the  electronic  ground  state  is  determined  by 
the  gas  temperature,  and  T  was  used  to  evaluate  Q(0H).  In  run  #  5  we  estimate 

o 

T  «  4100  K. 

E 

A  shown  in  Table  VI,  the  concentration  of  OH  is  practically  constant  in 
scans  #2  to  #5,  and  only  slightly  lower  in  the  center  of  the  high-current  arc* 
This  important,  albeit  rather  unexpected  result,  clearly  demonstates  that  the 
concentration  of  OH  in  the  plasma  is  not  determined  by  the  local  values  of  T 
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or  Tj^(exc),  and  can  only  be  understood  when  the  dissociation-recombination  pro¬ 
cess  in  the  plasma  is  analyzed  on  the  microscopic  level  (Sections  9A  and  lOB), 

With  IPg  «  ^.52  (see  Section  2B),  the  measured  log  [OHJ  is  then  of  the 
order  of  15.6  ,  For  an  equilibrium  water  vapor  plasma  at  the  same  pressure, 

this  concentration  of  OH  corresponds  to  a  temperature  of  about  3300  K  (compare 
with  Figo  4  showing  the  equilibrium  composition  of  a  water  vapor  plasma  at 
slightly  lower  pressure) . 

E,  Concentration  of  OH  in  the  Arc  at  Atmospheric  Pressure 

The  p(0,0)  correction  factor  obtained  in  the  preceding  Section  was  used 
to  obtain  the  gas  temperature  from  the  "apparent  rotational  temperature"  for 
the  0,0-band,  Even  in  the  peripheral  zone,  the  noise  level  was  too  high  for 
precision  measurements  at  high  spectral  resolution.  Based  on  the  two  line 
intensity  ratio,  the  peak  absorptance  of  the  02^*6  line  appears  to  be  less  than  G*10; 
T^(exc)  is  therefore  certainly  higher  than  T  *■  5OOO  K,  but  its  exact  magnitude 
could  not  be  determined.  The  shape  of  the  0^-6  line  was  measured  with  relatively 
wide  slits  (50/50^),  its  profile  is  Gaussian  and  consistent  with  a  Doppler 
temperature  of  the  order  of  50OOK, 

n*(A,0,0)  was  obtained  by  a  measurement  of  the  absolute  intensity  of  the 
o 

2 

P^-2  line.  The  population  density  of  the  X  IE  ground  state  in  the  arc  center  is 
then  estimated  by  setting  T_(exc)  «  9000  K  «  T  ,  since  the  population  distribu- 
tion  over  the  excited  levels  of  atoms  has  been  shown  to  be  close  to  an  equilibrium 
distribution  at  the  electron  temperature.  For  the  plasma  corresponding  to  the 
peripheral  arc  zone  we  calculate  P®  J  with  T^(exc)  *  7800  K  (i,e,  the  distri¬ 
bution  temperature  corresponding  to  the  linear  portion  of  the  atomic  Boltzmann 
plots,  and  presmably  the  electron  temperature),  and  also  with  Tjj(exc)  «  6OOO  K, 

We  then  obtain  the  occupation  of  the  OH  ground  state  by  using  again  log  P^  »  4,52, 
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In  equilibrium,  the  occupation  of  the  ground  state  of  OH,  of  0,  and  of  H 


are  related  by  the  mass  action  equation: 

n*(l)  ^2iTm„kTl6l 

n*  (1)  i  i?  171 


3/2 


exp  (-DykT) 


(5-8) 


where  D  =  35,501  cm“^  is  the  dissociation  energy  of  OH,  The  first  factor  on 
o 

the  right  side  of  Eq,  5-8  contains  the  translational  partition  ftmction,  and 

consequently  for  this  term  T  •  T  ,  The  gas  temperature  also  enters  the  expo- 

S 

netial  terra,  since  dissociation  and  recombination  are  predominantly  induced  by 
heavy  particle  collision  (see  Section  IOA),  Sq,  5“8  caun  be  written  in  the 
following  form,  convenient  for  numerical  coraputationi 


log  n*(X,0,l)  -  log  n*(l)  +  log  n*  -  20.24  -  1.5  log  Tg  +  20,200/Tg  (5-9) 


The  occupation  of  the  ground  state  of  hydrogen  and  oxygen  was  obtained 
by  extrapolating  the  respective  Boltzmann  plots  (Fig,  I3)  to  E  =*  0,  The  internal 
partition  function  •  i,e,  essentially  the  statistical  weight  of  the  ground  state  - 
times ri’^(l)  yields  then  and  [o J,  In  both  arc  regions  the  ratio  of  hydrogen 
to  oxygen  atoms  is  practically  1j1,  and  not  at  all  2tl,  as  expected.  For  argon 
at  1  atm  pressure,  saturated  with  water  vapor  at  room  temperature,  the  ratio 
[h]  :  [Arl  **  2  [0  I  I  =1:  21  is  considembly  higher  thaw  the  measured  ratio 

corresponding  to  the  axe  center.  Both  observations  are  strong  evidence  for  arc 
deraixing  due  to  diffusion. 

By  substituting  the  derived  occupation  of  the  ground  state  of  hydrogen  and 
oxygen  in  Eq,  5-9,  we  determined  the  equilibrium  occupation  of  the  OH  ground  state 
shown  in  Table  VII,  where  the  resulte  of  these  calculations  are  summarized. 
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TABLE  VII 


CONGERTRATION  OF  OH  IN  THE  ARC 

Are  region: 

AT  ATMOSPHERIC 

r  -  0 

PRESSURE 

r  =»  1*1  mm 

T-  (0,0) 

(K) 

4800 

4400 

Gas 

tenperature  T  (k) 

& 

5400 

5000 

log 

P°  + 

0 

log  n*jj(A.0,0) 

12.45 

12.93 

Distribution  temperature  Tj^  ( exc ) 

(k)  9000 

7800  6000 

log 

p®  + 
o 

log  n*jj(X.0,l) 

14.69 

15.51  16.29 

log 

),1)  (log  P®  -  k.5Z) 

10.17 

10.99  11.77 

log 

[0]  - 

log  n*(l)  +  0.96 

16.13 

16.85 

log 

[H]  - 

log  n*(l)  +  G.30 

16.11 

16.73 

1 — 1 

:  [Ar] 

a  [Oj  :  [ArJ 

1  :  100 

1  t  23 

log 

“oh  ^  ' 

8.89 

10.57 

log 

aOH^^»^ 

).i) 

1.28 

0.42  1.20 

In  the  last  row  of  Table  VII  we  show  the  overpopulation  of  the  OH  ground 
state  with  respect  to  an  equilibrium  occupation.  Notwithstanding  the  inherent 
ambiguity  in  extrapolating  the  population  distribution  over  excited  states,  the 
result  of  these  calculations  supports  the  conclusion  that  in  the  axe  at  atmos- 
I^eric  pressure  the  concentration  of  OH  is  also  higher  than  the  equilibrium 
concentration  corresponding  to  the  local  values  for  the  relevant  temperature  - 
the  gas  temperature  -  and  the  concentration  of  the  dissociation  products* 
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6.  MEASUREMENT  OF  THE  ULTRAVIOLET  BANDS  OF  OH 


A,  Experimental 

In  view  of  the  large  niuiiber  of  lines  to  be  measiired,  the  standard  procedure 
of  obtaining  the  integrated  line  intensities  by  planimetering  the  area  under  each 
recorded  line,  would  hare  increased  the  data  acquisition  and  reduction  time  beyond 
tolerable  limits*  Since  the  rotational  lines  are  quite  narrow  and  of  Gaussian 
shape,  the  intensity  integration  can  be  performed  directly  by  the  spectrometer 
without  loss  of  precision.  When  a  narrow  entrance  slit  is  used,  and  the  exit  slit 
is  much  wider  than  the  half  width  of  the  line,  the  integrated  line  intensity  is 
proportional  to  the  amplitude  of  the  recorded  line  trace  (peak  signal  V^).  We 
determined  experimentally  that  this  condition  was  fulfilled  for  an  entrance  slit 
width  of  0,1  mm  and  an  exit  slit  width  of  0,4  mm,  and  this  was  the  setting  used 
for  the  ”Long  Runs”,  By  slowly  scanning  over  various  lines  we  determined  the  slit 
function  at  different  grating  angles  for  the  two  grating  orders  used  (8th  and  9th). 
We  always  found  the  area  under  the  line  to  be  equal  to  V^  x  hw".  The  half  width 
of  the  line  trace  (hw* )  in  chart  units  (cm)  is  a  constant  for  each  slit  setting, 
and  does  not  depend  on  grating  angle  or  order. 

By  comparing  the  peak  signal  measured  at  the  scanning  speed  used  in  the  ”Long 
Riins”  (medium  speed  of  the  fine  grating  drive)  with  a  slow  scan,  we  established 
that  pen  overshoot  does  not  occur.  This  prooved  that  the  response  time  of  the 
detection  system  (meiinly  determined  by  the  time  lag  of  the  recorder  pen)  was  faster 
than  the  time  required  to  scan  over  a  line.  In  order  to  ensure  that  only  light 
emitted  by  an  homogeneous  plasma  saunple  falls  on  the  detector,  the  entrance  slit 
was  masked  to  a  height  of  0.25  m  aod  the  imaging  mirror  (M^  in  Fig,  8)  was  stopped 
down  to  f/l20. 
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The  scan  nas  started  in  the  9th  grating  order  at  the  hand-head  of  the  R^(l,0) 

®  0 
branch  at  2811  A  and  continued  in  the  8th  order  at  3l60  A*  The  strong  excitation 

of  the  Schunann-Riinge  bands  of  0^  precluded  the  neasurenent  of  OH  lines  beyond 
o 

3300  A,  The  weak  (0,l)  band  of  OH  could  thus  not  be  measured. 

The  scans  were  interrupted  at  regular  intervals,  and  selected  monitor  lines 
were  then  scanned  repeatedly  in  order  to  check  for  arc  stability.  The  intensity 
variations  were  typically  of  the  order  of  3%  with  no  notable  difference  between 
short-term  and  long-term  stability.  The  repeatability  was  of  the  sane  order  when 
the  airc  was  restarted  and  the  same  experimental  parameters  were  reproduced. 

With  the  aid  of  Ref.  12  all  lines  in  the  spectnm  were  identified  and  assigned 
to  the  rarious  transitions  belonging  to  fire  bands  of  the  OH (A  2  -  XU)  system. 
Only  those  lines  were  used  for  the  data  reduction  whose  peak-to-peak  separation 

O 

from  an  adjacent  line  is  greater  than  about  0.2  A,  ±.e.,  which  are  free  from  line¬ 
blending.  We  selected  a  total  of  390  lines,  distributed  as  follows i  0,0-band’ 

(135  lines)  »  1,0-band  (95  lines)  t  1,1-band  (73  lines)  j  2,1-band  (57  lines)  ) 
2,2-band  (30  lines).  In  all  scans  the  lines  belonging  to  bands  originating  in 
the  V*  *  3  level  -  the  3*2  €md  3*3  bands  -  were  missing. 

B,  Results  for  the  Hain  Branches 

The  measured  integrated  line  radiances  were  reduced  according  to  Eq,  2-25. 

il. 

was  divided  by  the  appropriate  line  strength  factor  S/k  ,  and  the  logiurithm 
(base  10)  of  this  quantity  was  plotted  vs,  the  energy  term  of  the  upper  rotatio¬ 
nal  level.  The  data  corresponding  to  the  R,  Q,  and  P  branches  of  each  band  were 
plotted  separately,  Srane  representative  Boltzmann  plots  are  shown  in  Figs.  17  -  25. 
The  quantity  represented  on  the  ordinate  -  called  log  y  for  simplicity  -  is 
defined  by  the  relation! 
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log  y  2  log  (n^P)  -  10*624 


(6-1) 


-3  -1  3 

if  n*  is  expressed  in  cm  and  P  in  s  fUi  . 

All  obtained  Boltzmann  plots  show  essentially  the  same  pattern.  Below  the 
onset  of  predissociation  the  data  points  lie  on  a  straight  line,  and  within  the 
experimental  precision  the  log  y  values  for  the  and  R^t  the 
the  P^  and  P^  branches,  respectively,  are  Identical*  In  many  cases  the  scatter 
is  too  small  to  be  shown  in  the  figure,  and  then  we  have  drawn  an  open  circle 
around  the  dot  marking  the  data  point*  For  the  levels  affected  by  predissociation 
a  marked  splitting  occurs  between  the  Boltzmann  plot  corresponding  to  the  two 
spin  components,  and  a  different  symbol  has  then  be  used  for  the  weadter  (i*e, 
underpopulated)  f^^  level  (see  Figs*  17  •  25)#  Predissociation  affects  the  occu¬ 
pation  of  states  and  is  then  better  discussed  in  relation  to  the  interaction  proc-^ 
esses  in  the  plasma  (Section  lOB),  than  in  the  present  context  whore  we  are  con¬ 
cerned  with  the  transition  probabilities  for  the  rotational  lines* 

With  the  possible  exception  of  the  (l,0)  band,  discussed  below,  the  log  y 
values  for  all  three  main  branches  corresponding  to  lines  originating  in  the  same 
rotational  state  (i.e*  also  the  same  spin  state  if  the  level  occupation  is  affected 
by  predissociation),  are  identical  within  the  precision  of  our  measurements*  This 
will  be  seen  by  comparing  Figs*  1?  and  18;  and  Figs  20  and  21,  where  the  same  Boltz- 
raemn  plot  has  been  drawn  through  the  data  points  corresponding  to  the  R  and  to  the 
Q  branches. 

If  a  divergence  between  different  branches  does  exist,  it  should  increase 
with  rotational  quantum  number  (l^,  14),  A  comparison  of  the  data  for  the  higher 
rotational  lines  is  therefore  especially  significant*  Unfortunately,  the  lines 
of  high  K  were  very  weak  because  of  predisssoclation,  and  the  experimental  noise 
is  proportionally  higher.  However,  a  data  analysis  showed  the  scatter  to  be  random 
(see  Fig,  19). 


0^  and  ^9  and 
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Systematic  discrepancies  between  the  Boltzmann  plots  belcn^iner  to  different 
bands  were  found  in  some  cases,  but  could  then  be  traced  to  self -absorption.  In 
run  D  of  Table  VIII  the  strong  Q  and  P  lines  of  the  (0,0)  band  were  significantly 
aiffected  by  self -absorption  (see  also  Section 

As  another  possible  source  of  systematic  errors,  the  variation  of  the  instru¬ 
ment  response  with  wavelength  has  to  be  considered.  The  response  curve  of  the 
detector  is  known  to  be  fairly  flat  over  the  wavelength  region  of  interest,  but 

increased  reflection  and  absorption  losses  can  be  expected  at  wavelengths  shorter 
o 

than  3000  A,  Although  the  optical  elements  are  coated  with  magnesium  fluoride, 
the  total  loss  can  be  substantial,  because  of  the  relatively  large  number  of  opti¬ 
cal  elements  in  the  optical  train. 

Of  greater  significance  yet  is  the  efficiency  loss  of  a  mating  used  in  higher 

orders  with  an  increasing  difference  between  grating  angle  and  blaze  angle,  because 

an  increasing  fraction  of  the  incident  light  is  then  diffracted  into  ether  orders 

o 

and  lost.  Thus  in  the  range  2811  to  2850  A,  the  intensity  calibration  showed  an 

increase  of  R(\)  of  more  than  1%  per  A,  which  is  7  times  higher  than  the  variation 

o 

of  the  response  in  the  range  2970  to  3IOO  A,  This  is  to  be  expected,  since  the 

o 

■blaze  wavelength  for  our  grating  in  the  9th  order  is  close  to  3000  A. 

The  rapid  variation  of  r(X.)  in  conjunction  with  the  previously  described 
difficulties  in  obtaining  an  accurate  intensity  calibration  at  the  shorter  wave¬ 
lengths,  greatly  increases  the  likelihood  of  error  in  our  measurements  of  the 
(1,0)  band  -  especially  the  R-branch  -  as  compared  to  other  bands.  It  must  be 
admitted,  therefore,  that  the  small  divergence  resulting  for  the  Boltzmann  plots 
of  the  different  branches  of  the  (l.O)  band  (see  Figs.  22,  23  and  24)  may  well 
be  only  spurious. 
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Bel OH  the  onset  of  predissoeiation ,  the  Boltzmann  plots  for  the  various 

a 

branches  are  well  represented  by /straight  line,  and  the  results  of  our  measure¬ 
ments  are  then  best  presented  in  analytical  form  by  giving  the  equation  of  this 
line,  obtained  by  a  least  squares  fit  to  the  data  points i 

logy  »  log  y^  -  nF(K')  (6-2) 

Hhere  F(K* )  is  the  rotational  energy  tern  in  en”^,  referred  to  the  rotationless 
state  of  vibrational  level  v' .  In  Table  VIII  the  slope  n  is  also  expressed  as 
an  “apparent  rotational  temperature" 


TABLE  VIII 

RESULTS  OP  MEASUREMENTS  OF  THE  MAIN  BRANCHES 

branch  log  y  m  x  10'^  (cm)  T;ot(K)  Fig.  # 


Run  At  P  "  1  atmt  I  -  4,1  Amps.}  Sample  gasi  argon  saturated  with  Hater  vapor 
at  room  temperature.  Arc  region  in  focus t  1.1  mm  off  the  arc  center. 


(0.0) 

2.304 

0.1415 

4400 

17,  18 

R(1,0) 

1.318 

0,1318 

4700 

Q(1,0) 

1.328 

0.1468 

4250 

(1,1) 

1.738 

0.1705 

3650 

(2,1) 

1.00 

0.177 

3500 

(2,2) 

1.04 

0,242 

2500 

,2  Torrj 

I  =  32.5  Amp, 

»  Sample  gas:  pure 

Hater  vapor. 

.ly  scanned  to  maximize  the  intensity  of 

Q^-6(0,0) 

(0,0) 

4,266 

0.1260 

4950 

R(1,0) 

3.399 

0.1090 

5750 

22 

Q(1,0) 

3.399 

0.1282 

4850 

23 

P(l.O) 

3.373 

0.1241 

5000 

24 

(1.1) 

3.?26 

0.1465 

4250 

(2.1) 

3.12 

0.226 

2750 

25 
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TABLE  VIII  (cont.) 


band 

logy^ 

m  X  10^  T'  . 

rot 

Fig.  # 

Run  G-1: 

P  «  18  Torrj 

I  =*  50  Amos . , 

other  parameters 

as  in  Run 

(0,0) 

4.122 

0.1185 

5250 

19 

(1.1) 

3.62 

0.134 

4650 

20,  21 

(2,2) 

3.04 

0.316 

2000 

Run  G-2i 

Repeat  of  Run 

G-1 

(0,0) 

4.101 

0.1135 

5500 

(1.1) 

3.58? 

0.133 

4700 

(2.2) 

3.05 

0.275 

2300 

Run  D-lj 

P  =  9.6  Torrj 

I  =•  18  Amps. 

,  other 

parameters 

as  in  Run 

(0,0) 

4.295 

0.1295 

4800 

(1.1) 

3.765 

0.1495 

4200 

(2.2) 

2.99 

0.256 

2450 

Run  D-2: 

Repeat  of  Run 

D-1 

(0,0) 

4.304 

0,1291 

4800 

(1.1) 

3.717 

0.141 

4400 

(2.2) 

3.01 

0.276 

2250 

Run  Bt  P  ”  60  Torrj  I  =•  50  Amps, 
Arc  setting  as  in  Run  B 

(0,0)  2.351 

(1.1)  1.839 

(2.2)  1.27 


Sample  gass  mixture  of  argon  and  water  vapor, 

0.108  5800 
0.119  5200 
0.188  3000 


In  run  E  an  attempt  was  made  to  lower  the  concentration  of  emitters  in 
order  to  minimime  self -absorption .  As  shown  in  Table  VIII,  the  intensity  of  the 
rotational  lines  was  of  the  same  order  as  in  the  run  at  atmosjrfieric  pressure 
with  an  argon  matrix  (run  A),  and  about  two  orders  of  magnitude  lower  thM  in 
the  other  runs  with  pure  water  vapor  at  reduced  pressure.  However,  the  compo¬ 
sition  of  the  sample  gas  could  not  be  well  controlled,  and  the  arc  stability 
was  adversely  affected. 

C.  Results  for  the  Satellite  Branches 

In  general,  the  weak  lines  belonging  to  the  R22^,  ,  and  satellite 

branches  are  too  close  to  the  corresT>onding  line  of  the  main  branch  for  an  inten¬ 
sity  measurement  with  the  wide  bandpass  used  in  the  "Long  Runs",  The  few  favor¬ 
able  cases  were  used  for  a  check  of  xhe  theoretical  Hoenl-London  factors  by 
comparing  the  intensities  of  satellite  and  main  line.  Although  the  effect  of 
self-absorption  on  integrated  intensities  is  smaller  than  on  peak  intensities, 
the  strong  lines  will  still  be  depressed.  Only  for  relatively  weak  maiin  lines 
will  the  intensity  ratio  be  practically  equal  to  the  line  strength  ratio.  The 
R^-5  and  R^-20  lines  of  the  (0,0)  band  fulfill  this  condition,  since  their  inten¬ 
sity  is  typically  5  times  lower  tham  the  intensity  of  the  Qj^-6  line. 

In  various  scans  we  measured  the  following  intensity  ratio  between  the 
and  linesj  3.58  -  3.69  -  3.69  -  3.^5  -  3.74  -  3.68  -  3.72  -  3.6o  -  3.70. 

The  average  of  these  measinrements  is  3.65,  which  is  very  close  to  3.63,  the  ratio 
of  the  Hoenl-London  factors  calculated  with  the  coupling  constant  a  ■  -  7.54, 

The  results  for  the  R^-20  (0,0)  line  and  its  satellite  are  shown  in  Table  IX.  In 
this  case  the  measured  intensity  ratio  is  also  in  good  agreement  with  the  line 
strength  ratio  calculated  with  the  correct  a(K  *  20)  »•  -  9.0  (see  Section  2B), 
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TABLE  IX 


CALGULAIED  AND  MEASURED  LINE  STRENGTHS  OF  R^-20(0,0)  AND  ITS  SATELLITE 

Line  Strengths 


RjL-20 

R21-20 

Ratio 

Calculated:  "  <’.55)  80.5 

2.08 

38.7 

(a  =  *  9.0  )  80.8 

2.72 

29.7 

Measured:  ^  ^ 

30.5 

(Run  #  D-2) 

29 

The  S(0,0)  branch  lies  in  a  wavelength  region  where  interference  by  other 
bands  is  small,  and  the  first  12  lines  are  free  from  line  blending.  The  higher 
lines  form  a  bandhead,  but  S-19  could  be  measured  again.  This  line  at  3024.21 
is  not  listed  in  Ref.  12. 

The  Hoenl-London  factors  calculated  with  the  proper  a(K)  taken  from  Fig.  2, 
were  used  to  reduce  the  intensity  data  for  the  S  branch  as  explained  in  the  case 
of  the  main  branches.  However,  the  representation  of  log  y  vs.  F(K*)  does  not 
follow  a  regular  pattern.  The  experimental  results  are  therefore  presented  in 
the  form  of  a  relative  P-f actor,  by  using  the  log  y  plot  for  the  R-branch,  i.e. 
the  smoothed  equivalent  data  obtained  in  the  same  run: 

log  P*(S)  2  logy(S,K*)  -  log  Hr, K')  (6-3) 
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c  < 


TABLE  X 


RELATIVE  P-FAGTORS  FOR  THE  S-BRANGH 
K*  S(K)x8  log  P'{S)  logF'(S) 


Run  A 

Hun  B 

Run  G 

1 

0.75 

0.040 

0.037 

0.027 

P 

0 

'-A 

2 

1.20 

0.005 

0.004 

0.004 

0.004 

3 

1.39 

0.004 

0.014 

0.013 

0.010 

4 

1.44 

0.027 

0.027 

0.035 

0.030 

5 

1.42 

0.024 

0.059 

0.011 

G.03I 

6 

1.38 

0.037 

0.041 

0.024 

0.034 

7 

1.32 

0.014 

0.034 

0.020 

0.033 

8 

1.26 

0.038 

0.057 

0,020 

0.039 

9 

1.20 

0.059 

0.053 

0.054 

0.055 

10 

1.14 

0.084 

0.074 

0.039 

0.066 

11 

1.09 

0.049 

0.064 

0.065 

0.060 

12 

1.05 

0.050 

0.078 

0.094 

0.074 

19 

0.89 

0.028 

-0.010 

0.009 

S(k)  =  Hoenl-London  factor 


D.  Ilie  Transition  Probabilities  for  OH  (A^E  -  ) 

It  follows  from  the  linearity  of  the  apparent  Boltzmann  plot  that  both, 

the  occupation  of  the  non-predissociating  rotational  levels  and  the  P-factors 

*) 

depend  exponentially  on  the  rotational  energy  term# 


n*(v’,K') 

n*(v*,K*-0)  exp  r-F(K*)/kT  ] 

(6-4) 

P(q' ,q"  )  = 

P^(q*,q")  exp  C-p(q',q’')F(K)] 

(6-5) 

*)  The  P-factor  for  the  S(0,0)  branch  appears  to  be  an  exception. 
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For  simplicity  we  use  the  notation  to  designate  a  particular  branch  of  a 

band; 

although  the  P-factors  appear  to  be  branch-dependent  only  for  the  (l,0)  band. 

With  the  experimental  data  expressed  ast 

log  y  2  log  -  a  F(K*)  -  log  y^  -  0.4343  F(K')/kr^^  , 

the  exponential  factor  is  given  t>yi 

p  -  l/kT*  .  -  l/kT  “  m/log  e  -  l/kT  (6-6) 

rot  g  g 

The  gas  temperature  was  determined  independently  in  the  measurements  re¬ 
ported  in  Section  51^,  and  by  using  Eq,  6-6,  we  obtained: 

p(0,0)  =  0,36  X  lo"^  ("cm") 

This  ralue  for  p(0,G)  and  the  measured  m(0,0)  were  used  to  obtain  the  gas  temp¬ 
erature  for  the  "Long  Runs".  The  exponential  factors  for  the  other  bands  were 
then  derived  form  the  measured  slope  of  the  apparent  Boltzmann  plots. 

The  concentration  of  OH  in  the  plasma  could  not  be  determined,  and  it  is 
then  not  possible  to  derive  the  absolute  values  for  the  P-factors.  However, 
we  can  still  obtain  relative  values  by  referring  the  P-factors  to  the  molecular 

constant  called  here  P®  and  defined  by  Eq.  2-31. 

o 

P;(q’.q")  =  Po(q*,q**)/P°  “  exp  [g(v' )/kT^]  (6-?) 

Where  G(v*)  is  the  vibrational  energy  term  referred  to  the  rotationless  state  of 
the  lowest  vibrational  level.  We  assume  here  that  the  occupation  of  the  non- 
predissociating  vibrational  levels  is  also  in  a  Boltzmann  equilibrium  at  the  gas 
temperature.  Hie  results  for  the  "Long  Runs"  are  summarized  in  Table  XI.  The  in¬ 
dicated  average  has  been  obtained  by  weighting  the  various  runs  according  to  their 
long-term  stability  and  the  possible  influence  of  self -absorption. 
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TABLE  XI 

NORMALIZED  TRANSITION  MOMENTS  FOB  THE  ROTATIONAL  LINES  OF  OH  (A^S  -  ) 


Run  T 


log  (1/p;) 

(1,0) 

Q  P 


a)  Measured  Values 

/p;) 

(1,1)  (2,1)  (2,2 


p  X  10  (cm) 

(1,0)  (1,1)  (2,1)  (2,2) 

Q  P 


A  5000  0.611  0.601 


0.191  0.57  0.53  0.04  0.49 


1.03  1.18  2.68 


B  5700  0.537  0.543  0.563  0.210  0.51 
c-1  6100  0.194 

C-2  6400  0.224 

D-1  5500  0.190 

D-2  5500  0.24? 

E  6800  0.236 


Average 


0.550 


0.200 


•0.03  0.41  0.34  0.83  2.66 
0.71 
0.82 
0.82 
0.62 
0.61 

0  0.45  0.34  0.75  2 


branch! 


b)  Comparison  with  Calculated  Vajnes 
Franck- Condon  factor  I  Exponential  correction  factor 


?  1  ) 

P 

X  10^  ( 

cm) 

This 

D*C 

James 

A&L 

This 

James 

A&L 

work 

(12) 

(12) 

(14) 

work 

(12) 

(lit) 

R(0,0) 

P(0,0) 

b 

1 

1 

0.99 

1 

1 

0.36 

0.30 

0.55 

0.48 

R(1,0) 

0.26 

0.47 

0 

-0.73 

0.13 

Q(1,0) 

■  0.28 

0,22 

0.46 

0.45 

0.20 

P(1,0) 

el 

0.25 

0.45 

0.34 

-0.53 

0.26 

In  Table  Xlb  the  results  of  the  present  work  are  compared  with  theoretical 
predictions,  Dieke  and  Crosswhite  have  calculated  Franck-Gondon  factors  neglecting 
the  effect  of  vibration-rotation  interaction  and  the  variation  of  the  electronic 
transition  moment  with  intemuclear  distance  (l2).  Anketell  and  Learner  present 
tables  of  calculated  relative  transition  probabilities  for  several  lines  of  the 
band  system  (14).  Although  not  noted  in  their  paper,  their  data  show  an  almost 
strictly  exponential  dependence  on  the  rotational  energy  term,  and  are  then  better 
represented  by  the  expression  of  Eq,  6-5.  James'  data  for  the  (0,0)  and  (l,0) 
band  (l/j)  show  only  slight  departures  from  an  exponential  dependence  and  can  be 
approximated  by  an  expression  of  the  same  form. 

In  run  D-2  we  measured  the  intensities  of  P2(2,l),  euid  also  of  R2(2,2)  emd 
Q(2,2)  lines  corresponding  to  high  K’ -numbers  (14-18).  Within  the  precision  of 
our  measiirements  of  these  weak  lines,  the  log  y  values  for  the  same  K*  were  identi 
cal  for  the  three  branches.  Since  uredissociation  weakens  different  lines  ori¬ 
ginating  in  the  same  upper  energy  level  by  the  same  amount,  it  must  be  concluded 
that  the  P-factors  -  and  therefore  also  their  dependence  on  the  rotational  energy 
term  -  are  similar .  This  is  contrary  to  Anketell  and  Learner's  predictions  (see 
Table  Xlb), 


7.  RATE  PROCESSES  IN  PLASMAS 


In  electric  discharges,  collisions  with  electrons  is  the  predominant  mecha¬ 
nism  for  electronic  excitation  and  for  ionization.  Provided  the  electrons  have 
a  Maxwell  velocity  distribution,  the  rates  for  a  process  induced  by  electron 
collision  and  its  inverse  are  related  by  the  principle  of  detailed  balancing  (2). 
Thus  with  the  collisional  excitation  rate  for  a  transition  from  energy  level  i 
into  energy  level  j  expressed  in  terms  of  the  rate  coefficients: 

Rgoi(i.j)  =  "(i)  "o 

the  following  relation  holds  rigorously  between  the  rate  coefficient  for  excita* 
tion  and  for  the  inverse  process  (de-excitation  from  level  j  into  level  i),  re¬ 
gardless  of  the  equilibrium  situation  in  the  plasmas 

X(,j,i)  -  n®^(i)  X(i,j)  (7-2) 

The  ratio  of  the  equilibrium  populations  of  level  i  and  j  is  given  by  the  Boltz¬ 
mann  factor  containing  the  electron  temperature  T  s 

exp  [-(Ej  -  E^)/kT]  (7-3) 

n  (i)  g(i) 

Similarly,  for  collisional  ionization  from  level  i,  and  three-body  recombi¬ 
nation  into  level  i  (N^  »  total  concentration  of  ions  of  same  species  and  next 
higher  ionization  stage,  c  «  continuum) s 

R^^l(i,c)  =  n(l)  Hg  S(i,c)  h  i7-i^,5) 

3(i,c)  /  S(c,l)  -  F(i,T)  (7-6) 

F(i,T)  is  the  Saha  equation  connecting  the  equilibrium  population  of  energy 
level  i  and  the  concentrations  of  ions  and  free  electrons,  and  Eq,  7-6  follows 
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from  the  principle  of  detailed  balancing.  Introducing  the  quantities  a(i) 
that  relate  the  occupation  of  energy  level  i  to  its  occupation  in  Saha  equi¬ 
librium  with  the  ions  and  free  electrons  (55) , 

a(i)  2  n»(i)  /  n*^*^(i)  ,  (7-7) 

we  have,  as  follows  from  Eqs.  7-4,  7-5,  and  7-6, 

a(i)  =  R  (l,c)  /R  (c,i)  (7-8) 

col  cnl 

In  absence  of  equilibrium  we  find  it  advantageous  to  refer  the  Saha  equa¬ 
tion  to  n*(oo),  the  occupation  of  the  last  bound  state  of  the  atom  whose  energy 
coincides  with  the  reduced  ionization  energy  .  At  the  low  electron  densi¬ 
ties  typical  for  non-LTE  plasmas,  the  lowering  of  the  ionization  energy  is 
negligible,  and  one  can  set  numerically  equal  to  the  ionization  potential 

* 

of  the  free  atom.  The  actual  value  of  i(oo)  does  not  concern  us  in  this  context. 

With  n*(l,l)  «  number  density  of  ions  in  the  ground  state;  Z  *  total  partition 

e 

function  of  free  electrons: 

n  n*(l,l)  /  , 

■  -  -  -  Ze  *  2  (2m^kT/h^y'^  =  4.82  x  10^^  (7-9) 

n*®^(i)  -  n*(o<»)  exp  [(E^  -  E^)/kT]  (7-10) 

The  true  i(«>)  state  is  a  hybrid  because  of  the  mixing  of  the  eigenfunctions 
for  the  bound  and  the  continuous  state  of  the  same  energy,  Clasically,  the  rate 
coefficients  pertaining  to  a  transition  from  i(<»)  into  the  continuum  will  be  very 
large.  Eq.  7-9  will  therefore  hold  fi.e,  n'*(<»)  >■  n*^ (<*»)],  even  when  the  occupa¬ 
tion  of  the  lower  bound  levels  is  not  described  by  Eq,  7-10  ([a(i)  ^  ij,  Tbe 
popiilation  distribution  over  bound  levels  (Boltzmann  plot)  must  then  converge  to 
the  n*(w)  that  is  in  equilibrium  with  ground  state  ions  and  free  electrons. 


*)  For  our  low-nressure  plasma*  «  50  (see  Eq.  6-38  in  Ref.  4) 
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In  principle,  the  rate  coefficients  for  the  various  processes  in  the  plasma 
can  be  obtained  by  averaging  the  appropiate  cross-sections  -  which  depend  on  the 
energy  of  the  colliding  electron  -  over  a  Maxwell  velocity  distribution,.  However, 
with  the  exception  of  hydrogen,  accurate  cross-sections  are  generally  not  avail¬ 
able  because  experimental  data  are  scanty  and  theoretical  calculations  are  beset 
by  considerable  difficulties. 

In  optically  thin  plasmas,  radiative  decay  and  recombination  are  irrevers¬ 
ible  processes  in  the  sense  that  photons  cam  escape  the  plasma,  and  deviations 
from  LTS  can  be  expected  if  the  rate  for  a  collision-induced  process  is  smaller 
or  comparable  to  the  rate  for  the  equivalent  radiative  process,  McWhirter  and 
Hearn  have  obtained  the  population  distribution  over  the  hydrogen  energy  levels 
for  this  collisional -radiative  plasma  model  by  finding  the  steady-state  solutions 
for  the  system  of  rate  equations  over  a  range  of  electron  densities  and  temper¬ 
atures  (25).  Their  calculations  for  an  optically  thin  hydrogen  plasma  with 
*  10^^  cm*"^;  T  *  l6,000  K  -  which  is  comparable  to  the  electron  density  and 
temperature  in  our  plasma  at  reduced  pressure  -  predict  significant  departures 
from  LTE  only  for  the  ground  state,  which  will  be  overpopulated  by  a  factor  of 
the  order  of  ten.  The  higher  energy  levels  are  in  Saha  equilibrium  with  the 
ions  and  free  electrons,  because  the  effectiveness  of  collisional  processes  rela¬ 
tive  to  radiative  processes  is  a  rapidly  increasing  function  of  the  principal 
quantum  number  (See  next  Section), 

Pig,  14  shows  convincingly  that  our  experimental  data  cannot  be  reconciled 
with  these  theoretical  predictions:  The  collisional -radiative  plasma  model  is 
clearly  not  applicable  to  arc  plasmas  at  reduced  pressure. 
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8,  ATOMIC  SPECIES 


A .  Analysis  of  the  Interaction  Processes  in  Hydroi^en  Plasmas 

The  excitation  rate  coefficients  for  hydroj^en  have  been  given  by  Drawin  (36), 
and  are  conveniently  expressed  in  terms  of  the  absorption  oscillator  strengths 
f^.  .  Disregarding  factors  of  the  order  of  unity,  and  with  u.  .  »  (E*  -  S  )/kT} 

®  ionization  potential,  of  hydrogen: 

X(i,j)  -  8.7  X  10"®  (Ej^/kT)^/^  Ejj/(E^-.E^)  exp(-Ujj)  (8-1) 

U  2  lnCl.25(l  +  a"^)]  /  (l  +  (8-2) 

The  collisional  excitation  rate  froa  level  i  into  level  i  +  n,  and  the 
de-excitation  rate  for  the  inverse  process  can  he  written  with 
G(t)  2  1.7  xio"'^  (Ej^/kT)^^^  » 

-  n*(i)  n^  C{T)  I(i.i+n)  exp(-u^^^^)  U  I  (8-3) 

-\o]^(i+n,i)  2  n(i+n)  n^  X(i+n,i) 


n*(i+n)  Hg  G(T)  l(l,i+n)  U  [u^ 


I(i.i+n) 


i,i+n 


(i  +  n)2 


1^  (i  n)^ 

n  (2i  +  n) 


(8-if) 

(8-5) 


®^Pr®ssion  ^iven  hy  Jnsold  ( 57 )  has  heeii  used  to  relate  the  hydrogen  oscilla¬ 
tor  strengths  to  the  principal  quantum  number  i.  c  is  a  slowly  varying  function 
of  i  and  n,  whose  value  ranges  from  1.41  for  to  I.85  in  the  limit  of  large 


We  Kish  to  compare  the  probability  of  a  collisional  transition  from  level  i 
into  the  two  neighboring  levels  i  ±  1  Kith  the  probability  of  a  collisional  tran¬ 
sition  into  more  distant  levels.  In  this  context,  Ke  can  take  c  to  be  effective¬ 
ly  constant  and  also  disregard  the  ratio  of  the  functions  defined  by  Eq.  8-2, 
which  is  close  to  unity.  Thus  we  obtain i 


(8-6) 

(8-7) 


It  follows  from  Eq,  8-?  that  a  collisionally  induced  de-excitation  from  the 
third  into  the  second  energy  level  of  hydrogen  is  about  20G  times  more  probable 
than  a  transition  into  the  ground  state.  Ibis  ratio  decreases  with  increasing 
principal  quanttam  number  of  the  level  involved,  but  even  in  the  limit  of  high  i, 
a  de-excitation  into  the  next  lower  level  will  be  at  least  l6  times  more  proba¬ 
ble  than  a  transition  into  the  level  i  -  2. 


The  ratio  of  the  excitation  rates  into  the  ad.iacent  and  the  next  higher 
level,  given  by  Eq,  8-6,  contains  the  electron  temperature  via  the  exponential 
factor.  For  the  higher  quantum  numbers  this  term  is  insignificant,  and  the  limit 
for  i->a»  is  the  same  as  for  Eq,  8-7,  The  lower  boxmd  for  Eq.  8-6  is  6,  auid 
corresponds  to  the  probability  ratio  for  excitation  from  the  ground  state  to  the 
second  and  to  the  third  energy  level  in  the  limit  of  infinite  temperatures.  At 
T  «  25fOOO  K,  this  ratio  is  larger  than  10  for  all  energy  levels  of  hydrogen. 

Since  transitions  from  level  i  into  levels  higher  than  i  +  2  and  lower  than 


i  -  2  are  still  less  probable,  this  analysis  shows  that  the  error  introduced  by 
considering  only  stepwise  excitation  and  de-excitation  is  never  serious. 
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The  collisional  ionization  coefficients  for  hydrogen  are  given  by  Drawin 

•8  <i6)l 

3(i,c)  -  6.7  X  10“®  (Ejj/kT)^/^  exp  U  (8-8) 

where  U  'tli®  sane  form  as  Eq.  8-2,  with  s  (E^  -  Ej^)/lcT.  The 

ratio  of  collisional  ionization  to  three-body  reconbination  for  each  level  i 
can  then  be  ealetilated  hy  using  Eqs.  7-8  and  7-10, 

Recombination  can  also  occur  by  a  two-body  collision,  in  which  the  atom 
formed  is  stabilized  by  enission  of  a  photon.  The  plasnas  considered  here  are 
optically  thin  with  respect  to  the  recombination  contlnuun,  and  the  inverse 
process  (photo-ionization)  can  surely  be  disregarded.  The  radiative  recombi¬ 
nation  rate  is  given  by  (p  =  free  protons): 

^rad^°'^^  “  "eCpl®(i)  »  (8-9) 

and  one  obtains  for  the  ratio  of  the  probabilities  of  radiative  and  three-body 
reccMBbination  into  level  ij 


.  ,  901 

a(i) 

F(i,T) 

(8-10) 

%  S(i.c) 

am 

a(i) 

2.3  X  10^°  ! 

(8-11) 

"e  [^1 

The  radiative  recombination  coefficients  are  only  weakly  temperature-dependent, 

and  were  obtained  by  extrapolating  to  25,000  K  the  a(i)  tabulated  by  Bates  and 
Delgarno  ( 38 ) . 

Lastly,  we  consider  depopulation  of  level  i  by  radiative  decay  into  lower 
levels: 

«rad^i.:5<i)  -  «  1.5  X  10^0  n*(i)  (8_i2) 
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The  preeciiag  eoneepts  aire  aoe  applied  nuaerleally  te  the  water  Taper 
piMwa.  generated  In  the  eore  ef  the  high  enrreat  are  at  reduced  pressure.  In 
Table  XII  we  eoapare  the  releraat  population  and  depopulation  rates  for  the  first 
six  lerels  ef  hydrogen,  based  on  the  aeasured  electron  density  (n^  •  8  x  10  ca  ■^) 
and  population  distribution  (Pig.  14).  We  also  use  T  25*000  K,  and 
n*(«o)  -  2  X  lo"^  cn“^  (see  Sections  8B  4  D). 

TABUB  XII 

POPULATION  AND  DEPOPULATION  RATES  FOR  HYDROGEN  LEVELS  IN  THE  LOW  PRESSURE  ABC 


4 

s^d.ja) 

1 

(c.'3 

1 

7  X  lO""^ 

15© 

e.ei9 

1,6  X  1©^^ 

2 

5  X  lO"^ 

0.81 

2.7 

©.175 

2,3  X  1©^^ 

3 

-4 

7  X  1© 

3.5  X  1©”^ 

3.5  X  1©“^ 

0.132 

2,7  X  10^ 

4 

4  X  1®"^ 

4.3  X  lo”^ 

2  X  1©'^ 

0.095 

21 

2.8  X  lO'^-^ 

5 

1  X  1©"^ 

8  X  lo"^ 

3  X  10^ 

0.087 

3.5  X  1©^^ 

6 

3  X  lO”^ 

2  X  lO"^ 

5  X  1©’^ 

0.071 

4.9  X  10^^ 

The  ratio  of  the  collisional  ionization  rate  for  atoms  in  energy  level 
i  to  the  collisional  recomhination  rate  into  this  level  equals  a(i)^  the 
overpopulation  of  level  i  (Eq.  7-8) .  As  shown  in  the  second  column  of 
Table  XII^  the  ground  state  is  thus  overpopulated  by  more  than  6  orders  of 
magnitude  with  respect  to  the  occupation  of  this  level  that  is  in  Saha  equili¬ 
brium  with  free  electrons  and  protons.  Three -body  recombination  into  the 
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ground  state  is  insignificant  compared  to  ionization.  Although  radiative  re¬ 
combination  is  larger  by  more  than  two  orders  of  magnitude,  it  still  is  a  very 
infrequent  process  compared  to  ionization.  The  probability  for  three-body  re¬ 
combination  increases  rapidly  with  principal  quantum  number  and  becomes  then 
much  larger  than  the  probability  for  the  radiative  process.  But  even  for  1=6, 
recombination  amounts  to  only  3^  of  the  ionization  rate. 

Radiative  depopulation  is  only  inportant  for  the  first  excited  level  of  hydro¬ 
gen,  and  radiative  decay  from  this  level  into  the  ground  state  is  about  three 
tiroes  as  probable  as  a  collisionally  induced  transition  into  the  next  higher  level. 
However,  Lyaan-a  will  be  practically  completely  absorbed  in  the  plasma,  i.e,  ihoto- 
excitation  will  balance  radiative  decay,  and  the  net  contribution  of  this  larocess 
to  the  depopulation  of  the  i  *  2  level  will  then  be  much  smaller. 

With  the  exception  of  the  levels  i  ■  2,  and  i  «•  3,  the  probability  for  an 
hydrogen  atom  to  bo  excited  into  the  next  higher  level  is  at  least  ten  tines  larger 
than  the  probability  of  a  direct  transition  into  the  continuum.  Although  this  nay 
not  be  entirely  correct  for  the  second  and  third  energy  level,  direct  ionization 
can  therefore  also  be  disregarded  as  compared  to  the  stepwise  excitation  auid  de- 
excitation  process. 

B,  Population  of  Hydrogen  Energy  Levels  by  Stetwise  Sxcitatlon 

The  steady-state  occupation  of  level  i  is  determined  by  equating  the  rates 
for  transitions  into  level  i  to  the  rate  of  transitions  out  of  level  i.  As  seen 
previously,  to  a  good  approximation  we  can  disregard  all  population  and  depopulation 
processes  except  transitions  from  and  into  neighboring  levels  induced  by  collisions 
with  electrons.  The  rate  equation  for  level  i  can  then  be  written j 
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n(i)  Cx(i,l-1)  +  X(i,l+1)]  (8-13) 


n(l-l)  n  X(l-l,i)  +  n(i+l)  X(i+l,l)  - 
o  c 

and,  lay  applying  the  principle  of  detailed  halaneingt 

n(i-l)  [  1  -  J  X(i-l,i)  -  n(i)  [  1  -  J  X(i,i+1)  (8-14) 

The  Quantities  a(j)  were  defined  in  Eq,  7“7#  It  should  he  noted  that  the 
ratio  a(j)/a(k)  is  a  function  of  the  electron  te»perature,  but  that  n*(«>)  cancels. 
The  expressions  in  square  brackets  in  Eq.  8-14,  for  which  we  shall  use  a  short¬ 
hand  notation  -  e.g*  B(i-l,i)  -  denote  the  fraction  of  exciting  collisions  that 
are  effective  in  populating  the  next  higher  level,  taking  into  account  the  inverse 
process  (de-exciting  collisions  into  the  lower  level). 

We  have  used  the  meastired  number  densities  and  the  appropiate  rate  coeffi¬ 
cients  calculated  for  T  *  18,000  ;  20,000  ?  and  25,000  K,  to  obtain  the  state 
occupation  of  the  next  higher  level  and  compare  it  with  the  measured  quantity. 

The  ratios  of  the  various  factors  in  the  expression  for  the  rate  coefficients 
(Eq.  8-1)  were  calculated  separately,  and  are  shown  in  Table  XIII,  Eq.  8-14  is 
then  written  as* 
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TABLE  XIII 


lALCULAlED  OCCUPATION  OF  HYDROGEN  ENERGY  LEVELS  IN  THE  GORE  OF  THE  LOW-PRESSURE  ARC 


log  log  log  log 

i-1  n*(i-l) - - 

■  1(1, 1-H)  e  i-A+1  B(i.l-H) 


log 

n*(i) 


log 

n*(i) 


a) 

Assumed 

electron 

temperature: 

T  =  18.000 

K 

1 

16.12 

-1.522 

-0.888 

-2.327 

0.013 

11.40 

0.82 

2 

12.22 

-0.927 

-0.436 

-0.444 

G.032 

10.44 

0.29 

3 

10.73 

-0.675 

-0.225 

-0.099 

0.081 

9.81 

-0.01 

4 

9.80 

-0.530 

-0.136 

-0.039 

0.088 

9.19 

0.04 

5 

9.23 

-0.  444 

-0.092 

-0.019 

(0.174) 

8.85 

-0.02 

6 

8.83 

8) 

Assumed 

electron 

temperature : 

T  -  20,000 

K 

1 

-0.883 

-2.095 

0.024 

11.64 

0.58 

2 

-0.418 

-0.309 

0.036 

10.57 

0.16 

3 

see 

a) 

-0.218 

-0.089 

0.080 

9.83 

-0.03 

4 

-0.132 

-0.035 

0.088 

9.19 

0.04 

5 

-0.080 

-0.017 

(0.178) 

8.87 

-0.04 

c )  Assumed  electron 

temt>erature  1 

T  =  25,000 

K 

1 

-0.862 

-1.674 

0.029 

12.19 

0.03 

2 

-0.389 

-0.248 

0.041 

10.66 

0.07 

3 

-0.195 

-0.071 

0.085 

9.87 

-0.07 

4 

-0.120 

-0.029 

0.082 

9.20 

0.03 

5 

-0.078 

-0.013 

(0.100) 

8.83 

0.00 

Note t  n*  is  the  measured,  and  n*  the  calculated  state  occupation.  The  occupation 
of  the  energy  level  is  given  by  n(i)  =  2i2n*(i).  Quantities  in  parentheses  are 

based  on  an  estimated  occupation  of  the  energy  level  i  =  7. 
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The  predicted  population  process  is  quantitatively  confirmed  by  the  excellent 
agreement  between  the  measured  and  the  calculated  number  densities,  as  shown  in 
Table  XIII,  It  should  be  noted  that  only  the  relative  occupations  of  the  second 
and  third  energy  level  depend  significantly  on  the  assumed  electron  temperature. 

The  agreement  is  best  for  an  electron  temperature  of  the  order  of  25,000  K,  The 
temperature  can  still  be  higher,  however,  if  the  assumption  of  complete  trapping 
of  resonance  radiation  does  not  hold.  Although  Eq.  8-15  will  then  yield  too  high 
an  occupation  number  for  the  second  energy  level,  this  results  fr<Ma  the  neglect 
of  radiative  decay  from  the  second  energy  level  into  the  ground  state. 

The  ratio  of  the  exponential  factors  approaches  unity  for  the  higher  energy 
levels,  and  the  increase  of  the  factor  containing  the  U-functions  is  approximately 
cancelled  by  the  slow  increase  in  the  ratio  of  de-exciting  to  exciting  collisions 
expressed  Yxf  the  B-factor,  The  state  occupation  pertaining  to  energy  level  i  is 
then  mainly  determined  by  the  I-factor  that  depends  only  on  the  principal  quantum 
number  and  for  i-^  ^  increases  as  i^.  Thus  n*(i)  will  fall  rapidly  with  in¬ 

creasing  principal  quantum  number  i,  leading  to  a  downward  curvature  of  the  Boltz¬ 
mann  plot.  The  failure  to  detect  experimentally  any  measurable  lines  in  the  series 
limit  is  consistent  with  this  pattern, 

C,  Population  Distribution  over  the  Oxygen  Energy  Levels 

The  oxygen  atom  has  a  complex  term  scheme,  and  the  collisional  excitation 
rates  are  not  available;  it  is  therefore  not  possible  to  obtain  the  number  densi¬ 
ties  by  calculation,  as  in  the  simple  case  of  hydrogen.  However,  the  pattern  of 
the  measured  population  distribution,  shown  in  Fig,  15,  can  be  explained  in  a 
qualitative  manner. 
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Gollisional  as  well  as  opticaQ.  transition  probabilities  are  determined  by 
the  matrix  element  of  the  dipole  -  in  general  miiltipole  -  moment.  If  a  tran¬ 
sition  between  two  energy  levels  is  forbidden  for  optical  dipole  radiation,  the 
probability  for  a  transition  induced  by  electron  collision  will  also  be  very 
small.  Intercombination  lines  between  terms  of  different  multiplicity  are  for¬ 
bidden  for  oxygen  (Russel -Saunders  coupling),  and  measured  collisional  cross- 
sections  involving  a  change  of  multiplicity  are  significant  only  at  energies 
close  to  threshold. 

A  splitting  between  the  Boltzmann  plots  for  singlet,  triplet,  and  quintet 

states  can  then  be  expected  in  a  non-equilibrium  plasma.  This  is,  in  fact,  a 

notable  feature  of  the  population  distribution  shown  in  Fig.  15,  Thus  the 
3  —1 

4p"^P  state  at  99092  cm  is  overpopulated  by  a  factor  of  the  order  of  three 
with  respect  to  the  adjacent  4p^P  state  at  99680  cm"^.  The  next  higher  state  - 
3s' at  101135  cm”^  -  lies  practically  on  the  sane  Boltzmann  plot  as  the  4p^P 
state,  2d  though  it  is  a  state  of  different  parentage.  Combinations  between 
states  of  different  parentage  but  same  multiplicity  are  not  only  allowed,  but 
can  be  very  strong. 

Analagously  to  hydrogen,  in  this  plasma  the  ground  state  and  lower  lying 
levels  of  oxygen  are  highly  overpopulated  with  respect  to  the  occupation  that 
is  in  Saha  equilibrixim  with  oxygen  ions  and  free  electrons.  The  triplet  ground 
state  of  0  I  is  populated  by  dissociation  of  the  X^JL  and  predissoeiation  of 
the  A  S  state  of  OH.  The  low-lying  metastable  singlet  state  of  0  I  (2p^  ^D) 
will  be  populated  mainly  by  dissociation  of  oxygen  molecules,  a  substantial 
component  of  the  plasma  (see  Section  IGA).  The  quintet  states  of  01  are  formed 
by  dissociation  of  the  higher  excited  electronic  states  of  0^  and  OH,  which  in 
this  plasma  have  very  small  occupation  numbers.  This  population  mechanism  can 
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lead  to  the  observed  overpopulation  of  the  excited  singlet  states  of  oxygen  with 
respect  to  the  triplet  states,  which,  in  turn,  are  overpopulated  with  respect  to 
the  quintet  states. 

A  shown  in  Pig.  15,  the  two  Boltznann  plots  for  the  doubly  excited  bound 
2  2 

levels  of  oxygen  of  P  and  of*  D  parentage  cross  undisturbed  the  energy  threshold 

A  ^ 

E^(  S)  for  ionization  into  ground  state  ions.  The  measured  distribution  temper- 

A  \ 

ature  is  actually  higher  for  the  energy  levels  above  B^(  S).  It  must  therefore 
be  concluded  that  the  probability  for  an  oxygen  atom  in  a  doubly  excited  state 
to  be  excited  into  a  higher  bound  level  of  same  multiplicity  is  much  higher 
than  the  probability  of  ionization  into  ground  state  ions. 

The  population  distribution  over  bound  states  of  S  parentage  must  converge 
to  the  occupation  of  the  last  bound  state  n*[E^(^S)j  that  is  in  Saha  equilibrium 
with  free  electrons  and  ground  state  ions  (Section  ?).  This  imilies  a  sharp  drop 

A  V 

in  the  Boltzmann  plot  for  the  quintet  states  in  the  energy  interval  between  E^(  S) 

and  the  last  measured  occupation  plotted  in  Fig.  15#  It  also  implies  a  splitting 

a 

between  the  Boltzmann  plot  for  the  higher  excited  triplet  states  of  S  parentzige 

2  2 

and  the  Boltzmann  plot  for  the  triplet  states  of  D  and  P  parentage,  since  the 
point  on  the  latter  corresponding  to  E^(^S)  is  higher  than  n*[E^(^S)J  by  two  orders 
of  magnitade, 

•Hie  splitting  between  Boltzmann  plots  corresponding  to  states  of  allowed 
combination,  can  be  explained  by  the  rapidly  decreasing  probability  for  tran¬ 
sitions  involving  a  quemtum  jump  over  more  than  one  principal  quantum  number.  The 
(^D)3s  ^  state,  whose  number  density  is  plotted  in  Fig,  15,  is  the  only  triplet 

4 

state  of  excited  ion  parentage  below  E^(  s).  1116  next  higher  doubly  excited  tri¬ 

plet  states  of  principal  quantum  number  3,  lie  already  around  113,500  cm“^.  The 
(  D)3s  D  state  is  followed  by  closely  spaced  triplets  with  the  valence  electron 


87 


in  levels  with  principal  quantua  number  4  and  higher.  Once  an  atom  has  been 

excited  into  these  levels,  it  has  a  much  higher  probability  of  being  further 

4 

excited  into  an  adjacent  S-triplet  than  of  going  into  the  high  triplet  states 
2  2 

of  D  or  P  parentage.  Transitions  from  these  doubly  excited  triplets  of  low 

4 

principal  quantum  number  into  the  relatively  underpopulated  S-triplets  of  high 
principal  quantum  number  will,  in  turn,  also  be  an  infrequent  process.  This 
will  lead  to  the  predicted  splitting  of  the  Boltzmann  plots. 

When  progressing  to  higher  energies,  the  Boltzmann  plot  for  the  doubly 
excited  triplet  states  will  again  show  a  drop  before  converging  to  the  ionization 
limit  for  ^  and  ^P,  These  excited  ions  can  be  formed  by  ionization  of  singlet 
and  of  triplet  atoms,  and  a  split  of  each  Boltzmann  plot  is  predicted  and  indi¬ 
cated  qualitatively  in  Fig,  15, 

The  population  distribution  over  oxygen  energy  levels  is  further  discussed 
in  Section  11. 


B.  Occupation  of  Energy  Levels  in  the  Ionization  Limit 

Since  the  concentration  of  molecule-ions  in  the  water  /apor  plasma  is  negli¬ 
gible,  and  second  ionization  of  oxygen  can  still  be  disregarded,  the  condition 
of  plasma  quasi -neutrality  requires  the  eiedtron  concentration  to  be  equal  to  the 
added  concentrations  of  free  protons  and  oxygen  ions. 


Thust  “  X  n 

e 

where  i  x 


CpJ  =• 
[o^J 

[0*]  +  [p1 


(1  -  x) 

0  <  X  <  1 
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Protons  do  not  have  discrete  energy  states,  but  in  order  to  obtain  the 
occupation  of  ground  state  oxygen  ions,  we  need  the  partition  function,  i.e, 
the  population  distribution  for  ions.  Next  to  the  statistical  weight  of  the 
ground  state  [g(l)  -4],  the  largest  coniribution  comes  from  the  low-lying  meta¬ 
stable  states, which  are  not  likely  to  be  in  a  Boltzmann  equilibrium  with  the 

ground  state  (see  Section  11 ).  However,  in  view  of  the  approximate  nature  of 

this  calculation,  the  possible  error  introduced  by  using  Qq+  “  QQ+(Eq. ,25000K) 

is  not  serious.  With  n  -  8  x  10^^  cm“^,  T  -  25,000  K,  the  occupation  of  the 

6 

last  bound  states  of  the  oxygen  and  hydrogen  atom  is  given  by  the  Saha  equation 
(Eq.  7-9)  asj 

log  n*(«>)  “  6.69  +  log  X 

log  n*(<»)  «  7.52  +  log  (1  -  x) 

Considering  the  bounds  for  the  ratio  x,  it  can  be  seen  that  even  for 

[p1  «  10^  )  ,  n*(»)  is  lower  by  an  order  of  magnitude  than  the  linear  extension 
0 

of  the  measured  population  distribution  for  the  quintet  states  of  oj^gen,  shown 
in  Pig.  15.  Conversely,  even  for  [p]  »  ,  "the  distribution  temperature  be¬ 

tween  the  sixth  energy  level  of  hydrogen  and  E  must  be  lower  than  the  last  meas- 

eo 

ured  distribution  temperature  of  2100  K, 

In  order  to  show  the  trenifor  the  Boltzmann  plot  at  high  quantum  numbers,  as 
predicted  by  theory,  we  have  arbitrarily  set  x  »  0,5  t  ba,sed  on  the  higher  diffu¬ 
sion  rate  for  protons  and  on  the  assumption  of  approximately  equal  rates  of  for¬ 
mation  of  oxygen  and  hydrogen  ions. 
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9.  MECHANISMS  LEADING  TO  DEPARTURES  FROM  LTE  IN  ARC  PLASMAS 


A.  Diffusion  of  Ground  State  Atoms  and  Aabipolar  Diffusion  of  Charged  Particles 
The  axial  zone  ("core”)  of  the  water  vapor  arc  at  high  currents  and  reduced 
pressure  is  characterized  by  a  large  overpopulation  of  ground-state  neutrals 
with  respect  to  an  occupation  in  equilibrium  with  their  ions  and  free  electrons 
at  the  electron  temperature.  The  excitation  and  ionization  rates  greatly  exceed 
the  rates  for  the  inverse  processes,  as  follows  from  the  principle  of  detailed 
balancing.  Ions  are  therefore  constantly  formed,  and  ground  state  neutrals  con¬ 
stantly  lost,  either  by  direct  transitions  into  the  continuum  or  by  the  step-wise 
excitation  process.  In  order  to  maintain  the  plasma  in  a  steady  state,  there 
must  be  a  corresponding  Influx  of  ground-state  neutrals  and  an  outflow  of  ions 
and  free  electrons. 


The  distance  an  atom  in  energy  level  i  can  diffuse,  on  the  average,  before 
undergoing  an  exciting  collision,  is  given  approximately  byj 


d(exc, i)  «  X 


t(exc , i) 


«  \fx  V  t(exe,l) 


(9-1) 


where  X  and  t^  are  the  mean  free  path  and  the  mean  free  time  for  elastic  col¬ 
lisions,  respectively,  v  is  the  root  mean  square  velocity  of  the  atom,  and 
t(ex,i)  the  free  time  for  exciting  collisions. 

The  pertinent  elastic  cross-sections  are  of  the  order  of  1G”^^  cm^.  The 
kinetic  temperature  of  heavy  particles  (gas  temperature)  in  the  core  of  the  arc 
at  reduced  pressure  has  been  determined  as  T^  »  5600  K.  The  diffusion  length 
for  hydrogen  atoms  is  then  estimated  byt 

d(exc,i)  «  0,15  1/10^  t(exc,  i)  (cm) 


90 


¥e  have  seen  previously  that  in  this  plasma  excitation  of  hydrogen  atoms 
occurs  mednly  by  a  step-wise  process.  The  average  lifetime  with  respect  to 
collisional  excitation  is  then  given  byi 

t(exc,i)  »  n(i)  /  (s)  (9-2) 

Based  on  the  measured  number  densities  and  the  calculated  collisional  excitation 
rates  (see  Table  XII ),  we  obtain  for  the  hydrogen  ground  state i  t(exG^l)  =»  1,6  x  10  ; 

d(exc^l)=^  0.6  cm;  and  for  the  first  excited  statei  t(exc^  2)  »  5,7  X  10  ^  3} 

d(exc^  2)  =  0. 01  cm. 

Since  the  diffusion  length  for  the  hydrogen  ground  state  is  larger  than  the 
radius  of  the  arc  channel  by  a  factor  of  the  order  of  3f  It  can  readily  be  seen 
that  diffusion  of  ground  state  atoms  from  the  arc  perij^ery  will  maintain  the 
high  overpopulation  in  the  arc  center.  The  concentration  of  atoms  in  the  first 
excited  state  will  not  change  significantly  over  a  radius  comparable  to  their 
diffusion  length,  and  the  level  occupation  will  then  not  be  sensibly  affected  hy 
diffusion,  Collisional  lifetimes  scale,  at  least,  as  the  occupation  of  each  level, 
and  diffusion  will  still  be  less  important  for  the  higher  bound  levels. 

In  evaluating  the  free  path  for  ions,  ion-ion  encounters  aLLso  have  to  be 
taken  into  account.  Although  in  this  context  one  only  has  to  consider  strong 
collisions  (i,e.  those  deflecting  the  scattered  particle  by  more  than  90^),  the 
pertinent  cross-sections  are  larger  than  ion-atom  cross-sections  by  almost  an 
order  of  magnitude.  However,  since  this  plasma  is  only  2%  ionized,  ion-ion 
collisions  can  still  be  neglected,  and  the  mean  free  path  for  protons  is  then 
practicably  the  same  as  for  hydrogen  atoms.  This  £w:tually  might  involve  an  under¬ 
estimate  because  of  the  correlation  between  the  diffusion  of  ions  and  of  electrons 
(ambipolar  diffusion),  that  leads  to  a  diffusion  rate  for  the  ions  twice  ais  large 
as  for  ordinary  diffusion. 
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We  have  used  E<|,  8-8  that  applies  to  the  rate  coefficients  for  collisional 

ionization,  in  order  to  derive  the  relation  of  Eq,  9-4,  and  it  is  then  only 

valid  for  recombination  into  excited  levels.  Recombination  into  the  ground 

state  occurs  predominantly  by  the  radiative  process  (see  Section  8a),  yielding 
—2 

t(rec,l)  «  10  s.  The  lifetime  of  a  free  proton  with  respect  to  recombination 
into  the  second  energy  level  is  of  the  same  order  of  magnitude.  Because  of  the 
small  probability  for  volume  recombination  into  lower  energy  levels,  a  steady- 
state  can  be  maintained  with  a  large  overpopulation  of  these  levels  compeured 
to  their  equilibrium  values. 

Even  for  high  quantum  numbers,  say  i  «  20,  we  obtain  with  t(rec,20)  «  lO"*^  s, 
a  diffusion  length  of  the  order  of  0,45  A  free  proton  formed  in  the  arc 

core  can  thus  be  lost  by  ambipolar  diffusion  to  the  arc  walls  before  recombining 
into  these  high  energy  levels, 

Eqs.  9-1  and  9-4  yield  the  following  approximate  scaling  lawi 

T  3 A  T 

d(ree,i)  o<  —^-7 - —  (9-5) 

^2-plf2 
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Eq..  9-5  has  been  derived  for  hydrogen  plasmas  in  the  limit  of  negligible 
ionization,  and  has  then  only  a  restricted  validity  for  other  plasmas,  but  is 
useful  for  an  order  of  magnitude  estimate.  Thus  one  can  expect  departures  fr<» 
the  Saha  equilibrium  in  a  plasma  due  to  ambipolar  diffusion  of  ions  and  electrons. 


ifi 


pI 


<  A 


(9-6) 


where  the  ratio  of  electron  temperature  to  electron  density  is  expressed  as  the 
square  of  the  Debye  radius,  is  the  total  concentration  of  particles,  and  L  is 
a  typical  linear  dimension  of  the  plasma  over  which  significant  changes  in  the 
plasma  state  occur,  i,e,  L  is  a  measure  of  the  spatial  inhomogeneity.  A  in  Sq.  9-6 
is  not  a  constant,  but  depends  on  the  plasma  composition  via  the  elastic  cross- 
sections.  Because  of  the  larger  ion-ion  cross-section,  A  decreases  with  in¬ 
creasing  degree  of  ionization. 

It  is  of  interest,  to  apply  these  concepts  to  the  low-current  water  vapor 
are  in  an  argon  matrix  at  atmosj^eric  pressure.  For  the  plasma  in  the  arc  core 
we  obtain  a  diffusion  length  for  proton  recombination  which  is  smaller  by  more 
than  two  orders  of  magnitude  than  the  equivalent  length  in  the  case  of  the  low- 
pressure  arc.  This  explains  why  only  mild  departures  from  LTC  were  observed. 

At  the  perij^ery  of  the  atmospheric  arc,  the  ratio  of  electron  temperature 
to  electron  density  is  practically  the  same  as  in  the  core  of  the  low  pressure 

-I/2 

arc,  and  the  respective  diffusion  lengths  scale  then  only  as  P  '  ,  or  1*6, 

However,  this  plasma  region  corresponds  to  an  arc  zone  of  very  steep  gradients 
of  electron  density  and  temperature  (i,e,  L  in  relation  9-6  is  very  small),  which 
will  lead  to  a  high  probability  for  the  ions  to  diffuse  outward  and  recombine  with 
electrons  via  wall -recombination.  Thus  we  observed  a  “dropping-off**  of  the 


Boltzmann  plots  for  the  excited  states  of  argon  and  hydrogen  entirely  equivalent 
to  our  measurements  of  the  center  region  of  the  arc  at  reduced  pressure, 

B,  Electronic  Excitation  hy  Heavy-Particle  Collision  and  Non-Haxwellian 
Velocity  Dlstrihution  for  Electrons 

The  author  is  not  aware  of  prior  similar  work  with  are  plasmas  at  high 
currents  and  reduced  pressure.  Our  Measurements  at  atmospheric  pressure,  how¬ 
ever,  can  be  compared  with  experimental  work  by  Kolesnikow,  who  has  reported 
departures  from  LTE  in  low-current  arcs  burning  in  a  mixture  of  argon  and  hydro¬ 
gen  at  atmosi*ieric  pressure  (22) •  He  also  observed  a  decrease  of  the  distribution 
temperature  over  the  higher  excited  states,  but  attributes  this  to  an  entirely 
different  effect, 

Kolesnikow  argues  that  at  low  electron  densities  encounters  with  atoms  - 
notwithstanding  the  much  smaller  cross-sections  -  can  be  more  effective  than 
electron  collisions  in  inducing  transitions  between  higher  levels.  The  gas  temper¬ 
ature  in  these  plasmas  is  substantially  lower  than  thedectron  temperature,  and 
because  of  the  threshold  condition,  the  effectiveness  of  excitations  induced  by 
collisions  with  heavy  particles  will  increase  with  a  decrease  of  the  energy  jump 
involved  in  the  transition.  Thus  excitation  of  atoms  in  the  lower  energy  levels 
still  occur  by  electron  collision  because  the  energy  gaps  are  much  larger 
than  the  mean  kinetic  energy  of  heavy  particles,  but  excitation  by  atom-atom 
collision  should  become  progressively  more  important  for  the  higher  energy  levels. 
One  should  then  observe  an  assymptotic  approach  of  the  distribution  temperature 
to  the  kinetic  temperature  of  heavy  particles. 

Kolesnikow’s  interpretation  of  observed  departures  from  LTE  in  arc  plasmas 
at  low  electron  densities  appeaue  to  be  generally  accepted  (^,^).  It  is  however, 
irreconcilable  with  the  results  of  more  refined  measurements. 
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In  the  peripheral  zone  of  our  arc  at  atmospheric  pressure,  the  measured 

14  -3 

electron  density  was  2,5  10  cm  ,  which  is  lower  by  a  factor  of  about  three 

thain  in  the  plasma  of  lowest  electron  density  investigated  by  Kolesnikow,  But 
even  at  the  high  ratio  of  atoms  to  electrons  of  the  order  of  5000,  the  observed 
departures  from  LTE  are  surely  due  to  diffusion  and  not  to  excitation  by  atom 


collisions, since  the  measured  distribution  temperatures  over  the  higher  excited 
states  of  hydrogen  and  of  argon  are  substantially  lower  than  the  gas  temperature. 
Atom-atom  collisions  can  be  more  effective  thaua  electron-atom  encounters, 
if  the  ratio  of  atoms  to  free  electrons  in  the  plasma  fulfills  the  following  con¬ 


dition  (41  )t 


(9-7) 


The  function  f  in  Eq,  9-7  is  determined  by  the  specific  energy  dependence 
of  the  pertinent  cross-sections,  and  is  difficult  to  estimate  for  the  range  of 
interest,  i,e.  for  a  mean  kinetic  energy  of  the  colliding  particles  much  laurger 
than  the  threshold.  Since  >  T^,  we  assume  that  in  this  range  f  is  not  sig¬ 
nificantly  smaller  than  unity.  Electronic  excitation  by  collision  with  atoms 

will  then  only  become  importauit  if  the  concentration  of  atoms  exceeds  the  electron 

4 

density  by  a  factor  of  the  order  of  5  10  in  a  hydrogen  plasma  -  the  most 

6 

favorable  case  -  and  of  the  order  of  10  in  an  argon  plasma. 

At  atraosiiieric  pressure  auid  belon,  the  electron  densities  required  to  main¬ 
tain  an  arc  discharge,  are  always  higher  than  the  critical  density  determined  by 
this  estimate,  and  it  can  be  concluded  that  for  arc  plasmas  eieetronic  excitation 
by  heavy  particle  encounters  is  probably  always  negligible. 
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Our  treatment  of  the  interaction  processes  in  the  plasma  is  based  on  the 
assumption  of  a  Maxwell  velocity  distribution  for  the  electrons.  Equilibration 
of  the  velocity  distribution  occurs  mainly  by  electron-electron  interaction  because 
this  process  is  very  efficient,  and  the  pertinent  cross-sections  are  large.  The 
electron-electron  relaxation  time  can  be  estimated,  according  to  Spitzer  (42) . 
with: 

«  0.03  (s)  (9-8) 

which  yields  t^^  «  1,5  x  s  for  the  plasma  in  the  arc  core  at  reduced  pressure 

and  t  «  0,7  x  s  for  the  peripheral  zone  of  the  arc  at  atmospheric  pressure 

and  loM  current.  These  times  are  very  short,  much  shorter  them  the  recombination 
times  (Sq,  9“^).  Since  electrons  amd  ions  diffuse  together,  and  the  ambi polar 
diffusion  rate  is  determined  by  the  slowiaoving  ions,  the  diffusion  length  for 
an  electron  in  time  t  will  be  small  compared  to  typical  arc  dimensions.  The 
velocity  distribution  of  electrons  will,  therefore,  be  indeed  Maxwellian,  and  a 
well  defined  electron  temperature  can  be  assigned  to  small  radial  zones  of  the 
arc, 

Schluter  has  attributed  observed  anomalies  in  spectras  of  high-frequency 
discharges  to  a  Druyvestyn  velocity  distribution  for  the  electrons,  caused  by 
diffusion  of  the  energetic  electrons  from  the  tail  of  the  Maxwell  distribution  (42)- 
The  theoretical  treatment  of  this  question  is  beset  by  considerable  difficul¬ 
ties  (26),  No  evidence  for  this  effect  has  been  found  in  the  present  work. 
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in  the  Flasna 


So  far  we  hare  only  discussed  departures  from  equilibritm  with  respect  to 
species  concentrations  and  population  distributions  over  quantized  energy  levels. 
We  now  consider  the  lack  of  equipartion  between  the  kinetic  energy  of  electrons 
and  of  heavy  particles  that  leads  to  the  observed  discrepancy  between  the  elec¬ 
tron  tenperature  T  and  the  gas  temperature  T  , 

If  the  energy  loss  by  inelastic  collisions  can  be  disregarded,  in  a  steady 
state  the  energy  gained  by  an  electron  from  the  applied  electric  field  in  t^  , 
the  mean  free  time  between  elastic  collisions  with  heavy  particles,  must  equal 
the  meam  energy  transferred  by  such  a  collision  (^). 

As.  -  tj  -  S.SxlO^E^t^  (ev)  (9-9) 

2n.  ^  * 

2i& 

A.E^(el)  -  I  k  (T^  -  Tg)  —  “  l.iK)  X  10“^  (T^  -  T^)/A  (eV)  (9-IO) 

ithere  temperatures  are  reckoned  in  K,  and  A  is  the  atcMiic  weight  (mass  ntnaher) 
of  the  collision  partner.  In  the  numerical  form  of  Eq.  9-9,  the  field  strength 
is  expressed  in  v/cm.  E  is  experimentally  determined  "by  meaisuring  the  Toltage 
drop  along  the  arc  channel  (see  Section  3G), 

By  setting  Ae^  »  Ae^  »  and  with  obtain  the  electron 

free  path  with  respect  to  elastic  collisions  with  heavy  particles  i 


1  «  8,5  X  10‘ 


T,  -  T  I  1/2 


(cm)  (9-11) 


-4 

Eq,  9-11  was  used  to  calculate  X  ■»  8,4  x  10  cm  from  the  measured  data 

e 

corresponding  to  the  plasma  in  the  core  of  the  arc  at  atmospheric  pressurei 


97 


X  »>  9300  K,  T  =  5^0  K,  E  ■  15  V/cm.  The  chemical  composition  of  the  plasma  is 
approximately  At  with  only  1%  H  and  \%  0,  but  one  has  to  use  X  »•  28  in 
Eq.  9-11,  by  properly  averaging. 

Since  this  plasma  is  only  0.2^  ionized,  we  can  disregard  collisions  with 

-1  -I*!  2 

ions  and  set  \  =•  cj  N.  ,  yielding  cr  =1, Ox  10  cm.  This  agrees  to 

6cL  w  Get 

better  than  20^  with  measured  eleistic  cross-sections  of  argon  for  collisions 
with  slow  electrons. 

If  the  long-range  interaction  between  charged  particles  is  taken  into  account, 
the  cross-sections  for  electron-ion  encounters  aire  much  larger  than  for  scattering 
by  neutrals.  For  approximate  calculations  one  can  use  an  expression  valid  for 
all  singly  charged  ions  (42): 

{Jg!  »  3  X  10"^  t'^  (cm^)  (9-12) 

Thus  at  T  «  26,000  K,  the  estimated  electron  temperature  in  the  core  of 
6 

the  water  vapor  arc  at  reduced  pressure,  «  5  x  10  cm  ,  which  is  50  times 

larger  than  the  cross  section  of  hydrogen  with  respect  to  elastic  collisions  with 

electrons.  At  2%  ionization,  the  freq^ueney  of  electron-ion  and  electron.-atom 

encounters  are  about  equal  in  this  plasma.  With  the  calculated  \  -  1,3  x  lO"^  cm, 

e 

we  obtain  t^  *•  1,2  x  10  s.  The  measured  field  strength  was  26  V/cm*  Thus,  as 

given  by  Eq,  9-9,  an  electron  gains  8*5  10  ^  eV  from  the  applied  electric  field 

in  the  free  time  between  elastic  collisions.  However,  with  (T  -  T  )  *  20,000  K 

®  S 

and  A  «  1.5  ,  the  mean  energy  transferred  to  a  heavy  particle  by  an  elastic 
collision  is  only  of  the  order  of  1.9  x  lo"*^  eV  (Eq.  9~10) .  This  is  less  than 
25^  of  the  energy  gained  from  the  field.  The  balance  must  then  be  dissipated 
by  inelastic  collisions. 
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In  diffusion-controlled,  non-equilibrium  plasmas,  the  electrons  lose  energy 
by  inelastic  collisions,  because  excitation  processes  are  not  balanced  by  their 
inverse,  i.e,  collisions  of  the  first  kind  are  much  more  frequent  than  collisions 
of  the  second  kind.  The  kinetic  energy  of  electrons  converted  per  unit  volume 
and  per  second  into  potential  energy,  by  exciting  the  atoms  in  energy  level  i 
to  energy  level  j,  is  given  by  the  product  of  the  energy  jump  As(i,j)  involved, 
the  eollisional  excitation  rate,  and  the  factor  B(l,j)  defined  in  Section  8B  as 
the  fraction  of  exciting  collisions  that  is  not  balanced  by  de-excitlng  collisions. 
As  has  been  sho'^n  previously,  in  the  case  of  hydrogen  no  serious  error  is  intro¬ 
duced  by  disregarding  all  other  eollisional  excitation  processes,  including  direct 
ionisation,  and  considering  only  step-wise  excitation.  This  can  be  virittenj 

Rgoi(^-.i+l)  -  const,  (all  i)  *  ,  (9-13) 

which  will  be  reccgninzed  as  Eq,  8-14  written  in  another  form.  Sincei 

oo 

E  iAE(i,i+l)  »  , 

the  average  energy  loss  per  electron  due  to  inelastic  collisions  in  time  t^  , 
can  thus  be  estimated  withj 

AE_(exc)  =  t^  /  n^  (9-14) 

In  the  arc  core  at  reduced  pressure,  B(l,2)  «  1,  and  we  can  set 
ef 

^col  ~  ^col^^’^^’  whose  value  is  shown  in  Table  XII.  This  yields: 

AE^(exc,H)  ^  4  X  10  ^  eV.  Assuming  that  Zi^^(exc,0)  is  of  the  satme  order, 
we  have: 

E  Ae^  =  z2iE^(el)  +  AE^(exc,H)  +  AE^(exc,0)  *  7.9  x  10*^  eV  =  Ae^  , 
which  comes  quite  close  to  the  Ae_^  calculated  before. 
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Lack  of  equipartition  between  the  kinetic  energy  of  electrons  and  heavy 
particles  is  generally  associated  with  plasnas  at  low  pressure.  It  has  been 
shown  here,  however,  that  even  at  atmos^eric  -pressure  a  substantial  divergence 
between  the  electron  and  the  gas  temperature  of  the  plasma  can  arise,  if  the 
degree  of  ionization  is  very  low.  Because  of  the  small  neutral  cross-sections, 
the  mean  free  path  is  then  so  long  that  the  electrons  can  gain  a  relatively  large 
amount  of  energv  from  the  applied  electric  field.  Ion  cross  sections  are  typical¬ 
ly  50  times  larger,  and  if  the  plasma  is  only  a  few  percent  ionized  ,  the  free 
path  will  be  substantially  decreased.  Elastic  collisions  with  neutrals  can  then 
be  disregarded  i”.  a  first  approximation,  and  by  combining  Eqs,  9-11  and  9-12,  one 
obtains  a  scaling  law  for  the  lack  of  equipartition: 


•r 

'e 


I 


(9-15) 


The  dependence  on  Tg/n^  is  thsrf fore  still  stronger  than  in  the  case  of  de- 
■partures  from  an  equilibrium  distribution  over  energy  levels  due  to  diffusion 
(Eq.  9-6). 

One  also  has  to  examine  if  equipartition  of  kinetic  energy  exists  between 
the  ions  and  neutrails  in  the  plasma.  In  plasmas  that  are  in  LIE  and  are  sub¬ 
stantially  ionized,  the  kinetic  temperature  of  ions  and  of  neutrals  will  tend 
to  equilibrate  via  the  ionization-recombination  process,  which  is  then  more 
®i*f®ctive  than  elastic  collisions  between  ions  and  neutrals  (^).  This  mechanism 
is  not  o-perative  in  diffusion-dominated  plasmas,  however,  because  of  the  low 
probability  for  volume  recombination  (Section  9A).  In  the  core  of  the  low- 
pressure  arc  irvestigrated  in  this  work,  at  2%  ionization,  the  frequency  of  elastic 
electron-neutral  and  electron-ion  collisions  are  about  equal.  At  slightly  higher 
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ionization,  most  of  the  kinetic  energy  transfer  to  heavy  particles  will  already 
occur  by  electron-ion  interaction.  Since  distant  encounters  between  charged 
paa^lcles  again  have  to  be  taken  into  account,  the  cross-sections  for  elastic 
ion-ion  collisions  are  much  larger  than  for  ion-neutral  collision.  Equi parti¬ 
tion  of  kinetic  energy  between  ions  and  atoms  depends  then  to  a  large  extent 
on  the  effectiveness  of  the  charge-transfer  mechanisms,  which  is  difficult  to 
estimate.  In  prior  work  with  high-current  argon  arcs  at  reduced  pressure,  we 
have  found  a  notable  divergence  between  the  kinetic  temperatures  of  ions  and  of  atoms 
(45);  although  the  ion  temperatures  reported  in  Ref.  45  are  too  high  because  the 

calculated  Stark  broadening  parameter  was  used  to  correct  the  half  width  of  the 

o 

Doppler-broadened  Ar  II  4348  A  line  (see  Section  3D). 


10,  MOLECULAR  SPECIES  IN  DIFFUSION-DOMINATED  PLASMAS 


A,  Excitation  Processes 

We  Hill  noH  apply  the  concepts  developed  in  the  preceding  Sections  to  a 
discussion  of  the  Interaction  processes  in  non-equilibriua  plasmas  leading  to 
the  observed  population  distributions  over  molecular  energy  levels. 

Also  for  molecules  such  as  OH,  electronic  excitation  and  ionization  and 
the  inverse  processes  are  dominated  by  colllsionaJ.  interaction  nith  free  elec¬ 
trons.  A  close  correspondence  will  therefore  exist  between  atomic  and  molecular 
species  with  respect  to  the  population  distribution  over  excited  electronic  states 
and  the  ionization-recombination  equilibrium. 

For  a  diffusion-dominated  water  vapor  arc  one  would  predict  a  population 
distribution  between  the  ground  and  first  excited  electronic  state,  described 
by  a  distribution  temperature  much  higher  them  the  gas  temperature  and  comparable 
to  the  distribution  temperature  pertaining  to  the  occupation  of  the  first  and 
second  energy  level  of  hydrogen.  The  distribution  temperature  will  fall  sub¬ 
stantially  when  progressing  to  higher  energies,  leading  to  a  decreased  occupation 
of  the  higher  electronic  states.  The  population  of  the  Rydberg  states  and  the 
concentration  of  OH  will  then  be  very  small. 

The  first  of  these  predictions  was  indeed  verified  by  our  experimental 
results,  and  also  no  trace  of  the  OHCb^S  -  A^S)  band  at  4336  A  was  found.  The 
search  for  the  OH  (A^n  -  X%)  bands  was  inconclusive,  because  these  bands  fall 

in  a  very  crowded  spectral  region  due  to  the  strong  excitation  of  the  Schumann- 
Runge  bands. 
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We  have  shown  that  in  the  plasma  the  distribution  over  the  rotational- 
vibrational  levels  of  0H(A  S),  that  axe  not  affected  by  predissooiation,  is 
close  to  a  Boltzmann  distribution  at  the  gas  temperature.  It  must  be  concluded, 
therefore,  that  the  rate  coefficients  for  excitation  of  nuclear  motion  by  electron 
impact  are  not  significantly  larger  than  the  rate  coefficients  for  the  equivalent 
process  Induced  by  heavy  particle  encounters,  so  that  at  low  electron  densities 
the  latter  predominates. 

Because  of  the  small  optical  transition  probabilities  for  the  transition 
A^S-»X^n,  radiative  decay  to  the  ground  state  is  negligible  compared  to  col- 
lisionally  induced  transitions  between  the  vibrational  and  rotational  levels  of 

p 

the  AD  state.  To  a  ^ood  approximation,  these  excitation  rates  are  proportional 
to  the  transition  moment  aoid  therefore  subject  to  the  same  selection  rules. as 
optical  dipole  transitions.  A  rotational  level  J  can  then  combine  only  with 
the  rotational  levels  ±  1#  and  rotational  excitation  and  de-excitation  will 
occur  by  a  step-wise  process.  The  threshold  is  smeill  compared  to  the  mean  kinetic 
ener^  of  heavy  particles,  and  the  exponential  factor  is  close  to  unity.  The 
magnitude  of  the  excitation  rate  for  the  transition  J  J  +  1  is  therefore  de¬ 
termined  by  the  transition  moment,  proportional  to  (J  +l)/(2J  +  l).  With  the  ex¬ 
ception  of  the  lowest  rotational  levels,  the  excitation  rate  is  thus  practically 
Independent  of  the  rotational  quantum  number.  This  can  be  contrasted  with  col- 
lisional  excitation  of  hydrogen  atoms,  where  transitions  into  other  than  neigh¬ 
boring  levels  are  less  probable,  but  not  strictly  forbidden.  The  rates  for  the 

4 

transition  i-^i  +  1  increase  as  i  for  the  higher  energy  levels  of  hydrogen, 

2 

The  occupation  of  the  A  x;  state  of  OH  is  mainly  determined  by  the  balance 

between  transitions  from  the  ground  state  and  into  the  next  higher  electronic 

2  2 

state,  induced  by  electron  collision.  The  excitation  process  X  I[-»A  2  occurs, 
of  course,  according  to  the  Franck-Gondon  principle;  moreover  the  electronic 
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excitation  rate  is  also  independent  of  the  nuclear  motion  of  the  excited  molecule. 


The  distribution  over  the  vibrational -rotational  levels  of  the  ground  state  is, 
therefore,  "frozen"  into  the  population  distribution  for  the  electronically  ex¬ 
cited  molecule. 

Since  dissociation  occurs  by  heavy  particle  collision,  the  equilibrium  occu¬ 
pation  of  energy  levels  with  respect  to  the  dissociated  state  is  determined  by  the 

gas  temperature.  Because  of  the  high  distribution  temperature  characterizing  the 

2  2 

relative  occupation  of  the  A  S  and  the  X  U  state,  the  overpopulation  of  the  elec¬ 
tronically  excited  state  will  always  be  larger  than  for  the  ground  state?  in  the 
low  pressure  arc,  in  fact,  by  several  orders  of  magnitude.  Experimentally,  it 
has  been  established,  however,  that  the  population  distribution  over  the  lower 
rotational  and  vibrational  levels  of  the  state  is  close  to  a  Boltzmann  dis¬ 
tribution  at  the  gas  temperature.  The  probability  for  a  eollisionally  induced 
transition  into  the  dissociated  state  must  therefore  be  very  small  for  these  states. 
This  is  discussed  more  fully  in  the  next  Section, 

Ihe  concentrations  of  0^  2uid  of  will  be  in  a  quasi -equilibrium  with  OH  at 
the  gas  temperature,  because  collisions  with  heavy  particles  are  effective  in 
Inducing  reactions  of  the  types 

OH  +  0  H  +  O2 

OH  +  H  ^  °  ^2 

whose  energy  threshold  is  comparable  to  the  mean  kinetic  energy  of  heavy  parti¬ 
cles,  Because  of  the  higher  gas  temjseratures  in  arcs,  the  pertinent  rate  coeffi¬ 
cients  will  be  higher  than  in  flames,  but  in  principle  the  same  considerations 
apply  The  experimental  evidence  for  a  substantial  concentration  of  0^  in 

the  arc  plasma  is  given  by  the  prominence  of  the  Schumann-Runge  bands  belonging 
to  transitions  between  highly  excited  vibrational  levels. 
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B,  The  Dissociation  Equilibrium  of  OH  and  Predlssoelatlon 

TTie  Saha  equation  is  only  a  specialized  form  of  the  general  mass  action 
equation,  and  the  conclusions  reached  in  discussing  ionization-irecombination 
equilibria  apply  also  to  the  dissociation-recombination  equilibrium.  Dissocia¬ 
tion  of  ground  state  OH  can  be  compared  with  the  previously  discussed  ionization 
of  ground  state  hydrogen,  with  heavy  particles  taking  the  part  of  free  electrons. 

The  transition  moment  for  direct  transitions  into  the  dissociated  state  is  small, 
and  the  threshold  energy  is  large  compared  to  thermal  energies.  It  follows  that 
the  free  time  of  a  OH(X^n )  molecule  with  respect  to  a  dissociating  collision  is 
much  lomger  than  the  free  time  between  elastic  collisions.  It  can  thus  diffuse 
over  a  distance  comparable  to  the  arc  diameter  before  being  effectively  dissocia- 
ted* 

Because  of  the  small  probability  for  volume  recombination,  ground  state  OH 
will  mostly  be  formed  by  the  effective  mechanism  of  wall -recombination.  The  con¬ 
centration  of  OH  in  the  arc  center  is  then  not  determined  by  an  equilibrium  rela¬ 
tion,  but  by  the  balance  between  the  rate  of  diffusion  of  OH  from  the  arc  peri- 
lAiery  and  the  excitation  rate  into  the  A  S  state,  induced  by  electron  impact. 

Thus  we  found  the  concentration  of  ground  state  OH  to  be  almost  constant  over  the 
entire  cross-sectiop  of  the  atrc  at  reduced  pressure,  and  higher  than  the  equili¬ 
brium  concentration  at  the  gas  temperature  by  almost  two  orders  of  magnitude.  In 
the  arc  at  atmospheric  pressure  an  overpopulation  of  OH(X  IE )  has  also  been  es¬ 
tablished  (Sections  5^  and  5®). 

Gaydon  and  Wolfhard  (^)  have  discussed  the  pertiirbation  of  the  A  state 

by  a  lower  electronic  state,  leading  to  a  predissociation  belonging  to  Herzberg's 

case  Ic  (48)*  The  perturbing  state  is  most  likely  a  S  state,  and  its  dissociation 

2 

products  are  normal  (i.e,  ground  state)  atoms,  whereas  the  A  2  state  dissociates 


105 


into  ground  state  hydrogen  and  excited  (^)  oxygen  atoms.  The  effective  potential 

energy  curves  (i.e.  including  rotation)  cross  at  about  E*  «•  42,000  cm*"^  for 

2 

^  “0  of  the  A  2  state  and  at  progressively  lower  energies  for  the  higher  vibra” 
tinal  levels.  An  unusual  characteristic  of  the  predissociating  rotational  levels 
is  the  splitting  of  the  Boltzmann  plots  for  the  two  spin  components  (see  Fig,  25), 
level  with  J  *  K  +  l/2  is  weakened  because  the  probability  for  a  transi¬ 
tion  into  the  perturbing  electronic  state  is  higher  than  for  the  f^  level  of  same 
K  (^). 

The  A  electronic  state  of  OH  is  a  state,  and  the  selection  rules  for 
Hund  s  case  b  apply,  i.e,  in  addition  to  the  strict  selection  rule  for  dipole 
transition i  J— »  J  ±  1,  we  also  have  the  selection  rule  K  — »  K  +  1.  Combinations 
between  the  f^  and  f^  levels  are  not  strictly  forbidden  and  correspond  to  the 
weak  satellite  branches,  whose  intensity  with  respect  to  the  main  branches  de¬ 
creases  with  increasing  J,  Similarly  for  collisional  processes,  if  the  f^  levels 
are  more  strongly  depopulated  by  predissociation  than  the  f^  levels,  the  splitting 
of  the  Boltzmann  plots  should  increase  with  J  because  of  the  decreased  probability 
for  f2->f^  transitions.  This  is,  in  fact,  what  is  observed  and  can  be  taken  as 
further  evidence  that  the  selection  rules  for  optical  dipole  transition  also  apply 
to  transitions  induced  by  collisions. 

As  is  seen  most  clearly  in  the  case  of  the  (2,l)  and  (2,2)  bands  (see  Fig.  25), 
after  the  onset  of  predissociation  a  short  transition  region  is  followed  by  a 
range  of  rotational  energy  levels  for  which  the  Boltzmann  plot  is  again  linear, 
although  its  slope  does  not  correspond  to  any  meaningful  temperature.  This  can  be 
explained  as  follows.  In  this  range  a(j)/a(j-l)  «  1,*^  and  de-excitations  can 


*)  a(j)  is  equivalent  to  a(i)  defined  in  Eq.  7-7 
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be  disregarded  in  a  first  approximation.  The  occupation  of  a  predissoeiating 
rotational  level  J  is  then  mainly  determined  by  the  balance  between  transitions 
from  J  -  1  to  J,  and  from  J  to  J  +  1  and  dissociating  transitions  via  the  per¬ 
turbing  state: 

n(J-l)  X(J-1,J)  «  n(j)  [  X(J,J+1)  +  D(J)  1  (10-1) 

Since  AE(J,J+1)  «  kT  ,  the  exponential  factor  is  close  to  unity,  auid  the 
specific  dependence  of  the  excitation  cross-section  on  the  energy  of  the  colliding 
particle  is  unimportant.  The  excitation  rates  can  then  be  approximated  by  an 
expression  of  the  form: 

X(J,J+1)  »  ^  ^  C(T  )  (10-2) 

2J  +  1  S 


and  Eq,  10-1  can  be  written; 

n*(J-l)  „  J  +  1  ^  2J  +  1  I)(J) 
n*(j)  J  J  C(Tg) 

J  »  1  }  D(J)  «  C(T^):  «  exp  [  2D(J)/G(T^)  ]  (10-4) 

®  n*(j)  ® 


The  Boltssfliann  plot  will  thus  be  approximately  linear  over  the  range,  where; 

D(J)  -  D(T  )  «  G(T  )  (10-5) 

O  O 

As  shown  by  Kronig,  the  transition  moment  for  predissociations  for  which  A  A  »  0 
-  which  applies  to  this  case  -  is  independent  of  J  in  a  first  approximation  (49) e 
The  predissociation  also  has  to  be  considered  as  weak  [[d(J)  «  G(T  )J,  because 
it  violates  the  K — >K  rule  (4?,  48)^ 

The  observed  sharp  drop  in  the  occupation  of  the  highest  rotational  levels 

occurs  because  of  the  now  significant  probability  for  direct  dissociation  of  the 
2  -1 

A  E  state  into  H(i)  +  (”predissociation  by  rotation”).  As  the  energy  term 
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of  the  vibrational -rotational  levels  approaches  the  dissociation  energy  for  the 
2 

A  E  state,  the  occupation  of  these  levels  will  tend  to  an  equilibrium  occupation 
with  respect  to  the  dissociated  state, 

11,  EFFECT  OF  THE  AUGER  PROCESS  ON  POPULATION  DISTRIBUTIONS 

The  population  distributions  over  the  energy  levels  of  oxygen  and  over  the 

2 

predissociating  rotational  levels  of  the  A  E  state  of  OH  are  both  affected  by  an 
Auger  process  of  small  probability  as  compared  to  the  excitation  rates,  and  a 
similarity  between  the  respective  Boltzmann  plots  can  be  expected,  A  splitting 
of  the  Boltzmann  plot,  which  is  readily  observable  in  the  case  of  OH,  has  also 
been  predicted  for  oxygen  and  is  caused  by  an  essentially  similar  mechanism  as 
discussed  in  Section  7G, 

The  Boltzmann  plot  for  oxygen  is  also  approximately  linear  over  a  consider¬ 
able  range  of  energy  levels  (see  Fig,  15),  and  the  same  explanation  appears 
plausible:  Predominance  of  step-wise  excitation  and  near-constancy  of  the  ex¬ 
citation  and  ionization  rates.  It  is  certainly  only  accidental  that  the  measured 
distribution  temperatures  for  the  quintet  states  of  oxygen  are  close  to  the  gas 
temperature,  because  excitation  by  heavy  particle  collision  has  been  shown  to  be 
insignificant  in  this  plasma  as  eompsured  to  the  same  process  induced  by  electron 
impact. 

The  dissociated  states  of  OH  -  oxygen  atoms  in  the  and  ^D  levels  -  are 
analagous  to  the  ionized  states  of  oxygen  -  ions  in  the  ^S,  and  levels  - 
in  the  sense  that  in  both  cases  optical  dipole  transitions  between  these  states 
are  strictly  forbidden  by  the  Laporte  rule.  Because  of  a  correspondingly  small 
collisional  interaction  among  these  states,  their  occupation  will  be  effectively 
uncoupled  in  non-equilibrium  plasmas. 
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12.  CONCLUSIONS 


It  has  been  shown  in  this  work  that  the  integrated  radiance  of  the  rotational 
lines  belonging  to  the  OH  (A^E  -  X^Il)  systefc  can  be  represented  by  an  expression 
of  the  form 

N,  .  C  n*(n',T*,K')  P°  P’(q’,q")  exp  [-p(q*,<l")  F(k)]  S/X^  (l2-l) 

ij  o  o 

The  integrated  line  absorption  coefficient  -  i.e,  the  quantity  obtained  by  an  absorp¬ 
tion  neasurenent  at  high  resolution  -  is  then  related  to  the  occupation  of  the 
lower  state  for  the  transition  by  a  siallar  expression,  with  the  "line  strength 
factor"  s/X  substituted  by  SA. 

Eq.  12-1  follows  from  the  linearity  of  the  "apparent  Boltzmann  plot",  and 
has  been  verified  here  for  the  main  branches  of  the  (0,0),  (l,0),  (l,l),  (2,l), 
and  (2,2)  bands  below  the  onset  of  predissociation.  The  molecular  constant  de¬ 
noted  by  P°  -  formally  equivalent  to  an  "average  electronic  transition  moment"  for 

the  system  -  is  proportional  to  the  radial  part  of  the  transition  moment  for  the 

P-1  lines  of  the  (0,0)  hand  (see  Eqs*  2-27  and  2-3l)  that  corresponds,  in  a  sense, 

to  the  fundamental  transition  for  this  system  (the  •'vibration  and  rotation-less 

molecule”).  The  numericatl  value  for  could  not  be  ascertained  because  of  the 

o 

difficulties  in  determining  the  concentration  of  OH  in  the  plausma,  noted  below. 

The  quantities  P^(q*,<l”)  -  Pq(v*,v”)  have  been  shown  to  be  constant  for  a  parti¬ 
cular  band,  i.e.  they  do  not  depend  on  the  specific  branch,  and  are,  therefore, 
enitrely  equivalent  to  the  concept  of  Franck-Condon  factors.  The  dependence  of  the 
radial  part  of  the  transition  moment  on  the  rotational  transition  is  contained  in 
the  exponential  factors  p(q*,q”),  that  also  have  been  found  to  be,  in  general, 
independent  of  the  specific  branch  within  a  band.  An  exception  has  been  noted  for 
the  (l,0)  band,  but  this  may  only  be  spurious  (see  Section  6b),  Because  of  the 
ambiguity  in  assigning  an  average  frequency  to  a  band,  the  representation  of 


109 


rotational  line  intensities  (or  absorption  coefficients)  in  terms  of  band 
transition  probabilities  or  band  oscillator  strengths  is  not  appropriate. 

Implicit  in  our  derivation  of  the  transition  aoaents  is  the  assuaption 
that  the  occupation  of  the  non-predissociating  rotational-vibrational  eneris?y 
levels  is,  in  effect,  ^iven  by  a  Boltsmann  distribution  at  the  gas  teaperature. 
However,  even  when  the  population  distribution  over  energy  levels  shows  an  expo¬ 
nential  dependence  on  the  energy  term,  this  is  not  sufficient  proof  for  aui  equi¬ 
librium  distribution,  as  has  been  discussed  in  Section  1GB,  Since  departures 
froa  LTE  in  diffusion-dominated  plasmas  lead  to  lower  distribution 
temperatures  when  progressing  to  higher  energies,  the  P^-f actors  (Franck-Condon 
factors)  derived  in  this  work  constitute  only  a  lower  bound,  and  the  exponential 
p-f actors  only  an  upper  bound  for  the  true  values.  This  is  especially  evident 

for  the  bands  originating  in  the  v*  »  2  level,  where  the  P*  factors  measured 

o 

here  are  probably  too  Ibn^  and  the  p-factors  obviously  too  high. 

-  o  - 

Even  at  the  lowest  electron  densities  at  which  a  steady  are  discharge  can 
be  sustained,  electronic  excitation  of*  OH  by  electron  impact  is  more  effective 
than  by  heavy  particle  collision.  The  distribution  temperature  characterizing 
the  relative  occupation  of  the  ground  state  of  OH  and  the  excited  electronic 
state  exceeds  then  the  gas  temperature,  that  is  the  relevant  temperature  for  the 
dissociation-recombination  equilibrium.  The  state  will  thus  be  overpopulated 
with  respect  to  the  dissociated  state,  and  the  population  distribution  over  its 
higher  rotational rvibrational  levels  will  be  determined  by  predissociation.  If 
the  electric  current  of  an  arc  at  atmospheric  pressure  is  increased,  the  degree 
of  ionization  of  the  plasma  rises  sharply,  and  equipartition  between  the  kinetic 
energy  of  electrons  and  of  heavy  particles  will  eventually  be  approached  by  the 
mechanism  discussed  in  Section  9G.  Based  on  the  present  investigation  it  can  be 
predicted,  however,  that  this  will  not  occur  at  temperatures  below  10,000  K, 
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Since  OH  is  almost  completely  dissociated  at  this  temperatiire ,  it  does  not 
appear  possible  to  obtain  arc  plasmas  with  an  equilibrium  concentration  of  OH, 
that  arc  suitable  for  spectroscopic  measurements. 

By  varying  the  pressure  and  the  current  of  a  water  vapor  arc,  the  occu- 

2 

pation  of  a  particular  rotational -vibrationaJ.  level  of  the  A  E  state  of  OH  can 
be  maximised,  as  has  been  done  in  this  investigation.  However,  this  occupation 
cannot  be  calculated  with  precision  because  it  depends  to  a  large  extent  on  the 
effectiveness  of  the  various  diffusion  processes  ocurring  simtiltaneously  in  the 
plasma.  Although  the  Fowler-Milne  method  is  then  clearly  inapplicable,  this  arc 
regime  still  offers  the  significant  practical  advantage  of  an  increased  stability 
against  small  are  fluctuations. 

Because  of  the  small  probability  for  direct  transitions  into  the  dissociated 
state  or  into  the  adjacent  rotational -vibrational  levels,  the  step-wise  excitation 
and  de-excitation  mechanism  in  conjunction  with  the  near-equality  of  the  transi¬ 
tion  rates,  effectively  uncouples  the  lower  rotational -vibrational  levels  from 
the  dissociated  state.  The  relative  population  distribution  over  these  levels 

will  therefore  approach  an  equilibrium  distribution  at  the  gas  temperature,  even 

2 

when  the  total  occupation  of  the  A  E  state  is  orders  of  magnitude  higher  them  its 
equilibrium  occupation. 

The  relaxation  times  with  respect  to  the  dissociation-recombination  process 
are  still  much  longer  for  OH  molecules  in  the  ground  state,  which  can  diffuse 
over  a  considerable  distance  before  equilibrium  is  established.  The  relative 
population  distribution  over  the  excited  states  is  essmtially  determined  by  local 
parameters ,  and  can  then  be  derived  in  principle ,  even  in  absence  of  equilibrium , 
if  the  appropiate  rate  coefficients  are  known.  The  concentration  of  OH  can  be 
experimentally  measured  in  high-resolution  work  if  the  molecular  constant 
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o  2  2 

P^(A2“Xn)is  known,  as  hM  been  shown  here.  It  does  not  appear  possible, 
however,  to  obtain  the  concentration  of  OH  in  ares  by  calculation,  which  requires 
not  only  detailed  infomation  on  the  entire  arc  profile,  but  also  on  the  boundary 
conditions  because  of  the  significance  of  wall-recombination. 

It  should  be  noted,  that  in  this  respect  the  situation  is  not  much  more 
favorable  in  diffusion  flames,  and  the  use  of  equilibrium  relations  in  order  to 
calculate  the  concentration  of  OH  in  these  sources  and  derive  absolute  transition 
probabilities,  has  to  be  viewed  with  caution. 

Much  work  remains  to  be  done  in  chaureieterizing  and  interpreting  departures 
from  equilibrium  in  flames.  Overpopulation  of  excited  electronic  states  can 
occur  not  only  via  chemiluminescence,  but  conceivably  also  by  electron  collisions, 
since  the  concentration  of  free  electrons  in  many  flames  is  much  greater  than  ex¬ 
pected  at  equilibrium,  and  it  has  also  been  suggested  that  the  electron  temper¬ 
ature  is  higher  than  the  gas  temperature.  It  appears  likely  that  observed  in¬ 
consistencies  in  measurements  of  flame  temperatures,  that  in  the  past  have  been 
attributed  to  incorrect  transition  probabilities,  were  adso  caused  by  non¬ 
equilibrium  effects. 
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APPENDIX  -  SIMBOLS 


A 

A(  J,l) 

a 

a(l) 
B(i, j) 
c 

D 

d 

E 

% 

El 

E 

00 

F(i,T) 

F(K) 

Fo 

f(i, j) 

G(v) 

g(l) 

h 

hw 

I 

I(i,j) 

% 

Ia 


area  under  recorded  line  trace  or  atomic  weight 
spontaneous  transition  probability  (s^^) 
coupling  constant 
overpopulation  of  energy  level  i 

ratio  of  effective  to  total  number  of  exciting  collisions 

velocity  of  light  in  vacuum 

dispersion  (mm/A) 

diffusion  length  (cm) 

electric  field  strength  (v/cm) 

ionization  potential  of  hydrogen 

energy  term  pertaining  to  energy  level  i  (cm“^) 

ionization  potential 

Saha  equation  for  atoms  in  energy  level  i 
rotational  energy  term  (cm“^) 
normalized  field  strength 
absorption  oscillator  strength 
vibrational  energy  term  (cm“^) 
statistical  weight  of  energy  level  i 
Planck  constant 
(full)  half  width  of  a  line 
current  ( amps ) 

part  of  excitation  rate  coefficient  that  depends  only  on  the  quantum 
numbers  of  the  two  levels 
integrated  line  radiance  (W  cm“^  sr“^) 
intensity  =  spectral  radiance  (W  cm""^  sr"^  M-”^) 

(principal)  quantum  number 
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J 

j 

K 

k 

k’ 

L 

m 

Nt 

n(i) 

n*(i) 

n*(oo) 

P 

P'(v',v'’) 

p 

Q 

<1 

R 

|R(j,i)P 


total  angular  momentum  q.uantum  number 
(principal)  quantum  number 

quantum  number  of  total  angular  momentum  except  for  snin  or 
equilibrium  constant 
Boltzmann  constant 
net  absorption  coefficient  (cm“^) 
length  of  optical  path  (cm) 

magnetic  quantum  number  or  slope  of  Boltzmann  plot  (cm) 
concentration  of  ions  (cm”^) 

integrated  radiance  of  optically  thin  line  (W  cm"^  sr~^) 
total  concentration  of  particles  (cm"^) 

spectral  radiance  of  optically  thin  line  (W  cra"^  sr~^  ia“l) 

electron  density  (cm^^) 

occupation  of  energy  level  i  (cm“^) 

occupation  of  degenerate  energy  states  belonging  to  level  i 
occupation  of  last  bound  energy  state  (cm~^) 

pressure  (atm)  or  radial  part  of  the  transition  moment  times  a 
constant  (s“l 
Franck-Condon  factor 

partial  pressure  (atm)  or  exponential  correction  factor  (cm) 
internal  partition  function 

specified  group  of  rotational-vibrational  levels 
response  of  detection  system  or  reaction  rate  (cm"^  s“^) 
transition  moment 
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s 


S(K’J',K"J") 

S(a) 

T 

% 

Te 

"g 

“^rot 

U[u(i^  j)] 
u(i, j) 

V 


X 

W 

y 

z 

a(l) 

“a 

€ 

^A 


so-urce  function  (w  cm"^  sr"^  M-~^)  ot  resultant  spin  quantum 
number  or  rate  coefficient  for  collisional  ionization/ 
recombination 

Hoenl-London  factor  (rotational  line  strength) 

Stark  broadening  parameter 
temperature  (K) 

distribution  "temperature"  (K) 
kinetic  temperature  of  electrons  (K) 

kinetic  temperature  of  heavy  particles  (gas  temperature)  (K) 

apparent  rotational  "temperature"  (K) 

function  entering  the  rate  coefficients 

exponent  of  the  Boltzmann  factor 

potential  (volt) 

recorded  signal  (relative  intensity) 

vibrational  quantum  number  or  root  mean  square  velocity 
(cm/s) 

rate  coefficient  for  collisional  excitation/de -excitation 

concentration  of  species  X 

reduced  intensity  of  rotational  line 

total  partition  function 

radiative  recombination  coefficient  (cm^  s“^) 

absorptance 

emissivity 

emission  coefficient  (W  cm“^  sr"^  4“^) 
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^  correction  factor  for  the  intensity  cf  the  continu-um 

A  quantnin  nnmher  of  the  resulting  orbital  angular  moment nm  about  the 

internuclear  axis 

A  wavelength  or  mean  free  path  (cm) 

V  frequency  ( s ) 

Debye  radius  (cm) 
o  cross  section  (cm^) 

SUBSCRIPTS,  SUPERSCRIPTS  MB  INDICES 
c  (ionization)  contin-um 

col  colli  sional 

Eq  equilibrium 

e  electron 

ef  effective 

exc  excitation 

o  line  center 

rad  radiative 

rec  recombination 
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Fig,  1.  Diagram  of  the  rotational  energy  levels  of  the  X^II  and  A^L  states 
of  OH.  The  twelve  branches  corresponding  to  the  allowed  transitions  to  a 
lower  rotational  level  K  >  2  are  shown.  The  branches  are  labelled  ass-uming 
Hand's  case  (b)  for  the  ground  state.  Herzberg' s  notation  for  the  satellite 
lines  is  indicated  in  parentheses. 
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Figure  3.  Calculated  equilibrium  composition  of  a  water  vapor  plasma  at  1  atm  pressure. 
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Fig.  4.  Calculated  equilibrium  composition  of  a  water  vapor  plasma  at  0.02  atm  pressure. 


Ig.  5.  Equilibrium  occupation  of  the  energy  states  corresponding  to  the  upper  energy  levels  for  the 
j_“2  (solid  line),  and  the  Rq-20  (broken  line)  transitions  of  the  (0,0)  band  of  OH  in  a  water -vapor 
lasma  at  1  atm  and  at  0=02  atm  pressure.  An  arrow  indicates  n*(max)  at  the  normal  temperature. 
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Fig.  6.  Photoelectric  record  of  the  hyperfine  structure  of  Hg  I  4358  A, 
scanned  with  the  five  meter  spectrometer  at  a  beam  aperture  of  f/80;  4  mm  slit- 
height,  and  0.030  mm  slitwidths.  The  B&L  replica  grating  (600  lines/mm;  blazed 
at  54°)  was  used  in  the  6th  grating  order  at  a  grating  angle  of  51°40’ . 

Source:  uncooled  low-pressure  mercury  discharge  (germicidal  lamp). 

Detector:  EMI  9558  QA  photomiiltiplier . 


Fig.  7.  The  measiired  chart  positions  of  the  hyperfine  lines  in  the  spectrum  of  Hg  I  4358  &  shown 
in  Fig.  6  have  been  plotted  vs.  the  increment  in  wavelength  units. 


Experimental  set-up  for  concurrent  emission  and  absorption  measurements  of  arc  plasmas 


Fig.  9.  Results  of  the  Ahel  inversion  of  side-on  measurements  of  the  R]_-20  (0,0)  line  of  OH.  The  radial 
distribution  for  the  emission  coefficient  has  been  obtained  by  Barr's  (^)  and  by  Bochasten’ s  (25)  method. 
The  source  was  a  water  vapor  plasma  at  50  Amps,  and  22  Torr.  =  radius  of  arc  channel. 


0009 


i^ig.  10.  Radial  temperature  profile  of  a  water-vapor  arc  at  50  Amps,  and  22  Torr^  as  derived  from  the 
radial  distribution  of  the  emission  coefficients  for  the  Pi-2  and  the  Ri-20  lines  of  the  (0,0)  band  of 
OH.  The  profiles  obtained  from  an  Abel  inversion  of  the  same  input  data  by  Barr’s  and  by  Bockasten's 
method  are  shown  to  the  left  and  the  right,  respectively,  of  the  arc  center. 
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-^ig,  15.  Measirred  occupation  of  energy  states  of  hydrogen  and  oxygen  in  the 
center  (r  =^0),  and  at  the  periphery  (r  =  1.1  mm)  of  an  arc  at  atmospheric 
pressure  and  4.1  Amps.  The  sample  gas  was  argon  saturated  with  water  vapor. 


134 


6400K 


Fig.  14.  Occupation  of  the  hydrogen  energy  states  with  principal  quantum  numbers  2  to  6;  measured  in  the 
center  of  an  arc  in  water  vapor  at  22  Torr  and  30  Amps.  is  the  distribution  temperature  for  neigh¬ 

boring  energy  levels.  The  equilibrium  distribution  is  indicated  qualitatively  for  an  electron  temper¬ 
ature  of  25,000  K. 


PARENTAGE 


Fig.  15.  Occupation  of  oxygen  energy  states  in  the  center  of  an  arc  in  water  vapor  at  22  Torr  and  30  Amps 
The  predicted  population  distribution  over  the  higher  states  corresponding  to  singlet,  triplet,  and 
quintet  energy  terms  is  indicated  qualitatively  by  broken  line. 
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Fig.  16.  Photoelectric  record  of  a  scan  in  the  9th  grating  order  of  the 
Q  e(0,0)  line  of  OH  at  3087  X  and  of  its  satellite,  emitted  in  the  end-on 
mode  by  an  arc  in  water  vapor  at  50  Amps,  and  20  Torr.  The  spectrometer  was 
focused  on  an  arc  region  1.1  mm  off  the  arc  center.  Slit  widths:  0.025  mm* 
slit  height:  0.25  ram;  aperture:  f/120. 

AXi  is  the  spectral  slit  width;  hw  the  half  width  of  the  line  corrected  for 
self-absorption.  The  calculated  profile  for  a  Gaussian  of  half  width  hw'  and 
normalized  at  the  line  peak,  is  indicated  by  dots. 
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that  for  the  R(0,0)  branch,  shown  in  Fig.  17. 
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Fig.  21.  "Apparent  Boltzmann  plot"  for  the  Q(l,l)  branch  in  run  C-1.  The  plot  drawn  is  identical 
to  that  for  the  R(l,l)  hranch,  shown  in  Fig.  20. 


'Apparent  Boltzmann  plot*'  for  the  R(1,0)  branch  in 
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Fig.  23.  "Apparent  Boltzmann  plot"  for  the  Q(l,0)  branch  i 


Fig.  24,  "Apparent  Boltzmann  plot"  for  the  P(l,0)  branch  in 
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Fig.  25.  "Apparent  Boltzmann  plot"  for  the  P(2,i; 
splitting  of  the  Boltzmann  plots  for  lines  with  an 
and  an  upper  spin  state  (cross)  is  shown. 
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